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PEEFACE. 


My time for some years was so muoh occupied by 
administrative duties in Trinity College that I was unable 
to make any attempt to carry out the intentions stated 
above fourteen years ago in the Preface to the first volume 
of this treatise. 

I have now to some extent accomplished what I then 
proposed. I came, however, to the conclusion that a 
ohapter on Conjugate Functions was not suited for such 
a treatise as the present, and that to a student having a 
limited amount of time at his disposal some account of 
Maxwell's Theory of Light would be more interesting and 
instructive. This theory is not of course part of the 
Theory of Attraction, but is so intimately connected with 
bhe properties of magnetized bodies, electric currents, and 
lielectrics treated of in the present volume that its 
introduction does not seem unsuitable. 

1 should recommend a student reading this book for 
the first time to omit the whole of Chapter Ylll after 
Article 146. 
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Preface, 


Of the more recent developments of the electro- 
magnetic theory of light I have not attempted to give 
any account. So far as I can judge some of these rest 
on insecure foundations. I imagine, however, that before 
studying the most recent investigations a preliminary 
knowledge of Maxwell’s theory is required, and I trust, 
therefore, that my chapter on the subject will not be 
entirely useless to the student. 

I have to thank Mr. S. B. Kelleher, f.t.c.d., for his 
kindness in reading the proof-sheets of this book, and 
furnishing me with many valuable corrections. 


FllANClS A. TAKLETON^. 


Tkinity Colleqb, Dublin. 
Aprils 1913. 
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CHAPTER VIIL 


SPHERICAL AND ELLIPSOIDAL HARMONICS. 


Seohon I . — Spherical Surfaces, 


137. Expansion of Potential in Herles of Solid 
Harmonics. — It was shown in Art. 78 that the potential V 
at a point P, more distant than iany point in the attracting 
mass from the origin, can be ! expanded in a series of 
descending powers of r, where r denotes the distance of P 
from the origin. 

In this case, the series for the potential is of the form 


^ + T 


where M denotes the attracting mass, and Pi, ^ 2 , &c., are 
functions of 0 and 0, the singular coordinates of P, and of 
constants depending on tjh:6 attracting mass, but independent 
of the position of P. 

Since V* F" = 0 for all positions of P outside Af, the 
coefficient of each power of r in V* F must vanish separately, 
and therefore 

Usiijg for V* the expression given, equation (17), Art. 48, 
we obtain 

^ /I ^ ^ Fr, . ' t \ Tr A /i\ 

_ (1 - -jp- + • (» + 1) I . . 0. (1) 
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Spherical and ElUpeoidal Harmonics, 


If JTbe a function of d, and 0 satisfying (1), it is easily 
seen that V*(r”F)=0; accordingly, if we have 

V* Vn = 0, and Vn f homogeneous function of y, z of the 
degree n satisfying Ld). ^ace^s equation. Such a function is 
called a spherical solid K rmonic of the degree n. 

It appears from whay has been said that if denote a 
solid harmonic of the degree w, then Vn is &o a solid 

harmonic whose degree is - (;^ + 1). 

The function is termed a spherical surface harmonic of 

the degree w, and is what has been denoted above by F„. 

In the present case, hy considering the expression from 

whose expansion was obtained, it is easy to see that Vn 

is a rational and integral function of x^ y, s. In what follows, 
Vn will be termed a solid, and ¥» a spherical Ijarmonio. 

138. liaplaee's and liegendre^s Coefficients. — If 
the attracting mass be concentrated at a point Q whose polar 
coordinates are /, O', <p\ and whose distance from any point P 
is r, we have 


r r 


r r 
1-2A- + 






where X = fi// + 1 - fi‘ - ju'* cos - f'). 

In this case, if P be farther than Q from the origin, 

/n 

V = if2X„ — ; 
and if P be nearer the origin, 

F=J/SX„:^. 


The eoefficieuts Xi, Lt, &o., in the development of r'* are 
called Laplace's Coefficients. They are obviously spherical 
harmonics of a particular kind. They may be denied as 
the e^ffifdents of the successive powers of h in the expanmm if 
(1 - 2\A + A*)-». 

These coefficients are plainly symmetrical with respect to 
the angular coordinates of P and Q. 



Properties of Complete Spherical Harmonics, 3 

If the point Q be on the axis froia which 0 is counted, 
X = )Le, and Laplace^s coefficients beooi^e tAe coefficients of the 
successive powers of h in the development of 

(1 - 2ph -f h^)f. 

In this case, these ooeflSioients 4re functions of p solely, 
and are called Legendre* s Coefficients, They are usually 
denoted by Pi, P 2 , &c. ) 

It is plain that P« satisfies the{' equation 

( 2 ) 

In general, a spherical harmonic of the degree w, whicli 
is a function of p solely, satisfies (2), and is called a zonal 
harmonic. 

139. Properties of Complete filpberical Har- 
monics. — A spherical harmonic which when expressed as 
a function of the coordinates is finite and single-valued for all 
points of space, is said to be complete. If and Yn be 
complete spherical harmonics of diflfereiit degrees, 

|'|j"'r„r„rf/urf0 = o. (3) 


This may be proved as 
137 that 




follows: — It appears from Art. 


uiul 


^,n+l 


satisfy Laplace’s equation ; and by (6), Art. 58, if we take 
ns the field of integration the space outside a sphere 8 of 
radius a described round the origin as centre, we have 



Hence 

nnA fViArAfnrA. iitiIarr m s n. eouation (3^ must hold irood. 
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Spherical and Ellipsoidal Harmonics. 


If Yn be a complete spherical harmonic of the degree w, 
and Ln a Laplacian coefficient of the same degree, 


r+i rzir 

J J L^Yn d/j. 




2n + 1 


- Y' 

t n* 


(4) 


To prove this, take as the field of integration the space 
outside a sphere 8 whose centre is at the origin, and whose 
radius a is less than r\ the distance of tlie point Q from the 
origin, and let r denote the distance of any point from Q ; 
Tn 

by (10), Art. 69, we have 


then the function satisfies Laplace’s equation, and therefore 


but 


J di 


r J r'-*' 


« 4- 1 F„ 


i . 2£. ’ 


/n +1 


at all points for which r < r', whence at the surface 5 we 
have 


iL 

dr 


- ) = ■ 


,.n-l 

Tn+V 


also dS = a- dp d(p^ and 

TjfnYndpd<p = 0 , 

unless m = n. Hence we obtain 


1 . 




from which equation (4) follows by transposition. 

^ If two series of spherical harmonics are equal for aU 
values of p and each harmonic of one series is equal to 
the harmonic of the other series whose degree is the same. 
Here 

Yq + Yi + jP^ + <&o. ■■ Eo + El + E2 + &o. 
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If each side of this equation be multiplied by and 
integrated, since 

I Ln F„, dfi = 0, 

by (4) we obtain 

47r 47r ^ 

c\ r ^ n — ,T " n y 

2n -H 1 Zn + 1 

and as this equation holds good for all values of p and (p\ 
we get Yn = Zn. 

Any function of u and 0 which is finite and single- valued 
can be expanded in a series of spherical harmonics. 

The method of arriving at this result is suggested by 
what has been already proved. If it be possible to express 
in the form SF«, we must have 

ATrf{pff^') = 47r2 Y'n « 2 (2n + 1) [ [ Ln Yn dp dtp 

J-iJo 

= 1 /(a(^)S(2m+ (5) 

J-lJ 0 


Now ( 1 - 2 X/i + A*)-i = SZ„ A" ; 


■whence, differentiating and multiplying by 2//, we have 


2(X/i- A’) 

(I - 2AA + A’)^ 


S2»Z„A”, 


then by addition to the former equation we get 

1 - 


(1 - 2XA + A’)t 


= S(2» + ])Z„A”. 


Aoijordingly, if the supposed expansion be possible, we must, 
when A = 1, have 

(/9') = [_^ £ / (1— JaaWi 


and conversely, if this equation be true, the expansion is 
possible. 



6 Spherical and Ellipsoidal Harmonics. 

That equation (h)'i8 true can be shown in the following 
manner : — 

Let Q be a point outside a sphere S whose centre is the 
origin and whose radius is and let r denote the distance of 
any point on S from Q. Then 

r* = ^ - 2\ar\ 


where r' denotes the distance of Q from the origin ; and if 

we have 
r 

1 « /;2 (,/2 ^ ^ 2 ) 
r^ 

As in Art. 42, we have 


and therefore 

J r" 



1 V M 

r Vr' -a r' a) r\r^ ~ «*) ’ 


also dS^d^dpd^j and accordingly 

r+i C2. I _ ^2 


The value of the definite integral above is therefore 
independent of r' ; but A = X when r' *= a, and in this case 
each element of the integral in (6) is zero, unless r be 
infinitely small, in which case p « and ^ Hence, 

when h « 1, we have 

-/("'rtlld _ 2 a'/vA’)*'''"'* 



Application of Spherical Harmonics. 7 

140. Application of Spherical Harmoiiics. — When 
the potential is due to mass on one side of a spherical surface 8 
and is given at each point of the surface S itself, the potential 
at any point on the side of 8 remote from the mass can be 
represented by a series of solid harmonics. At the surface 8 
this series becomes a series of spherical harmonics representing 
the known value of the potential at the surface. Hence by 
Art. 139 each harmonic in this series is determined, and con- 
sequently so also are the corresponding solid harmonics 
representing the potential on one side of 8, 

If the potential be due to a distribution of mass on the 
surface 8 whose density is given, the potential outside 8 can 
be represented by the series 

y ^ ^+1 ’ 


and at any point inside by the series 



At all points of the surface these two expressions must 
be equal, whence by Art. 139, Zn = Tn- Again, if V and V' 
denote the potentials outside and inside the surface, we have 
by Art. 46 at the surface 


that is. 


whence 


dr 


dV' 

dr 


+ Airtf = 0 , 


M Y 

~ + 'S (2« + 1) — = 47r<r = An'S.Sn ; 
' 'a* 


M = 47ra’/S'a, 


Y = S 
" 2» + l 


TT" i 4 ^ ^ 

r=4w«*2^, 


r'=4wa»2 - 




and 
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141. Potential of Homogeneous Spheroid. — the 

surface of a solid dijBfers but little from a sphere whose centre 
is at the origin, the radius vector r is given by an equation 
of the form r = a{l + ay), where « denotes the radius of the 
sphere, y a function of the angular coordinates p and and a 
a small constant whose square may be neglected. 

The potential V at any external point is the sum of the 
potential due to the sphere and of that due to the shell whose 
thickness at any point is aay. Hence if p denote the density 
of the spheroid, /x', 0' the coordinates of a point on the surface 
of the sphere, and r the distance of this point from the point 
r, ju, 0 in external space, we have 


4 TTpa^ 

r = H + « 
3 r 


' pa^l/ dp 

r 


but by Art. 139, y and therefore by Art. 138, and 

by (3) and (4), we get 


4^ 

3 r 

4 Ttpa^ 
" 3“ 


+ apa® 


\L„Y'„a”dfii 


+ ^airpa^ 2 


«”F„ 




(2«+l)r"+' J 


( 7 ) 


For the potential at a point inside the sphere, by Art. 42, 
we obtain, in like manner, 

V = 2*rp«* - f Trpr* + 4an-p«’ 2 

142. Potential of Heteroseneons Spheroid. — If 

a spheroid be composed of homogeneous layers comprised 
between surfaces given by equations of the form 


r = a(l + aXY„), 

where is a spherical harmonio which varies with the 
surface, and a is a variable parameter, we have, for the 



Potential of Heterogeneom Spheroid. 9 


potential SV of a single layer at a point outside, the 
equation 


Sr=4irP-^ + ia7rp'2 




and at a point inside, 

SF^ 4vpa da + 4a;rp 2 


Hence for the potential F" of a heterogeneous spheroid at 
a point outside it, if ai denote the parameter of the external 
surface, we obtain the equation 




rai 

47r pa^ da 
Jo 


+ iair 2 


pdia-^^Y,) 

J_o 

(2n 4. 1 );•«+! 


(9) 


For the potential of a heterogeneous shell comprised 
between surfaces whose parameters are ai and at an 
internal point, we get 


r«i [ 

F = 4ir pada + 4air 2 rv— ' 

J (^»» + 1) 


( 10 ) 


By combining the expressions given by (9) and (10), we 
find for the potential of a heterogeneous spheroid, at an 
internal point lying on a surface whose parameter is a, 
the equation 


4„ ra 

F = - (1 - aSF„)[ pa’rfa + W 

... 

+ 47r pada + 4o7r 2 — Tin 

J« ^ (2m+1) 


( 11 ) 
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143. Homogeneous Hass of Revolving Fluid. — 

If a homogeneous mass of fluid revolving with a uniform 
angular velocity be in a state of relative equilibrium under 
its own attraction, its external surface, if it be nearly 
spherical, must be an ellipsoid of revolution. 

This may be proved as follows: — 

By Ex. 5, Art. 24, at the free surface of a liquid in 
relative equilibrium, if V denote the attraction potential, 
which in this case is a force function, and if the axis of 
rotation be taken as the axis of 2 , we have 

dV + (i)^{xdx + ydy) = 0 ; 

whence, as in Art. 81, we get 

V + (1 -y?) = constant. (12) 

The last term on the left-hand side of this equation must 
be small, as otherwise the surface of the liquid could not be 
approximately spherical. In this term, therefore, we may 
put r = fl, and substituting for V from (7), we get 

- aSFn) + 4airpa^ 2 ^ + ^^(1 - /u*) = constant, 

o Zn + 1 Z 

(13) 

In order to make use of this equation, we must express 

/tt* by means of spherical harmonics. Since /u* = jj, it is 

plain that the solid harmonic corresponding to the spherical 
harmonic of highest degree in n* must he z* + ki^, where k is 
an undetermined constant. To determine k, we have 

V* {s’ + A: (a:* + y* + s’) ) = 0. 

Hence A: = - J, and /a’ - ^ + i is the required expression 
for /i*. 



Figure of the Earth. 
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By Art. 139, the sum of the spherical harmonics of each 
degree above zero in (13) must vanish separately. Hence 
Frt = 0 if n> 2, and 


8 

Lo 


airpci^Yi - 




Putting ^ = q, we get « (1 - 


Hence the equation of the free surface is of the form 


r = «{1 


which represents an ellipsoid of revolution nearly spherical 
whose ellipticity is i q. See Art. 81. 

144. Figure of the Earth.— On the hypotheses that 
the Earth is composed of homogeneous layers bounded by 
similar surfaces nearly spherical, and that it is covered with 
liquid in relative equilibrium, it is easy to show that the 
external surface of the liquid must be an oblate ellipsoid of 
revolution wliose axis is the axis of rotation. 

The attraction potential V of the Earth is given by (11). 
At the surface of the liquid, (12) must hold good. Hence, 
by substitution, we obtain 


47r 

(h 


(1 - a2 





pa^ da + 4a7r 2 ■■ 


( 2 // + 


+ (J - fi^) = constant. (14) 


Since the surfaces of equal density are similar, Yn does 
not vary with a, and as Oi is the greatest possible value for a 
if n be not less than 2, we have 

, (n + 3) «"« 

(2» + !)«,« 

Hence the multiplier of Y„ in (14) cannot be zero, and 
therefore if « > 2, we have Y„ = 0. 
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If n - 2, we obtain 

47r J pa^da - j pa^da\^ aYi = • (15) 

By Art. 78, when the centre of inertia is the origin, the 

coefficient oi — in V is zero. Hence, in the expression for 

the potential of a spheroid given by (9), if the surfaces of 
equal density be similar, and if the centre of inertia be the 
origin, we must have Zi = 0 ; and in the present case the 
form of the external surface is determined by the equation 
r = (1 + a ri), where Pi is given by (1 6). The external surface 

is therefore an oblate ellipsoid of revolution having the axis 
of rotation as its axis. 

It seems improbable that the hypothesis made above with 
respect to the form of the surfaces of equal density should be 
correct. In order that it should be true, it is necessary that 
these surfaces should have been formed under similar con- 
ditions ; but, unless the Earth were of uniform density, this 
could not have been the case, since the equatorial centrifugal 
force due to rotation varies as the distance from the centre, 
whilst the attraction of the sphere having this distance as 
radius varies in a different manner unless the density be 
uniform. 

A more probable hypothesis is, that the surfaces of equal 
density are represented by equations of the form 

r = a(l + ahy)y 

where A is a parameter varying with a, but constant for each 
surface, and y a function of p and 0, which is the same for all 
the surfaces. 

145. Clalraut's Tbeorem. — Whatever be the internal 
constitution of the Earth, if it be covered with liquid in 
relative equilibrium whose external surface is an ellipsoid of 
revolution nearly spherical, the elliptioity, 6, the ratio of the 
centrifugal force at the equator to gravity, q, and the diffe- 
rence between polar and equatorial gravity divided by the 
latter, y, fulfil the relation y + c «= -ly. 



Clairaut's Theorem. 
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This equation was proved in Art. 81 on a particular 
hypothesis as to the internal constitution of the Earth. 
Any hypothesis of this kind is, however, unnecessary, as 
was first pointed out by Sir Or. Stokes. 

At the external surface of the liquid, the Earth’s poten- 
tial V must satisfy (12); but as this surface is nearly spherical 
and the term in (12) due to rotation is small, the variable 
terms in V must be small. Hence, if M denote the mass of 
the Earth, we may assume 


-f a 2 “ 
r 




where a is a small constant. Again, by Art. 81, the form of 
the external surface is represented by the equation 


Substituting in (12), we get 

— { 1 + € (/t* - J)J + « 2 ~ ^ (J ~ m') “ constant. 

U CL 


Hence = 0, unless n = 2. If n = 2, we have 

where 

q = ttra / — • 

^ / a* 

Accordingly, 

„ M a^Yt M Ma^fq \f 1\ 

— = (16) 


IE" 

dr 



3a 


r* 


and 
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Hence, if g denote the acceleration due to gravity at any 
point on the Earth’s surface, in the same manner as in 
Art. 81, we find 



- (1 - ; 


that is, 


g 





Hence 


and therefore 


Qp - (/e 

7 == 

Qe 



7 + 


2 


(17) 


146. Tangential Comitonent of Attraction. — If P 

denote the component of the Earth’s attraction perpendicular 
to the radius at any point on its surface, by (16), we have 


P 


dF 

rdB' 


— T (S' 


dn 

de 



sin 2X, 


(18) 


where X denotes the latitude of the place. 

If we compare (16) with (2), Art. 78, we get 

ilf«’(2£ = iil -(A + B+C). 


Hence 


C~A=-^(2e-q). 


(19) 


The equations proved above were arrived at before in 
Art. 81 by means of a special hypothesis with respect to the 
internal constitution of the Earth. The facility with which 
these results have been obtained in the present and preceding 
Articles without any such hypothesis illustrates the power of 
the Laplacian method. 



Legendres Coefficients. 
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147. Itegendre’s Coefflclents. — The definition of these 
ooe£Soients given in Art. 138 enables us to see that they are 
rational and integral functions of ju. A general expression 
for these ooefiSoients cannot be readily obtained by the usual 
methods of expansion. If we put 

fJz 

(l-2yaf + y*)-i = — , 

we can, by integration, get rid of the negative index ; and 
thus we obtain 

(1 - 'iyx + w*)4 

s = - — + constant. 

y 

If we take i for the constant, we get 

y 

=l-2ya: + y*; 


whence we obtain 


s = a! + I (s’ - 1). 


( 20 ) 


We have now an expression for s suitable for the applica* 
tion of Lagrange’s theorem (Williamson, Differential Cakulm, 
Art. 125) by which we obtain 

i»*-l y^ dfc^-iy , 

z-x + y 2 +1-2 2 l'2-3 5 

whence 




d fx'^- r 

+y — I -?r- I + 


Hence, if 
we get 


dx\ 2 

(1 - 2fih + A»)-4 = 1 + 2P„A’*, 


f (x^- 1\2 
l-2rfa!»\ 2 j 


D„ = 


d» 


1-2-3 . . . n-2” dfi 




( 21 ) 


The development of P„ in powers of n is most easily 
effected by means of the differential equation satisfied by 
zonal harmonics. 
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If 8„ denote a zonal harmonic of the degree n, we may 
assume 

Sn = a, p’ + a,_i ju*'* + &o., 

where s, &o., must be positive in order that Sn should be 
finite at every point of space, and Sn must satisfy the 
equation 

+ + = ( 22 ) 

Hence we get 

s{s - l)ff, ju*"* - «(« + l)fl, p' - (s- 2) (s - l)a,.t + &o. 

+ n{n + 1) [Uap’ + «Ste.} = 0. 

Here the coefficient of p‘ is 

{«(« + !) -s(s + l)j«,, 


and that of /u’"* is 

«(s - l)a, + {«(» + 1) - (s - 2)(s - 1)} a,.j ; 


accordingly, as each coefficient must vanish separately, we 
obtain 


(« - «) (s + n + 1) = 0, 


*-*“(«-s + 2)(« + 8-1) “ 


From the first of these we get s = n, or « = -(« + !); 
and as the negative value for s is here inadmissible, we have 

- « (w - 1) - (« - 2) (» - 3) 

8 = «, «„.2 - _ 1) - 4(2» _ 3) 


and in general 


„ / ■,,^ «(«-l)(«-2)(n-3)...(«-2g+2)(«-2g + l) 

2.4...2g,(2n-l)(2»-3)...(2n-2j + l) 

( 23 ) 
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It is plain that the terms in (22) resulting from 
/u*"’ + aa-3n‘'^ + &c., 

must vanish independently of those arising from the series 
already considered, and that we get for the first term the 
equation 

(s - n - 1) (s + n) a,.i = 0, 

Hence &c., must each be zero, and we obtain 


“r 2(2«-l)'^ 2.4.(2//-l)(2«-3)'^ ® y 

(24) 

Hence, zonal harmonics of the same order can differ only 
in the constant factor, and we may write 

S„ = aP„, (25) 

where a is an undetermined constant. 

It is easy to see that the coefficient of ju” in (21) is 


2M(2n- 1)(2»-2) . . . (h + 1) 2" . [ « . 1 . 3 . 5 , . . (2» - 1) 

— ; or ^ — - 

2 " 


2 “ . 


and therefore 

1.3.5.(2«-1), „ . , 

j-^ 

It follows from the definition of P„ that when /* = 1 the 

voliiA P„ ifl nnifv- 
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148. Spherical Harmonics. — Since the expansion of 




oontaius only rational and integral functions of Xj and 
the coefficients &o.y must be rational and integral 

functions of sin^ and cos0, in which each power of sin0 
and cos^ is multiplied by the same power of -v/(l -- /u*). 
Hence, as 

47r 

the spherical harmonic must be a rational and integral 
function of sin 0 and cos ^ of the degree in which each 
power of sin 0 and cos ^ is multiplied by the same power of 
v^(l - fx^). If each power of sin 0 and cos ^ be expanded in 
a series of sines and cosines of multiples of fj>, we see that 
finally Yn is reducible to the form 

S {As Mg cos + BsNs sin 

where As and Bs are undetermined constants, and Ms and Ns 
functions of y. 

If we put — = 2>, and /it* - 1 = equation (1) becomes 
1 (PY 

Du D Yn + - ~ n{n + 1) F,, - 0. (27) 

Since the coefficient of the sine or cosine of each multiple 
of 0 must vanish separately in (27), we have 

EuDM, - + «(« + l)j M, = 0. (28) 

Again, since cos 80 and sin 80 can result only from 
(cos0)% (sin0)*, (cos0)*+*, (8in0)*^’, &o., 

M, must contain w’ as a factor ; and the other factor must be 
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a rational and integral function of p. Accordingly, we may 

put „ i 

Ma =* , 

where v denotes a rational and integral function of p. 
From (28), we have then 

5 s a 

V + sti^r + ( 2 .<? + 2)ti^ fiDv + 

'1 ^ 

- H^ir V - n (n + 1) r = 0 . (29) 

a 

Since /< = /Lt* - 1, equation (29) is divisible by and we 
get 

hD^v + (s + 1) Du Dv + - n (n + 1) t? = 0. (30) 

If we assume v = D^w^ equation (30) becomes 

(uDw) - n(n + l)Hw = 0. (31) 

Since t? is a rational and integral function of ju, it is plain 
that, with the exception of a constant factor, it is completely 
determined by (30). Hence any rational and integral function 
of fi wliich satisfies (30) or (31) must represent v. Equation (31) 
is satisfied if w satisfy 

Du Die - it {ti + 1 ) ?«? = 0 ; 

but this equation is the same as (22). 

Hence we may put w = P„, and we l»ave 

V = Ms = u'^D^Pn, 

It is plain that the equations by which iV, is determined are 
the same as those for Mg. Accordingly, these two functions 
can differ only by a constant factor, and we obtain 

£ 

Y„ = S (A, cos + B, sin »^) 

C2 


( 32 ) 



20 Spherical and Ellipsoidal Harmonics. 

The part of F„ depending upon that is, 

s 

(A, cos + B, sin «^) 

is termed a tesseral harmonic of degree n and order s, and we 
may write 

F„ = S T„s (As oos sijt + B, sin S(j>). (33) 

If we substitute for P„ its value given by (26), since As 
and Bs are undetermined constants, we have 


r...=(iV)* 


X 



(«-s)( w-s-l ) 
2(2»-l) ^ 


}(n-s-2)(n-s-3) 

2.4.(2»-l)(2«-3) 


-&o.j. (34) 

149. liaplaee^s Coefficients. — Laplace’s coefficients 
are, as we have seen in Art. 138, a particular kind of 
spherical harmonic ; and as they are functions of 0 and 0' 
through being functions of cos - 0'), and are symmetrical 
in p and p\ we must have 

Ln *= cos S (0 - <j/) (35) 

where Os is a definite function of n and s, which may be 
determined in the following manner; — 

By (4) we have 

*2»r / ~ 1 \ 

|5««(coss^ cos«^' + sin 5^ sins^Ow^w'^D'PwZ^'^P'nj 

( i 

X I %{A, cos + Bs sin si>)u^D^Pn 

2 ^^ ^ £ S (^As cos 50^ T Bs sin 50^) B'^P'^, 



^ dp dip ^ 
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It is plain that the only part of the multiplier of oos 
in the left-hand member of this equation which does not 
vanish after integration is 

a^As D'^P'n I U^jyPnY oos* 50 (Ifi d(^. 


Hence, if s be not zero, we get 

= (36) 

and if s = 0, we have 

The first term in Ln is a^P^P^ny when a^'= 1, all the 
other terms vanish, and P'n = l. Hence, in this case, 
Ln - a^Pn ; but Ln becomes P„ when fi = 1, and therefore 
r^o = 1. The remaining coefficients can now be found by 
means of (36) and (37). 


Let 


u^iL^PuYdfi = A.,, 

-1 


then Ai = J u(l)P„ydix ; 

also by (2), we have 

u (« + 1) Ao = / F„D{uDP„) dfi ; 


whence 


Ai + n(n+ 1)A„ = | {P«D{mDP„) + «PP„PP„} d/x 
= I ^ D{uP„DP„)dfi = 0, 


since u Vlanishes at each limit of the intefiral. 
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It is now easy to see that an equation similar to that 
-obtained above holds good for any two suooessive integrals 
of the series. In fact, hy (2), we have 

iy*\uDP„) - «(« + i)iyPn = 0 ; 

whence, remembering that = 2, D^n >= 0, by Leibnitz’s 
theorem, we have 

uIP*^P„ + (s + l)DulJ^*'P„ + s{s+l)I)‘P„ = n(n + l)D‘P„, 
and therefore 


(« - s) (« + s + 1 ) . (38) 

Hence we have 
A,+i + (n - s) (n + s + 1) A, 


-J 

-1 



P(u‘’*^l)^P„lP*^P„)dfi = 0. 

-1 

(39) 

Accordingly, A,+i « - (n - 8)(u + s f 1) A, ; 
and therefore, by (36), if s be not zero, we have 

(40) 


(1^ 

(n-s){n + .S + 1)’ 

(41) 

and, by (37), we get 



~ «(« + !)’ 

(42) 


Hence, as </o = 1, we obtain 

-2 

«(«+!)’"'’ «(« + ])...(» + s)(« -1)(« -2). . .(ji -s + 1)’ 


, . 2(- 1 )*!«-. s- 

that 18, ff, = — - — L= ( 43 ) 

I M + « 

and 

i„ = P,F„+25(-l)‘(l-M7(lV*)’i>‘-P»0'*i’'.. co8s(^-f). 


( 44 ) 
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150. Complete Harmonics. — The definition of solid 
and spherical harmonics in general has been g^ven in 
Art. 187 ; but the properties of spherical harmonics proved 
in Art. 189 have been obtained on the hypothesis that these 
functions are finite and single- valued for every point of 
space, and in that Art. m and n denote integers. 

If Yi denote a function of p and ^ satisfying the equation 


d dYi 1 d^Yi 

dp^ dp 1 - m’ 


+ i{i+ i ) t ;- = 0 , 


(45) 


where i denotes any real numerical quantity, corresponding 
to Yiy there are two solid harmonics, viz. r* Yi and Yi. 

As i is real, one of the quantities ^ and - {i + 1) is 
negative. 

Hence, selecting the two solid harmonics of negative 
degree which correspond to Yi and Yj, we see that when Yi 
and Yj are finite and single-valued, equation (3) holds good, 
unless « or i = - (/ + 1). Again, if Yi be finite and 
single- valued, by a process similar to that employed in 
proving (4), Art. 139, we have 

f+i r2ir / a / a 

J Jo = , (46) 


if i be positive, the coefficient of X„ being n - i if i be 
negative, and that of Y'l being 47r 


( 0 - 


Accordingly, by (3), F'* = 0, unless i = or i = + 

In either case Yi is a rational and integral function of )u, 

v/l - cos and v^l - sin ^ 


of the degree n. 

Hence we conclude that the degree of a complete spherical 
harmonic must be a positive integer j and that the correspond- 
ing solid harmonic of positive degree must be a rational and 
inte^al function of y, and z. 

This last result is usually expressed by saying that every 
complete solid harmonic is a rational and integral function 
of Xy y, and «, or can be made so by multiplying by a suitable 
power of r. 
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151. Reduction of a Function to Spherical 
Harmonics. — It was shown in Art. 139 that a finite and 
single-valued function of p and ^ can always he expressed 
by a series of complete spherical harmonics. If this series 
be finite, so that 

Fn, 

we have 

,, I ^ ,.2 + &c. ) . 


Hence r^f = /w + rfn^i where /» and /n-i denote rational, 
integral, homogeneous functions of the coordinates ; and it 
appears that if a function of p and 0 can be expressed by 
a finite series of spherical harmonics, the corresponding 
function of the coordinates must consist of a rational, 
integral, homogeneous function, together with another such 
function multiplied by r. Accordingly, the problem to 
express a given function of p and ^ in a finite series of 
spherical harmonics, when soluble, is reduced to that of 
expressing a rational, integral, homogeneous function of 
the coordinates in a series of the form Vn + Vn ~2 + 

This is effected most easily by means of Laplace’s 
operator. 

In fact, by Leibnitz’s theorem. 


Vm = + 2prP'^ 


X 


dr dV,n 
dx (fx 


^drdVm drdVm\ 

dy dy dz dz ) 




(47) 


but 


rf* d 

=P-T 
dor dx 


= + {p -2)1^1^], 


and W =!>(/>+ ; aUo, V*F‘m = 0, and 

Accordingly, 


dVm dVm ^ dVm 

dx dy dz 


WFm = {j»(i» + l) + 2pm\t-f'*Vm 


(48) 
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From (48), we get 
= ^2 Vn ^2 + a^r^ Vn-i + &c., 
&e., 


&C., 

&c., 


where Uzy « 4 > &c., &q,, are known numerical coefficients. 

If this process be repeated sufficiently often, we find 
ultimately 

^^%-hV,y or 

according as n = 2qy or = 2g + 1, the coefficients h and k 
being known numbers. In fact. 


h = + 1, and k = 


By the equations previously obtained, we can then determine 
the other solid harmonics. 

As a simple example, let 


/ = + Bif + Cz^ + 2Fyz + 20zx + 2Hxify 


then / = F 2 + r* and V*/ = h Fq, where A = 3 . 2 . 1 = 6 ; 


but V V’ = 2 (-4 + i? + C) ; hence Fo = J (-4 + JS + C)y 
V 2 = Ax^ + + Cz"^ + 2F^z + 20zx + 2Exij 

-\{A^B^C) {x^ + + s*), 

and 

/= ^ {2^ - J? - C)x^ + {2B - C-A)f + {2C-A- B)z^\ 
+ 2Fyz + 2Gzx + 2Hxy + {^A + J5 + O) r^. 


Again, let / = x^ + + x^y + y'^Xy 

then f ^ F 3 + r^ Fi, and V*/ = k Fi, where A* = | . 3 . 2 = 10 ; 
accordingly, 10 Fi = ^{x + y), hence Fi = |^ (ir + y), and 
Fa = i (x^ + // + x^y + y^x) - ^ z\x + y). 

The method originally given by Laplace for reducing to 
a series of spherical harmonics a function of p and ^ corre- 
sponding to a rational and integral function of the coordinates, 
differs somewhat from that given above. 
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A rational and integral function of the coordinates 
corresponds to a rational and integral function of 

/I, \/ 1 — ju* cos {py and y/l - ju* sin 0. 

If the various powers of cos tp and sin tp be developed in sines 
and cosines of multiples of 0, the series multiplying 

f 

(1 - /U*)* OOSStp 

will contain all the powers of not exceeding n - s, where n 
is the degree of the given function of the coordinates. 

If we collect together the terms containing the highest 
power of jj. in each series, we obtain an expression of the 
form 8_ 

2(^5 QOSS(p + Bs sin S(p) (1 - ju"’* + 

the function Yn may then be determined by taking its 2n + 1 
arbitrary constants, so that the terms of the above form may 
be equal to those in the expression given above. If we 
subtract Yn thus determined from /, we get a function, 
f - Yn of the degree n - 1 in 

p, - ju* cos <py and sin <p. 

The harmonic can then be determined in a way similar 
to that employed in finding F,*, and so on. 

When the original function of the coordinates is trans- 
formed into a function of r, /a, and <p, the various powers of r 
are in / regarded as constants. 

It is plain that the total number of terms or of independent 
constants in / is 1 + 3 + 6 . . . + + 1, that is, {n + 1)*. 

This is also the number of arbitrary constants in the series 

Fq + Fi . . . + Fn. 


152. Hietliods of forming Complete Siolid 
monies. — A complete solid harmonic of positive degree is, 
as we have seen. Art. 150, a rational and integral function 
of the coordinates. A solid harmonic of the degree n, since 

it is homogeneous, contains, therefore, ^ terms. 

2 
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The coefficients of these terms are, however, not all inde- 
pendent ; for, if Vn denote the harmonic, V*Fn must vanish 

for all values of the coordinates, and therefore — 

equations must be satisfied by the coefficients of Vn- 

Accordingly, Vn contains 2n 4 1 independent arbitrary 
constants. 


Since VM - ) = 0, we have 


\dgj \diS J r 


where /, y, k denote any integers. 


Hence 


and therefore 




/ dyfd\JfdY 1 


^ - = 0 , 


dxj \dt/J \dz) r 

d Y / y Y Y ^ 

dxj \dyj \dz) r 


is a solid harmonic of the degree -(e 4/4^4 1). If 
t-f,/ 4 A = Hy the number of different combinations of the 

which can be formed is ; but all the different 

^ 1 

functions which result by the use of these operators on - 

are not independent. In fact, ^ 

Y [AY [AY 1 = 0 

\dx^ ^ dy^^ ^ dz^j \(hj \dy) \dz) r ’ 

where i' 4 / 4 V = n - 2. There are equations of 

this form which must be satisfied identically by functions 
of the form 

] \dy ] j r * 
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where / + y + A = n. Consequently, of these latter functions 
there are only 2 n + 1 independent. Hence, every complete 
solid harmonic Vn of the degree n is given by the equation 


Vn = 


_ . / d\i f d\j f d\h \ 


where i +j + A = n, and where there are 2 n + 1 independent 
functions, and consequently 2 n + 1 independent arbitrary 
constants. 

Another method of forming complete solid liarmonics 
depends on the consideration that, if cii, j3i, 71 be the 
direction cosines of any line, 






7« 



1 

r 


satisfies Laplace’s equations, and more generally that this 
equation is satisfied by 



It follows from this that the function 



satisfies Laplace’s equation ; and as it is^ a rational, integral, 
homogeneous function of the degree^, containing 2 n + 1 
independent arbitrary constants, every comlplete solid harmonic 
of the degree can be expressed in this! form. 

It is not, however, obvious that a set lof real values of the 
coefficients ai, /3i, 71, &c., corresponding \ to any given com- 
plete solid harmonic always exists, and th^mt in general there 
is only one such set. 
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This proposition, which is necessary to complete Maxwell's 
method of representing solid harmonics, was proved by 
Sylvester {Phil, Mag,, October, 1876), in the following 
manner : — 


It has been shown above, that by tlie solution of linear 
equations for determining the coefficients, we can reduce any 
complete solid harmonic to the form given by (49). 

We have now to show that any rational homogeneous 
function of the degree of the symbols of differentiation 


operating on - can be reduced to the product of n real linear 


factors of the form 


d ^ d d 


dx 


dy ^ dz 


Since the symbols of differentiation obey the same laws as 
quantities, and since 


y/a?* ^ dy'^ dz^j r ' 

the theorem just stated is equivalent to asserting that any 
ternary quantic {x, y, 2 )”, whose variables are subject to the 
condition + y* + 2 ’ = 0, can be reduced to the product of n 
real linear factors, and that this reduction can be effected in 
only one way. 

The equations (a?, y, s)” = 0, ir* + y’* + 2 * = 0 may be 
regarded as representing plane curves having 2n points of 
intersection. If these points be joined in pairs, we obtain 
n straight lines, the coordinates of whose points of determi- 
nation are obtained by solving for x:y:z the simultaneous 
equations (a?, y, 2 )” = 0, and «* + y* + 2 =* = 0. 

X t/ 

If - be real, the corresponding value of ^ given by the 
latter equation must be imaginary, and so also, therefore, 
that of - • Hence of the three ratios, a? : y : 2 , two at least 

y. 

are imaginary. 
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The 

.< //, 2 ', 


equation of the straight line joining the points 
and a;", s" is 


0 ? (/s" - s'/') + // (sV' ^V') + s ' - //V') = 0, 


:o be each imaginary, and select for 
imaginaries, each term in the equa- 
tion of the straight line contains \/-i as a factor, and the 
line is therefore real. y 

If the equation of the degree 2n for determining have 

2m imaginary roots, there are 2m imaginary values of either 

^ X , 

- or - corresponding, and therefore m real straight lines. 

** y y 

Corresponding to the 2 (n - m) real values of there must 

X X ^ 

be 2(n - m) imaginary values of - and and therefore 

{n - m) additional real straight lines. Hence in all there 
are n real straight lines passing through the points of 
intersection of {x, y, z)*^ = 0 and x’^ + y* + = 0. 


y 

If we suppose and -r 
y" x" ^ ^ 

Y? and the conjugate 


y 

There are no more. For if we seek the values of - which 

z 

satisfy the equation of a real straight line, and the equation 
+ y* + s* = 0, these values must be real, or else conjugate 

imaginaries ; and in the former case, the values of - must be 

conjugate imaginaries, and also those of ~ . Hence, to obtain 

a real straight line, each imaginary value of one of the ratios 

V 

- , &o., satisfying (a*, y, zy^ = 0 and a?* + y* + js* = 0, must be 
z 

combined with its conjugate; consequently there are only 
n such lines. 

Let L - 0 denote the equation of n straight lines passing 
through the 2n points of intersection of {x, y, s)** = 0, and 
«’ + y* + 2 * = 0 ; then, whenever («, y, 2 )'* and, a:® + y’ + s* both 
vanish, so must Z, and therefore 

Z = JC (x,y,%Y + Y {a? + y* + a®). 


( 50 ) 
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From the degree of the various functions in this equation 
we see that X is constant, and F of the degree n - 2. * Since, 
in general, a ternary quantio of the degree contains 


[n + 1) (n + 2) 


constants, and the equation of n straight lines 


contains 2n + 1 constants ; and since 


{n + l){n + 2 ) 
2 


2w + 1 + 


{n - l)n 


• in ““ W n 

it ifl plain that the ^ constants of Y can be so deter- 
mined that the right-hand side of (50) shall represent n 
straight lines. It has been proved above that for one of 
these determinations the n straight lines are real. If 
aix + + yiZ = 0, &o., represent these real lines, then 

(aiX + /3iy 4 yis) . . . (a„x + j^„y + y„s) 

= A{x,y,zy' + F((>* + / + 2“). 

Applying the theorem which has been proved for the 
quantities x, y, s to the symbols of differentiation, by (49), 
we get 


r„ = 


d\i f d\J fd\k 


-■'^\dxj \dyj \dzj 

iLYY— VY— YY— — 

ix ) uly ) Ws ) \rfir* di/^ ^ dz^J) r 




(51) 


where V„ denotes any solid harmonic of degree n, A, A^fk, 
and Lify constants, and ai, /3i, yi, &o., the direction cosines 
of straight lines, and where 

i +j + k = n, i' +/ + A' = w - 2. 
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If lines be drawn from the origin, each in one direction, 
having ai, j3i, 71, &c., as their direction cosines, these lines 
meet a sphere, having the origin for centre, in n points 
which are called the poles of the corresponding spherical 
harmonic. 

The mode now described of forming spherical harmonics 
was given by Clerk Maxwell in his treatise on Electricity and 
Magnetism. 

Maxwell’s method of representing spherical and solid 
harmonics admits of an interesting physical interpretation. 

If hi denote a line whose direction cosines are ai, / 3 i, 71 
drawn through the point y, is, and h\ the parallel line 
through the origin, 

( a ^ rf \ 1 d I d 1 
\ ^ dx ^ dy ^ dz ) r dhi r dh\ r ’ 


and 


Vn - 


jL a. 

dh \ din 


d 1 
dhn r 


( 52 ) 


m 


Again, — expresses the potential of a mass m at the 


origin, and 


- mdh\ 


d 1 
dh\ r 


expresses the potential produced by superimposing on this 
mass another negative mass of equal magnitude, situated at 
a point at a distance from tlie origin infinitely small in the 
direction h\. If this system be displaced through the 
distance dh\y reversed, and superimposed on the former, 
the potential becomes 

mdh\dh\ 777- — r and so on. 

dh 1 dh 2 


The repetition of this process n times leads to the 
potential f 7 n, where 

dh\ dh 2 dh n r 

dhi dlh ’ ' ' dhn, r 


7 ^ 

,72 w +1 ’ 


A = mdh\ dh% . . . dhn = M. 


provided 
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If -^1 be a finite constant, m must be an infinitely great 
quantity of the order. 

As an easy example, illustrating the foregoing theory, 
we may consider the question to express in Maxwell's form 
a solid harmonic of the second degree containing only the 
squares of tl\e variables. 

Here, by Art. 161, the solid harmonic 

Fa = ax^ + bif - {a -vb) 
if (P 


( 



iP (P 

ay dz 





and therefore 




^ (P\l 

dx^ ^ dif j r 


= 3\x^ + Spjf - (X + ju) = (*2X - /i) .r* + {2fx - A) y* - (X + ju) 


Hence, comparing with the former expression for F, we 
have 2X - /x = «, b; whence 

X = la -h III ^ la lb. 


To reduce ^X -- + p ^ j - to Maxwell's form, we must 

consider the relative values of X and p. 

If X and p have different algebraic signs, and p^- p\ then 


Fa = 


d /_ d 


- d 
'^^dx 



If X and n have the same sign and X he the greater, 

r 




dt/^ ' dzy) 

- p ^ - y p^)-- 
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The most general solid harmonic of the second degree is 
reducible to the form considered above by a transformation 
of the axes of coordinates. 

In general, for the second degree, the reduction of a solid 
harmonic to Maxwell’s form is mathematically the same 
problem as the determination of the planes of circular 
section of a quadric surface. 

163. Incomplete Harmonics. — We have seen, Art. 
150, that if Yi be a complete spherical harmonic whose degree 
is real, i must be a positive integer, and Yi a rational and 
integral function of /i, -^/(l - /i*) cos 0, and v^(l - /i*) sin 0. 

If i be a negative integer, -(/+ 1) is zero or a positive 
integer. 

If i be real but not an integer, it is easy to obtain expres- 
sions for Yi whicli satisfy the differential equation (45) ; but 
these expressions become infinite at certain points on a sphere 
surrounding the origin, or alter in value after having passed 
continuously through a complete circuit surrounding the axis 
of z. In the latter case, accordingly, they are not single- 
valued. 

If we assume pi = + ^ 2 ju* + + &c., 

qi = ai/m + a^jY + a^p^ + &c., 
and substitute pi in tlie equation 

in order that each power of p should vanisii, we find that 


<"^n+2 — 


{i - n)(i + ?< + !) 


(54) 


Equation (64) is fulfilled also by two successive coefficients 
in the series denoted by qi, provided qi satisfies (53). Hence, 
if we assume 


Pi 


Qi^ (lx 


^•(^ + l) ^ /(e-2)(e + !)(/+ 3) 




T.2 ^ 

(iziX'll) 

2.3 


tx^ + 


1 . 2 . 3. 4 
(»-l)(/-3)(j+2)(/+4) 
2.3.4. 5 


- &c. 




(55) 
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each of the series and qi satisfies the differential equation for 
a spherical harmonic of the degree whatever be the value of i. 

If i be an integer, one of these series terminates : the other 
contains an infinite number of terms. 

If / be not an integer, both series contain an infinite 
number of terms. 

The sum of each series is finite so long as ju < 1 ; but if 
^ = 1, either series, if it contains an infinite number of terms, 
becomes infinite. In fact, (54) may be written 

(n - %) (w -f i + 1) 

(« + !)(« + 2) 


and, accordingly, as n increases without limit, all the terms 
become of the same algebraical sign, and the value of 

dn 

tends to become unity. Hence (Williamson, Differential 
Calculus^ Art. 73), if ju < 1, the series is convergent. 

In the case of the more general spherical harmonic F,, 
whatever be the value of we may assume 

» 

Yi = S {As cos 50 + Bs sin 50) (1 - [Jb^Ypu ; 


then, as in Art. 148, equation (30), we find that must 
satisfy the equation 

(ff - 1) D^p + 2/x (5 + 1 ) 2>i> - (^ - 5)(^ + 5 + 1)^ .= 0. (56) 

This equation is satisfied by the series 

and by the series f/iju + + &o., 

provided that in each series 


{i - 5 - n){i + 5 + w + 1) 

(« + 1 )(« + 2 ) 


Henoe, if we assume 

li-s){i + a + l) 
.astapr — 

teristi^'* *■ * ^ 


+ &o. j, I 


f -8-l)(t+8 + 2) , „ \ i 


(57) 


(58) 
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Yi = S {As cos 6*^ + Bt sill (1 - p^YiPh + ?«), ( 69 ) 

where i and s have any values whatever, Yi will be a spherical 
harmonic of the degree i. 

If i - 6* be a positive integer or ^ + s a negative integer, 
one of the series p* and j,-, terminates, and the other contains 
an infinite number of terms. 

In any other case, both series contain an infinite number 
of terms. 

When the number of terms is infinite, 

(1 - p^YPi* (1 - 

are each finite if /i < 1 ; but if ju = 1, each of these expressions 
becomes infinite. 

In order to prove this, we observe that ( 57 ) may be 
written 

+ 2(.s- l)n + (s + « + 1) 
an+2 = (h9) 


When n becomes very great, the ratio of to an tends 

towards 1 + ^ • 

n 

Again, if we put (1 = 1 + + 64^* + &o., we find 

that 


+ 2 4 - 2 (A - 1) 

;rr2 — 


As n becomes very great, the ratio of 6n+2 to hn tends towards 

n 

Hence, as p approaches 1, the functions pu and qu tend 
to become quantities of the same order as (l-/ii*)”*; and 
therefore, if 8 be positive, 4 _ , 

1 ^ 

(1 - p^Ypis and (1 - p^gis 

are finite so long as p<l; but it p^l, they become infinite. 
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The same thing is true if s be negative. In this ease, 
$ 

becomes infinite when ju = 1, and the products 

s f 

(1 ~ /u^Pis and (1 - 
become infinite as before. 

When i = s, we have pis = 5 when i = s + 1, the 

series gi, = 

It appears from what has been said that if we assume 

s* 

Yt = {A cos icji + B sin {(p) (1 - ju*)*, 

where i is positive, Yi is always finite; but if i be not an 
integer, Yi is not single- valued, for when 0 increases by 
27 r the functions cos /0 and sin i<p do not return to their 
original values. 

If ^=3^, Ti = (1 - ju»)* sin 


and e = ns, where n is any integer, and 5 = - , the function V 

satisfies Laplace’s equation, vanishes at infinity, and is zero 
at the planes for which 0 = 0 and 0 = a. At the surface of 
a sphere of radius a we have 

a'*' 

Thus on this sphere F is a function of fx and 0, which vanishes 
at each of two great circles, and is finite and single-valued 
for the intercepted portion of the surface. By bringing in a 
sufficient number of terms, and properly determining the 
arbitrary constants, it may be possible to make this function 
equal, at least approximately, to an assigned function having 
the characteristics above. 

In some oases, a function satisfying Laplace’s equation, 

find fiilfillinoi' nArf.fi.in hnnndarv nnndifinns. nan Via found bv 
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means of spherical harmonics of imaginary degrees. We 
have seen already that whatever be i and s, 

X 

S {As cos 5^ + Bs sin S(p) (1 - jx^) ^(pu + qis) 

is a possible form of a spherical harmonic of the degree t. 
If « be imaginary in order that pia and qia should be real, it 
is necessary only that 1) should be real. If we put 

i (i + 1) «/, we get / = - J ± + i) ; wlience, 

if i be imaginary, / must be negative and greater than 
J in absolute magnitude. Accordingly, putting /=-A;, 
we obtain 

i ^ + y ri y (A - ^), /' = - j - y ri y (/, . j) ; 

then, since Yi depends only on the value of + 1), wo have 
Yi = Fi', and both these functions are real. 

If we now assume 


r=r^F, + ri'Fi', 

the function V satisfies Laplace’s equation, and we have 

F- 2F 

F = {r» + r»') Yi = — (^x + ^xVT) =. —I oos x- 

Vr y/r 

If we assume \/-l F'= r*Yi-r*' Yvj we get in like manner 
2Yi . 

F' » sin x> where x “ i) order that 

y/r 

F or F' should vanish at a sphere of radius we have only 
to assume 

2n + 1 


2 


TT, or Xa ** 
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In the first case, we have 
h 

and, in the last, 


Incomplete spherical harmonics are liere briefly described 
in order to give the student an idea of their nature and of 
the kind of conditions which they can be made to satisfy. 
Tliey are useful in some departments of mathematical physics. 



Section II . — Ellipsoids of Revolution, 

154. Solutions of Dllferentlal Equation. — When 
the surfaces with which we have to do are not approximately 
splierical, the expansions for the potential which liave been 
investigated are of little use. In the case of ellipsoids of 
revolution, equations (35) and (41), Art. 98, enable us, by 
an extension of the theory of spherical harmonics, to arrive 
at suitable forms for tlie potential. 

Equation (41), Art. 98, if we write ^ instead of x> 
becomes by transposition 

d dv 1 d^r d dv 1 cpv 
dV^ dK ^ 1-5* rff “ dV dZ df^ 

( 1 ) 

If such a form be assigned to V as to make each member 
of this equation equal to the same quantity, the equation is 
satisfied ; but, by Art. 148, if P,, satisfy equation (22), Art. 
147, then 

s 

(1 - ooastp + B sin s^) 

satisfies equation (27), Art. 148, if substituted for Y„. 
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Hence 

(1 - (I)' Pn{%)Pn{t.){A 008 8^ + 5 sin .90) 

must satisfy (1), and if V satisfy Laplace’s equation through- 
out the region inside a prolate ellipsoid of revolution, we may 
put 

V = ^TnB{^)Tna{K) {^s OOSS^ + B, sill S(^). (2) 

The value of V given by (2) becomes infinite along with Z 
at points at an infinite distance from tlie centre of the 
ellipsoid. Accordingly, (2) does not give a suitable form 
for V in the space outside the ellipsoid. 

It appears, however, from Art. 147 that there are two 
solutions of equation (22) of that Article. One of tliese 
is Pn] til® other, which may be denoted by Qn, contains 
only negative powers of p* Accordingly, when t becomes 
infinite, Qw(?) becomes zero. 

Hence, if we put 

3„= (^^y(^~Jp„(|)(Q)„(Oi>4,coss0+£.8i.> 80), 

( 3 ) 

and V - S2)nj (4), we see that this form of F satisfies 
Laj)laoe’s equation tliroughout the space outside the ellipsoid 
and is zero at infinity. 

If we denote by Um the function corresponding to Tna 
in equation (33), Art. 148, we have 

(5) 

and 

Un = STns{K)J7nt{Z) (^s COS + Bg 8111 S(p). (6) 

At the surface of the ellipsoid, wliere 2 is constant, 
§)« becomes a spherical harmonic F„. 
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155. Determination of tl^e Function Q^. — The 

differential equation (22), Art. 147, being of the second 
order, has two partioular integrals ; one of these is Pn, the 
other Qn- Putting Z instead of fn, by Art. 147, we have 


^ 2(2n+3) ^ f2.4(2//+3)(2;i+5) ^ 

f 

By Art. 147, if as be tlie oHefficient of in this series, 

“ la oT / . . o\ /" T\ ^ 


■3s-}n + ^ 

Hence, as - s increases, t^e ratio a ,.2 : Og tends to become 
unity, and if ^ > 1, the seric;|j is convergent; but if ^ = J, it 
is divergent. Hence, in tip space inside the ellipsoid of 
revolution, 2 )m is not a suitiible form for the potential. 

156 . Dc^Aerinlnatlon Q,» as an Integral. — If we 

write y for 6n and ,v for ju, e4nation (22), Art. 147, becomes 


2^7 (///I u(n + l) 


da} 1 di 

which is of the form 




07^-1 


^y = 0, 


d’^H aV V c\ 




( 9 ) 


(10) 


If we put y = v>/t, equ^on (10) becomes 

\ dx 


\ dr 


'(S + -f- s * I * i ■ »■ 


If yi be a solution of (10), we get 


( 11 ) 



dr 

dx 



Xidx\ 


( 12 ) 
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whence, by integratiot, 


v^€,+ C, 


^ — J A 1 dx 


. //«' 


dx. 


and if t? = 0 when x = a?,,, have 

% 

rx dx 

V =. CA — — dx. 


I aro 


In tlie present case, 
2a; 


xr 1 X I -Xl dx 2 1 

A , = > and ^ = I??' - 1 ; 


also, Pi » Pn- Accordingly, ^ 


t) = O', 


f* eo^-r 


1) P\ 


(13) 


(14) 


By putting a? = expand iig the exiiression under the 

integral sign in ascending poweVs of and integrating, it is 
easy to see that when x= oo, or s = 0, we liave v = 0. Hence 


Qn = 


O^P„ 


r sujt 


If we choose - 1 for tlie vajnfi./(tf C\, we make Qn perfectly 
definite, and we obtain 


Id 

Q«= Pn 

J * 


dx 




(15) 


y 157. Kxpreseion of Q„ by means of a Finite 
llerles. — In order to express Q„ as a finite series, it is 
necessary first to prove some relations which exist between 
successive coefficients of Legendre and the functions obtained 
from thm by differentiation. 
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If we put - 1 = 


= D, we have 
dx 


Du = 2x, D^u = 2, D\i = 0 ; 

and we get 

2 ^n+ 2 ^«+l ^ 2 )n+l 

= (w 4- l)(I>c/0”+i + (n + l)D'uD^Hi*^) 

= 2(w + W” 2(« + 

Substituting for by a formula similar to that just 

obtained, we get ^ 

2)n+2^n+i « 2(fi + l)ir { 24 k«? D" + 2m* ) 

i + 2 (m + 1)*D" ; 

but by (21), Art, 147, we haj/e 

1 

and therefore, dividing by j^(M + 1), we get 
DFn+i = X^DPn^i nxPn-i + (m 4 1)P„, 


and subtracting i?Pn-i> we have 

PPn+l - PP «-1 = uDP,,.^ 4- HxPn-i + (w + l)Pn. (16) 

The right-hand side of (16) can be expressed in terms of 
P„, for we have 

!)»«» = Z>”M M”-! = + 2«i > • ««-» + ~ ; 

. 1 . 0 . T 

=» D^’^Du^ = P”’^M u^'^Du - 2 Ma?P”“^M””^ 4- 2 m(m 


Comparing the two expressions for we find 

= m(m - 1)^”“*^”*'. (17) 

Equation (17) shows that iPu^ satisfies (22), Art. 147, 
a result which has been already proved in Art. 147. If we 
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now substitute for n{n - 1)2)”“*?/**"^ in the first of 

the expressions for given above, we get 

D^u'^ = (18) 

and dividing by 2’* I w , we ^ave 

1 

Pn = - uDEr.i + ooPn^i, (19) 

whence w2)P„-, + iixP^^i = «Pn- 

Substituting in (16), we obtain* 

DPn,i - BPn-i r {2n + 1) P„. (20) 

From (20), we get immediately 

BPn * (2w “■ IjPfi.i + {2u — 5)Pn-3r • » • + {2n — 4»s + 3) P»-*j+i + &o# 

• ( 21 ) 

We can now express Qn as ti finite series by treating the 
equation 

D(uDij)-n{r.^l)y^Q (22) 

in a manner somewhat different from that previously 
employed. If we put y^vijx-w^ and substitute in (22), 
we get 

0 -= ^(ir*-l)^^- ;a(w + l)?7x- (P(wPw?)-w(« + l)e^} 

- ' Is j * S’ s 
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If we next suppose y\ to be a solution of (22), and 
determine v in such a manner as to satisfy 


la-rrn- = o 

dx ^ ' dx ’ 

we get , T. r. / 1 M o / , 1 \ 

{JDul)-n(n+l)Uvj2{x^-l)——< 

From (23), by iiitegratiou, we obtain 

(l-x'^)-!^ = constant. 
dx 


(23) 

(24) 


If we choose 1 as the valne of this constant, wo get, by 
integration, 



II 

1 + 

»— * h-a 

(26) 

and (24) becomes 




/ -.V r. 

(26) 

If we assume 



W = 

AzPn~i ... 1- A'tg-iPn-zs^i + &0., 



and make t/i = Pn, by (21) of the present Article, and (2), 
Art. 138, we get 

2 (2n - 4s + S) 

= At,-i(n(n + 1) - (« - 2s + !)(;< - 2,i + 2)} 

- {n* + »/ - [«“ - (4s - 3) » + 2 (s - 1) (2s - 1)] ) 

- Au.i { (4s - 2) « - 2(2s - 1) fs - 1) 1 = 2 (2s - 1) {« - s + 1) 


whence 


2n - 4 s + 3 
" (2s-l)(«-s + l)’ 


(27) 


2m - 1 jj 2 m - 5 

~T7n 3(m-1) 


2w - 9 
5 (n - 2) 


•Pn-5 


+ & 0 . 


( 28 ) 


to = 
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Thus we obtain 

(?„ = r7jp,41ogl^- «’j, (29) 

Where w is given by (^8). 

We have now arrived at ^iree expressions for Qn of which 
one is perfectly definite, and^he otlier two contain constants 
which can be determined by comparison with (15). 

From (29) and (16) we get 


C-I^Pnlog 


X - I 



dx 


whence, by division and differentiation, we get 


that is. 


i __ 

1 1 - 


dw dPn \ 
dx ^ dx ) 1 

" (1 - 


0 1- 


1 - 

~p 2 
n 


/„ .IP. 

V" T. - " rfi 


1_ 

P 2 > 
^ n 


but when X = 1 we have P« == 1, and therefore (7=1. 

In order to find K in (7) we put « = - , and identify the 

X 

coefficients of the lowest powers of z in the expressions for 
Qn(i^) given by (7) and (15). 

By Art. 147, we have 


Pn = A Ix^- 


n{n - 1) 
2(2/^-!)'^ 


+ &o. 


* Axe'll + ascending series in »}, where A = 
therefore, 


L« 



dx 


{X? - 1)P„’ 


= Ax*‘{l + &o, 


I 

•) 

J { 


dz 


(I + &0.) 


s’"( 1 - s*)‘‘ (I + &o.y*dz. 

0 
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Here the lowest power of % after integration is plainly 

and its coefiGicient is 1) ' 2 ’ " (2w + 1) Jl * 

that is, . ' . 

/f = • ■ (:i0) 

For Q{?) we liave, then, three expressions given by the 
equations 

' 1 .3.5 . ! r(2«TI) 

r 2(2»+3) ^ 2.4.(2»+3)(2«+6) ^ )’ 

(31) 

= (32) 

Q„(^;=^P„(01og|^ 

j2«-l 2»-6 2?« - 9 . ) 

» 1 . // 3 . (« - 1) 5(u - 2) ) 


158. — Analogues of Tesseral llaruionics. — We saw 

s 

in Art. 148 that the multiplier of cos 5 ^ in the spherical 
harmonic Yn must satisfy equation (30) of that Article. This 
equation has two particular integrals, D^Pn and IPQn] but, 
by means of (30), Art. 148, the latter can be expressed in 
terms of the former. In fact, (30), Art. 148, may be written 

(*» - 1) g + 2(s + l)x I + {«(« + 1) - n (« + 1)) y = 0. (34) 

If we compare this with (10), Art. 156, we see that 
^ 2 (5 + 1 ) IT 

and consequently 

^-SXidX ^ _ 1 )'(»+>). 
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Accordingly, if yi be a particular integral of (34), the other 
particular integral is given by the equation 


Vi = Cijx 


dx 


<1 


and therefore, adopting the notation of Art. 154, if we put 

s 

T„, » we have 





dx 

y,\x^ - ir' 



dx 


(35) 


We may regard TJ^t as defined by (5), Art. 154, and 
consider as given by the equation 

It is now easy to show that in (36) the value of C is 


(- ly 


n s 


For 

7;. = = 


~ .V 

1.3.5... (2n-l) - 


L!! 


(a;” + &o.) 


u^(x^ * + & 0 .) 


1.3.5. . .(2;^- 1) 


]^l - 8 

and if we put 

1.3.5.. .(2w-l) 
[^w - s 


+ &c.); 


« and “ = 2 , 

X 


we get Tn, « Nx^{l + Z\ where Z denotes a series in 
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(*“0 ({% 

Nx (1 + Z)^ ^ z^x^{l-z‘‘')N^oi?\\.+ 


Zf 


C«o 


but 

and 


^(1 + Z) I _ ,»(1 - ^)-(l H. So., 

("36) 

L!i 


Qn = 


+ &C.), 


1 . a . 5 . . . (2n + 1) 
i - \n (n-\-l)(n+2 ) . . . (n+s) 

Uns = U^D^Qn = - 2')' 

l.d.o . . . (2/^^-l) 

I ^ + .9 

= (- ly L== 

^ ^ 1.3.6 . . . (2n + l) 

Hence, equating the coefficients of the lowest powers of z in 
the two expressions for Uns, we get 


■=(- ly 


L 


n s 


(2n+l)]Sr ' '' 1.3.5 . . . {2n +l)^ 


whence 


\ 71 + S 

c = -(-iyh=; 


and as «o = therefore Xq - oo ^ we have 

f n + S r" ,7^ 




(37) 


It is to be observed that, in order to avoid the introduc- 
tion of imaginary quantities, Tns (i*?) is regarded as having a 
somewhat different signification according as x <1, or x>l. 
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In fact, in Art. 148, 
Tnsip) - 


and, in the present Article, 




(38) 

but 


(39) 


It is obvious that, whichever signification be attributed 
to Tns{x\ it satisfies the same linear differential equation. 

159. Expansions for External and for Internal 
Potential. — We can now write down the series express- 
ing the potential, inside and outside a prolate ellipsoid of 
revolution, due to a distribution of mass on its surface. 

Let V denote the potential inside, and F' that outside, 
the ellipsoid whose semi-axis major is ; then we may put 

+ s(^„, cos + Bn. sill ^ ’ 

and (40) 

2 cos S<p + Bn, sin S<j>) Tn,(X) ^n* (?) (C)l*| * 

(41) 

At the surface, where 2 - ?o» these two expressions become 
the same series of spherical harmonics which can be made 
equal to any assigned function of ^ and (p which is finite, 
and single-valued. 
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160. Surfoce Distribution corresponding to 
Potential. — If the internal and external potentials, V 
and V', due to a surface distribution whose density is a, 
be given by the equations 

r= ^ooss^rc?) r(C)j* 


V' - A cos ^ 




dt. 


<T can be found from the equation at the surface of the 
ellipsoid So* 

If dsi be an element, drawn outwards, of the normal to 
the surface, equation (12), Art. 46, becomes 

dV dV' . 

dsi dSi 


but by Ex. 2, Art. 75, dsi = 


ldZ,dl 

P ’ 


and therefore we have 


dV dV' . 


P 

dTC < dT[ 

A 


, fr(y\ 1 

-+ i ((,(,> ,y-j- 


{W-i)\T{WS 


T{^) cos ; whence 


^ Ap ^ _ AT(^) 00 Bs<t> 

(42) 

When the density of the surface distribution is given, 
the expressions for the potential inside and outside the 

ellipsoid may be determined by expanding expressed as 

Jr , 

a function of ^ and 0 in a scries of spherical harmonics, 
and determining each of the functions 2 (£) and the corre- 
sponding constants by means of (42). The potentials 
and F' are then given by (40) and (41). 
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161. Potentials of Homoeoid and Fo4salold. — If 

the surface distribution be homoeoidal, the density varies 
as p, and the multiplier of p on the right-hand side of (42) 
must be constant. Hence, V = constant, 




1 1 i 

2 ^3 S’ 


When Z becomes very great, tends to become — , 
M 

that is j-yi where M denotes the total mass of the homoeoid. 


Hence C? = -y- , and 
k 


„ *■, «.+i 


For a fooaloidal distribution the density varies inversely 

as jt), and - varies as -4, that is, as ~ Accordingly, - 
p p p 

is of the form AP^ + BP^ (^), and 




” dZ 

i (S’- W-iy 


It is easily seen that 

1 V 1 1 ^ 

S*-i 2VS-1 S + V’ 

and that 

1 9/1 1 1 




Hence, by integration, we have 


r = ^log^ + j(r-i)iog|^ - 2S 
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By a method similar to that employed in the case of the 

M 

homceoid, we find that where M is the mass of 

the fooaloid. 

It follows from (42) that the coefficients of Po and (S* - J) 

in 47rA*So(^o* “ 1 ) ~ are rTI T i “ varies as fo*- 

that is, as i)> therefore we must have 

A 2 = - Aq, Hence 



(46) 


By Art. 83, the potential of a fooaloid in external space 
is the same as that of the solid ellipsoid of equal mass of 
which it is the boundary. Hence (46) may be verified by 
comparing it with (2), Art. 78. This verification is readily 
effected by taking a point on t^e axis of revolution. Here 

5=1, and r = kti\ then putting ^ = z, from (46) we get 



(47) 

Again in (2), Art. 78, for points on the axis of revolution, 
J= A ; and, since C = B, we have 


A + B + C-U 
2 


= 5 - 


A = 


6 


Hence, putting A? for r, we get 


r= 


M M 


which agrees with (47). 
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162. OMate SSIlipsotd of Revolution. — When we 
have to do with oblate ellipsoids of revolution, Laplace’s 
equation takes the form given by (35), Art. 98. If we put 
: = -rv/^ in this equation, and write 0 instead of 
we get an equation in C, and which is the same as (1). 
Hence, in the case of the oblate ellipsoid, we may put for V 
and V' the expressions given by (2) and (4), or by (40) and 
(41), provided we put instead of t* In order to determine 
completely these expressions for V and V\ we have then to 
put 1 for and accordingly we have to find what P{K'\ 

T(0, become when ^^-1 is substituted 

for r. 

163 . Retermination of p (Z) and q{Z)> If we put 
instead of Z in the expression for P« (2^), given 

by (26), Art. 147, we get 


(- If A 


n {n - 1) 

^ ^2(2w-l)^ 




and we may put 




Also, putting \/-l = i, we have 

PM = 

In like manner, from (31) and (32), we have 

ll 


(48) 

(49) 


(50) 
dK 


1.3.6 -(2« 


1 ) {PM)^ 


( 61 ) 
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It is plain that the right-hand members of (50) and (51) 
can differ only by a factor which is some power of t from the 
expressions for Q(t?I), given by (31) and (32); and, as the 

term involving the lowest power of ^ is the same in the two 

forms of qn{Z)i they are consistent. 

From (31), we see that 


Qni^r,) = 


(52) 


In order to find a third expression for 5^n(?), we must 
consider what log becomes when $ is changed into 


If we put 3 = 2 = tan 6, we have 
4 


, ■ 1 , 1 - tz , cos 0 - I sin 0 _ . 

log ITTi = log = log t-t; = - 2i0. 




cos 0 + ( sin 0 


Hence, 

and 

also, 


i log tan-‘z, 

log = c”^‘p„(Q tan-> ^ ; 


P„..(«?) = 
Accordingly, by (33), 

1 /2»-l 


Qn(i^) = jij»(?) tan-‘ I - + & 0 . jj ; 

but qn(0 = i*"** = (- 1)“ ; and therefore 

?«(?) = (- 1)" jp«(2)tan-»i-^^p„.,(?) + ^^^p„j(^)-‘«-«l 
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164. Analogues of Tesseral Harmonics. — When 
we put c? for 2 hi T„,(Z), we get 

S 

+ lyiypniZ), 

and we naay write 

= (V + ifii’Pn (?), TM = ( 54 ) 

Also, we may write 

s 

+ (55) 

whence Z7«s(<2) = «"(“*’>«(„»(?). (66) 

Another form of «„,(0 is obtained from (37) from which, 
by means of (56), we have 

Mi»(0 = 1)* 


Ll 


whence 


Liz. 




-idZ 


Uns(Z) = (- 1 ) 


+ S 
f H — S 


fn,a) 


f + 1) MZyr 

dl 






(67) 


165, Expression tor Potentials. — If V and V' denote 
the potentials inside and outside an oblate ellipsoid of revolu- 
tion, due to a distribution of mass on its surface, we may 
write 


+ S (A„, cos + B„, sin s<p) J 

F' = ^A„p„{K)pn(Z) 

+ S(-4„, oos s0 + B„, sin s^) 7’„,(0 


dK 


( 68 ) 


dZ 

(^ + 1){M0}*' 

(69) 
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It is plain from what precedes that V and V' satisfy each 
Laplace's equation, that Fis finite when ^ = 0, and V' is zero 
when 2 = 00 , and that V and V' are identical at the surface 
of the ellipsoid. Hence they satisfy all the conditions 
required. 


166. Surface illtiitributlon corresponding to 
Potential. — Here we may proceed as in Art. 160. If the 
internal and external potentials V and V' be given by the 
equations 


F'-^™,0no<(o|"(jrnf7(oTv 


(60) 


and dsi be an element of the normal to the ellipsoid, 

dsi — j 

p 


but in this case, by Art. 98, we have X* = A:® (5* + 1), and 
therefore 

dS\ — • 

p 

Accordingly, as in Art. 160, we get 




Ap AT{^)ooas^ 

+ l)t{Q “ ktlQ^ ((^o^+lXC+rM’ 


When the density of the surface distribution is given, 
the potentials may be determined in a manner similar to 
that described in Art. 160. 
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167. Potentials of Homceold and Focalold. — 

These may be obtained from (60) and (61) in a manner 
similar to that employed in Art. 161 ; but if the expressions 
given by (44) and (46) are already known, we oan got from 
them the corresponding expressions for an oblate ellipsoid. 

Putting for 2 in (44), we get 


M 

k 



Hence we may conclude that V\ the potential of an oblate 
homoeoid of revolution, is given by the equation 

(62) 

k Z 

It is easy to see that this expression for F' must be 
correct, since it satisfies Laplaoe^s equation, vanishes at oo , 

is constant at the surface, and tends towards the value ^ 

at points very distant from the centre when r tends to 
become equal to kZ* 

To get the potential of a focaloid in external space, we 
may put iZ for Z in (46), and we get 





since < = 4 ®, and £* = - 1, this may be written 



Hence we may put 





( 63 ) 


V' 
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This expression for V' satisfies Laplace’s equation, vanishes 
at infinity, and at points very distant from the centre tends 
M 

towards the value ~ ; but to prove that it satisfies all the con- 
ditions of the question, we must show tliat the corresponding 

distribution of mass varies as - • 

P 

If the external potential be given by (63), the internal 
potential F, due to the surface-distribution producing F', 
is given by the equation 


F = 





> 


(64) 


since this expression for V satisfies Laplace’s equation, 
remains finite inside the ellipsoid, and is equal to V' at the 
surface. We have, then, 


P 

M ( 



dV 


1 





(65) 


Differentiating and reducing, it is easy to see that 

U ^ L /^tan- 1 - V ^ li 

Hence, 

dV dV' MS \ \ M 

d\;~'d% - k li + ^‘'(i + r)(r + J))"* (i + f^Xr + i)’ 
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and at the surface we have 

dV dV M + C 

p dti - *(l4CoO(?o* + i) f' 

where Q denotes a constant. 

C 1 

Accordingly a = . and therefore the distribution 

Lq p 

of mass producing the potential is fooaloidal. 


Section III* — Ellipsoids in QeneraL 


168. Ellipsoidal Harmonics. — When the surface, at 
which the potential or mass-distribution is given, is an 
ellipsoid not of revolution, the preceding methods are inap- 
plicable. The most general method of determining solutions 
of Laplace's equation which can be made use of in questions 
of this kind depends on the employment of functions called 
ellipsoidal harmonics. 

We have seen, Art. 92, that if X, ju, v be the primary 
semi-axes of the three confocal quadrics passing through a 
point, Laplace's equation may be written in the form 


V V d^V 


where a, j3, y are given by (17), Art. 92. If «, 6, e denote 
the semi-axes of an ellipsoid of the confocal system, we may 
change the variables by assuming 

X* = a* + = a* 4 - ij, i;* = rt* 4 - 2^. 

If we put 3£ = (a* 4* 5)^ (6* + (c* + by (17), Art. 92, 

we have 

d dK d% d 21ji d 
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If we assume 

2J = (a» + ri)^ {V‘ + (c> + r,)K 3 = («’ + (c’ + 

in like manner we get 

d 2g)v/^ d /ON 

d/3 k dy\ k dX} 

and Laplace’s equation becomes 

{(, - ?)(*!)■+ <? - 5)(30 * (I - .)(3|)‘l r - 0. 

(4) 

Following the analogy suggested by the methods em- 
ployed in the case of ellipsoids of revolution, we may 
suppose V to be the product of j&, a function of and of 
H a function of ij and H' a function of If these functions 
be such that 

= {ml ^j)E, (2)0H = (»^„ +yjH, 

and ^3^j H' = (m? j)H' , (5) 

where m and j are disposable constants; we* may put 
V = CEHH\ where G is an arbitrary constant, and we 
have 

K+y)} F = 0, 

where A denotes the other factor of V* F, since the expression 
inside the bracket vanishes identically. 

We have now to find forms of the functions H, and 
H which will satisfy (5). 

If we suppose £ to be a rational and integral function 
of we may put 

E--^^ + + + 
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Operating on we get 

M((«» + £)(6^ + g) + (6»+g)(c» + g) + (o» + g)(«» + g)} 
2 


= n {n 4- J) + n^{a^ + 6* + c*) + n {n - + (?V) 

+ n (n - 1) 

Hence, in order that 

(x^jE = {mE+J)E, 

we must have the series of n + 2 equations 
n (n + J) = m, (n - l)(n- ^)pi + n^ (a® + + c®) =y 

w (n ~ I) (a®&® 4- 6 V + cV) + (n - 1 )® (a® + 6® + c®)^i 

+ (n - 2) (w - 1) jh - &o., 

2a}b^c^Pn-2 + J + 6®c® 4* c^a^)pn-x = 

The first of these determines m, the second as a linear 
function of/ By substituting this value of px in the third, 
Pi is determined as a quadratic function of / and so on. 
Thus, finally, an equation of the degree ?/ + 1 is obtained 
to determine/ Each root of this equation corresponds to a 
set of values of joi, j» 2 , &c., which furnishes a function of the 
required form for E, There are thus n 4- 1 functions of tlie 
degree n in \ which are of the required form. It is plain 
that I is of the second degree in the coordinates a?, y, z. 
Hence the forms of V corresponding to those found for E 
must be of even degree in each of the coordinates a?, y, z. 
To determine forms of odd degree in these coordinates, we 
are guided by the formulae for expressing the Cartesian in 
terms of the elliptic coordinates of a point. 
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In fact (Salmon’s Geometry of Three Dimensions), 

> = (<»'’+ + K + ?) , ^ (6^+g)(y+ii) + 

(6* - a*) (6" - c») ’ 

(c* - a*)(c* - 6") ^ ' 

If we consider only the factor containing we see that x 
corresponds to 

+ $ , y to v^6“ + and z to -v/c“ + S ; 
and we are thus led to consider whether 

and y(a* + S)(i* + |)(c’ + f)», 

where iS is a rational and integral function of are possible 
forms of E. 

Operating on 5 we get 

(x^)V(«* ■* ^ ^ 

+ [«V + (n + + [n {n - J) + c*) ^ * (7) 

+ w(w + ^) 

Hence, we see that if i? = ^ iSn» where ffln is a 

rational and integral function of S of the degree w, we have 

and therefore, that by properly determining m,y, and the 
coefficients in fin? we can make 
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It is easy to see that the final equation for determining j 
is of the degree w + 1. 

If we operate on the functions 

+ + and 

we find that 

+ [(w 4- {a^ + 6*) + {n + lyc^"] + [n(« - I) j ’ 

+ n {u + J) (a* + ¥)c^'\ + n (n - 1) ) J 

( 8 ) 

and that 

(3e|)Vl («’+?)(*’+ 5)(c‘+S)K'‘ 

= X { (n + !•) (n + 2) + (n + 1)* {«* + 6* + c^) + n ' 

X (n + (a*6* + iV -f c^a^) + ;/ (w - 1) ^ 

Hence we conclude that 

(*J)V(«’ + 5)(i’+4)s« = v^(fl* + ?)(i'>TI) s„«, 

and that ^3£ ^ J 3E S„ = 3eS„« , 

and, accordingly, that ^/l(a^+ 5))Sn and 3£Sn are 

possible forms of E. 

If E be of the degree v in £, and v be an integer, the 
forms we have found for E are 

{(^* + 5)) iSv-ly {(i *+ + f)} Si/.i, 

v/{(c®+ S)(a* + 5))S,.-i. 

We have found also that there are v + 1 different functions 
of the first type, and v of each of the others ; so that there 
are 4y + 1 in all. 
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If V = n + where n is an integer, the forms found 
for E are 

Also, there are w + 1 functions of each of the first three 
types, and n of the fourth, so that there are 4n 4 3 ; that is, 
4v + 1, in all. 

Hence, in any case, there are 4v + 1 determinable func- 
tions of 2 of the degree v, any one of which may be taken 
for E in order to satisfy (5). 

It is plain that if H be the same function of and H' 
of as E is of the product CEHH\ where C is an 
arbitrary constant, will then satisfy Laplace^s equation, 

169. JEllfpsoidal Harmonics which vanish at 
Infinity. — The functions considered in the preceding 
Article do not vanish at infinity, and are therefore unfit to 
represent the potential of a finite mass in space outside itself. 
The form of the differential equation for E enables us, 
however, to obtain another function which will fulfil this 
condition. 

In fact, if 

y ^ {ml+j)y, (10) 

we have 

(^y dp we+y 


but, as was shown in Art. 156, by assuming y 
be a solution of the equation 


then 


dx^ * dx 


+ X^y = 0 , 



0-^X.idx 

yi 


dx 


yiw, ilyi 


is also a solution. 
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In the present case, J Xidx = log H ; and, therefore, if E 
be a solution of (10), so also is 


If now we take for E one of the forms found in the last 
Article, by writing E as the product of 5*' and a series of 
descending powers of it is plain that 


vanishes when | is infinite. 

Hence we see that, if V denote the potential inside an 
ellipsoid, whose semi-axes are a, b, <?, of a distrib of mass 
on its surface, and V' the potential of the same c 'tion in 
external space, and if F = CEHH\ then 

r * 

.J 

where 


^ Jo 


170. Kllllpsoidal Harmitaics expressed as fiine- 
tions of Cartesian Coordinates. — If Sin be a rational 
and integral function of S whose factors are S - S ~ ag, 
5 ~ a«, and E = Sny then 


CEHH- («*+ ai)(J*+ ai) (c*+ ai) 
(a*+ag)(i*+ag)(c*+a2)| 


' a?* 


s* 

L«* + ai 

+ + 
6* + oi 

C*+Oi 



2 * 

a® + ag 

6® + oj 

c® + ag 



JS* 

a*+ ttn 

+ — + 

^5*+an 
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For the expression 

(a* + m) (J* + u) (c® + u)( — + r — — + 

' ' ^ ^ ^\a*+ u h*+ u 


c* + M 




vanishes when « = |, or m = »». or « = ?, where »», f 
are the elliptic coordinates of the point whose Cartesian 
coordinates are x, y, z. 

Hence, whatever he u, we have 


{a’‘+u){b^ + u)(c*+u)f -^ — + -j^ + -r l) 

' ' '\a^ + u + M C* + M / 

= ir,(g 

and therefore, 

(a* + oi)(6* + ai)(c’ + a,) f 1'] 

+ ai 0® + ai c® + ai J 

= jr,(5-a.) 


where Ki and denote quantities independent of the 
coordinates. 

Hence, if wo denote the right-hand member of the 
equation (' ' Q, we have 

Q = C{P . (5 - a„) (17 - ai)(jj -a,) ... (ti- a„) 

(? - «,)(? - a,)... (H- an) = 


issions for », y, and z given by (6) we 
lat when 


Md — ^ V. ^ 

CEHH' = v/ {(a* - - c®J ) 

^=v'{(«’ + 5)(i* + 5)}S» 

C'.SS’J?' = («> - i®) v' { (i® - c®) (c* - «>) ) «yQ 

£ = v/{(«»+5)(6»+|)(c®+?))ffl„ 

C'Sfl'H' = (a® - 6®) («® - c*)(6® - c®)y ^ 


( 12 ) 

(13) 

(14) 


We have seen that there are 4v + 1 ellipsoidal har- 
monics of degree v in that is, of degree 2v in Cartesian 
coordinates. 


F2 
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A rational and integral function of a?, y, z, of the 
degree ^ contains 


0-+l)(i4-2)(/ f 3) 
6 


constants ; 


but if it satisfies Laplace’s equation, these constants must 
satisfy 


equations, 


and therefore such a function contains only (i + ly indepen- 
dent constants. Now, if we take all the different ellipsoidal 
harmonics from the degree 0 up to the degree i in a?, s, 
or \i in 5, we have 1 + 3 + 6 . . . . + 2^ + 1 in all ; the sum 
of this series is {i + 1)*. 

Hence, as each harmonic may be multiplied by an 
arbitrary constant, we can express any rational and integral 
function of y, z of the degree /, which satisfies Laplace^s 
equation, by a series of ellipsoidal harmonics, whose degrees 
in «?, y, z range from i to 0. 

At the surface of an ellipsoid of the confocal system any 
rational and integral function of a?, y, z can be expressed as 
a series of ellipsoidal harmonics. 

For, if J, c be the semi-axes of the ellipsoid, at its 
surface 

a?* y^ z^ 

9 O > 


and therefore, by multiplication, the degree of any function 
of 0 ?, y, z can be increased by 2 without altering its value. 
Hence a rational and integral function of the degree i can 
be reduced to two homogeneous functions of the degrees i 

andt-1. Of these, the first contains — — ^ ^ inde- 


pendent constants, and the second 


*(*• + 1 ) 
2 


Hence the two 


together contain (t + 1)* independent constants, and can 
therefore be expressed as a series of ellipsoidal hairmonios. 
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171. liurface Integral of Product of Harmonics. — 

If Ti = PiHiH\, and ^ V, = and S and S' 

denote two confooal ellipsoids of the system, whose normals 
drawn into the space between them are v and v', by Green’s 
theorem, we have 


I '.S' . |f. *\vfpS’. 

also where p is the central perpendicular on the 

di 

tangent plane to S, and dv - - and therefore 


2/ 


i( 




that is, 


] R,E\H,H\pd8 


At corresponding points on the surfaces S and S' the 
coordinates ij and 2 are the same, and therefore, so also are 
the values of H^ H\, J7a, H' 2 , Also, by Ex. 7, Art. 90, 
the volume elements pdS and p'dS' are proportional to the 
products of the semi-axes of S and S', that is, to J and 3E'. 
dE 

Henoe, as E and ^ are constant over the surface 8, we have 
dEA 


1( 


r 

1 


“ *■ -“I ~ ~ -“.-TP- 






EiEiE\H'tpd8 = 0. 


If we equate to zero the first factor of the left-hand side 
of this equation, we get 

Wrf(j).3ew(j!). 
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Since one surface may be taken as fixed and the other as 
variable, this equation is equivalent to 

(^) 

whence E 2 = C'Ei. 

Accordingly, either E2 and E^ differ only by a constant 
factor, or E2 is the external harmonic corresponding to tlie 
internal Ei. In either case Hi is the same as Hz^ and H\ 
as H\. If we reject the alternatives considered above, we 
must have 

IHiHzH'iH'zpdS^O, 

and therefore, 

iViVzpdS^O. 

Hence we conclude that the surface integral of the pro- 
duct of two ellipsoidal harmonics and the central perpendicular 
on the tangent plane, taken over an ellipsoid of the confocal 
system, is %ero, unless the two harmonics have a constant ratio 
to each other, or be a corresponding pair of harmonics, one 
internal and the other external. 

172. Identity of Terms In equal Series. — If two 

series of internal or of external harmonics he equal to each 
other, each harmonic of one series must be identical with a 
corresponding harmonic of the other. 

To prove this, let the series 

Fo + + Fa ... Fn = ; 

multiply each side of this equation by V^pdS, and integrate 
over the surface 8, then all the integrals on the left-hand 
side vanish except J Vm^pdS, and on the right-hand side 
they all vanish, unless Vm = CVmy in which case we have 

iU„,Vn,pd8^ ClV^fpdS. 

Hence C -1, and a harmonic of the right-hand series is 
identical with V^- 
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If two series of harmonics be equal throughout the whole 
of the space inside or outside an ellipsoid, both series must 
be composed of harmonics of the same kind, either internal 
or external, since an internal harmonic becomes infinite in 
external space at an infinite distance from the centre, and 
the differential coeflScient of an external harmopic becomes 
infinite at the focal ellipse in the plane of xy. 

It is easy to show by multiplication and integration over 
the surface S that, if two series of harmonics be equal at the 
surface of the ellipsoid /S whose semi-axes are a, h, c, and one 
series be composed of internal harmonics, the other of external ; 
then, if a term Vm of the first series be given by the equation 

' m m -*-*■ m » 


there must be a term Um in the second such that 


where 









0 


dl 


= 1 . 


173. Density of liurface Distribution producing 
given Potential. — If V denote the potential inside the 
ellipsoid a, 6, c of a distribution of mass on its surface, and 
V' the potential in external space of the same distribution, 
and if 


f" 


we have seen, Art. 169, that 

V' = CEjm' 


dj 


In this case, if a denote the density of the distribution, 
we have 

dV dV' , 



rfv = - 


2p’ 


but 
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whence 

2jr<r 

~F 


that is, 


dV 

dl 


dV' 

dS 


— fi: 


d^ 

d^ 

lE^ 


+ C 


HR' 


27rff = 


C 

abc {E)^ 


pRH'. 


(15) 


If the potential due to the surface distribution whose 
density is <r be the sum of a number of harmonics Fo, Fi, Fj, 
&c., it may be shown in a similar manner that 


where 


P (“l)o (^2)0 




F, = CxEiH.H 


' f” ^ 

‘Jo W 


When a is assigned, (16) enables us to determine the 
funotions ffi, B.\y &c.j and from thence Fi, &o. 


174. Potential of Homceold. — As an example of the 
mode of procedure described in the preceding article, we 
may find the potential of a homoeoid. Here o* varies as p, 


and - is constant, whence 
P 





dK 

$ ' 


To determine we have 


Hence, 


dV' ^ dV' 2pCo 

Oft 

^JpdS^jNdS^^irM, 


where M is the total mass of the homoeoid, but jpdS = 47r^/6c; 
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and therefore 2(7o= accordingly, V\ the potential of the 
homoeoid in external space, is given by the equation^ 

a J, I 

This agrees with the result found in Ex. 3, Art. 75. 

175. Ellipsoidal Harmonics of the S^econd Degree 
In the Coordinates. — The forms of jB which correspond to 
functions of the second degree in the coordinates are 

and (? - a). 


There are two functions of the last form, 5-ai, and 
5-02. We proceed to determine the values of oj and 02 . 

By (6), Art. 168, we have 
+/)(? -a) = {^ ^ - aj 

= f 5^ -t- (a* + 6* + c’*) 5 + ; 

whence 

m = -I, y - tna = a* + -/o = J (a^b^ + ¥c^ + c^ar). 


Eliminating m and/, we obtain 

3a* + 2(a* + 6* + c*)a + + ¥c^ + £?*«-» = 0. (17) 


Hence, 




- (a* -f ¥ + c*) + */{a^ + 6* + c* - a*6* - - c*a*) 

3 

3 


( 18 ) 


a% « 
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177. Potential of a Poealoid. — We have seen, 
Art. 83, that, for a focaloid distribution of mass, the surface 
density <t is given by the equation 


where 


^TTa =* 


2K 


P 




rt*6*+ product its value 

d,n,ity of Iho 

hyp. ByAr^/i 

space is ^ ^ ^ j 2 )j 


^ X v/ { (6* + f ) {b^ + n) (b^ + Z) (e* + 5) (0* + ij) (c* + 2) ) ; 

and the remaining two products are reduced to expressions of 
a similar kind. 

We may next assume 


aisc^ + a^y^ + 

= Oo + A'i($ - ~ ai)(2 - ai) + (^ - a2)(»? - agX^ - 


where ai and 02 have the values given by (18). In this 
manne; * ^et 

Of procev 

find the po* = Oq + (7i(«*+ ai)(i* + ai)(c* + oi)' 

^ ' ‘-‘•^.stant, ( ^ ) 

^ 7T7. " 

',= I + C',(«* + a2)(J* + ns)(c® + a») j 


X 


a;* t/* 2 * . 

1 L — 1 

a* +02 + 02 C* + O 2 


Hence we have 


0^ - (?i(a*+oi)(J*+di)(6**+0i)— (72(a*+a2)(J®+02)(c*+a2) *0, " 
Oi s= 0\{b^ + oi) (c* + Oi) + C%(b^ + 02)(<'* + 02 )) 
at « i7,(c* + ai)(a* + oi) + Ct{o^ + 02 )(a* + 02 ). 


( 19 ) 
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By substituting the equivalent expressions in and % 
for the functions 

(S - ai)(i? - ai)(^ - ai) and (g - (n - a^) (K - az), 
the potential U becomes of the form 

Uq + + fhy^ + UzZ^y 

where &o., are constants. 

.y _ aid^ - a 2 b^ + of Attraction of Focaloid and 

6*) (c» + ai) (ai - 02 ) ’ the form of U it follows 

''I space due to a 

At the surface of the ellipsoid «, b, c, the f<>roe 
function aix^ + azy^ + OsZ^ can be reduced to ellipse 
harmonics, whatever be the values of ai, az, and « 3 . Fo*^ 
in this case, we have 


x^ 

a* ^ 6* ^ c* 


-1 = 0 ; 


and therefore we may substitute for the given function the 
expression 


Ui + 


* (‘ 


+ +l«»+7f)2 -*• 




We can then determine e so as to satisfy the equation 


«l + ^ + «2 + ^ + «3 + ^ = 0, 


and putting fli + 4 
a* 


for ai, and 



for az, proceed in 


the same manner as before. In this case, the right-hand 
side of the first of the equations (19) is not zero but - e, 
so that 

Cq = ’Ci{(i^ + oi)(J* + + a2)(6* -i 02) (c* + 02) — 6. 
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In like manner, 

r* rfP f" dt. 

+ ?)' 

At a great distance from the centre, X, and iZ’ tend 

towards the values 

MM . M 

-rr, and 

where M denotes the mass of the ellipsoid. 

If we expand in descending powers of 5 the functions 
nnder the integral sign in the expressions for X, X, and X, 
and integrate, we find that X tends towards 
that Y tends towards and Z towards 

Hence we have a = = y = and we get 

179. Potential of KUIpsold in External Space. — 

If V' denote tke potential of the ellipsoid in external space, 

V' = - j [Xdx + Ydy + Zrfz) + constant. 
Integrating by parts, we find 

— d- 

f , r d^ ^‘r dK da> ^ 

(«* + ^j2e 2jf(«» + 5)ai J2(a“ + ?)I 

" 2 ]^ {a^ + m^i 2 («* + ?)$’ 

where in the second integral x is to be regarded as a function 
of p, and s ; and p and s are to be looked upon as constant 
in the integrt^tion. 


( 22 ) 
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Similar results hold good for the functions contained in Y 
and Z, Hence we have 


r 




, - a* Wg 
+ g c»+?/ as 


No constant is to be added, since the right-hand side of 
this equation vanishes at infinity. The three integrations 
involved in the last integral on the right-hand side are to be 
performed on three different hypotheses ; but, as 








= 1 , 


we have, finally. 





dl 1 





(23) 


We see that this result agrees with (11), Art. 87, remembering 
that in that equation cc, z are not regarded as functions of u. 


180. Potential of Plllpsold in its Interior.— The 

components of force due to a focaloid at the outside of its 
surface, by (22), are 

that is, if we put 

they are Aid, By, and Cz. 
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Hence, at any point inside the fooaloid, if these components 
be denoted by X, Y, and Z, we have 



For these expressions satisfy Laplace’s equation throughout 
tlie interior of the fooaloid, and take the proper values at its 
inner surface. 

Hence the potential U of the fooaloid is given by the 
equation 

-^(^Ax^ + By^+Cz^) + + A,, 


where Aq denotes an undetermined constant. 

By Art. 83, if V denote the potential of the ellipsoid 
inside itself, 




t. 



Hence, 




At the surface of the ellipsoid, V' = V, and therefore, 
by (23), we have 


ir+ Ao = 

rd^ 

lo 

accordingly, 



d^ ..r d^ j 

o(6’ + ^)3e "Jo(o» + £)*) 


( 24 ) 



.ungnet of Finite Dimensions. 


81 


CHAPTEE IX. 


MAGNETIZKD BODIES. 


Section I. — Constitution and Action of Magnets. 


181. magnet of Finite Dimensions. — A magnetized 
body is composed of elements each of which is a magnetic 
particle (Art. 17). When such a body is placed in a 
uniform field of magnetic force, each particle is acted on 
by a couple, and the resultant of all these couples tends to 
bring the body into a position in which a certain line in the 
body is in the direction of the uniform force. When the 
body is in the position in which this couple is the greatest 
possible, its ratio to the force is the magnetic moment of the 
hody. 

If a body be composed of a number of infinitely thin 
parallel bars, magnetized at their extremities so that the 
pole strength of each bar is proportional to its orthogonal 
section, it is plain that the magnetic moment of the body 
is proportional to the sum of the products obtained by 
multiplying the length of each bar by the area of its 
orthogonal section — in other words, to the volume of the 
body. Hence we may assume that the magnetic moment 
of an element of a magnetized body is proportional to the 
volume of the element, and we may denote this magnetic 
moment by the expression Jr/©, where «?© denotes the 
volume of the element, and I the intensity of magnetization. 
This latter is defined as the ratio of the magnetic moment of 
the element to its volume. Magnetization is a directed quantity, 
and its direction is that of the parallel bar magnets which 
are regarded as composing the element whose magnetic axis 
is a line in this direction. 


G 
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182. Potential of Magnetized Bot^Vrr% (^8), 
Art. 64, the potential of one element of the body is 

I cos 

o 

r> 


lix, z denote the coordinates of the element ; r/, ^ those 

of the point at which the potential is required ; X, jm, v the 
direction-cosines of the magnetic axis of we have 


cose = 


\ dr dr 
^ dx^ ^ df/ 


+ 



The quantities IX, Ifi, and Iv are termed the components of 
magnetization, and may be denoted hy A, B, and C. If V 
denote the potential of the magnetized body, we have, then, 


V 





If I, m, ?. direotion-coRW o# normal to 
the surface 8, ■which is the boundary of the body, we get, by 
integration, 


V 



+ niB + mC) 


r 


dA dB dC\ d® ^ 
dx dy dzjx 


( 2 ) 


Hence the potential of a magnetized body is equivalent 
to that of a volume-distribution, tliroughout the space occupied 
by the body, of mass wliose density is 


dx ^ dy^ dz )* 


together with a distribution on the surface bounding the 
body whose density is 


lA + mB + nC. 
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183. Poisson’s Pquatton. — From the expression for 
the potential given by (2) we have 


V^r-47r 


'dA dB 
^ dx dy 




d% ) 


0 


( 3 ) 


at any point inside the body. 

In space outside the body Laplace’s equation obviously 
holds good. 

184. Examples of Magnetized Bodies. — As an 

example of a magnetized body, we may take a sphere 
magnetized in a uniform direction so that the magneti- 
zation at any point is a function of its distance from the 
centre. Here, if / denote the distance from the centre of 
any point of the sphere whose coordinates are of, y\ z'; 
F the magnetic potential of the sphere at an external point 
whose coordinates are a?, y, s ; and r the distance between the 
points xy% and we have, the direction of magnetization 
being parallel to 


/<-■’> ;s(;) 7, 




where Jf denotes the component of force due to a sphere 
whose density at any point is /(r') Hence 


F = 



r^dr\ 


where a denotes the radius of the sphere, and therefore the 
magnetic action of the sphere at an external point is the 
same as that of a small magnet at the centre whose magnetic 
moment is expressed by 

4 ^ 1 “ r'^S{r’)dv'. 


If the magnetization be of uniform intensity J, the magnetic 
moment becomes 


47rJ 


G2 
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non 


— (28) 

If an ellipsoid be uniformly magnetized in the dire . . 
of its longest axis, the potential F, at an external point xgz^ 
is given by the equation 


V^l 


{±(1 

Jffe'Vr 




where X denotes the component of force of a solid homo- 
geneous ellipsoid whose density is 1. Hence by (12), Art. 87, 
we have 


V = 27rlnbcx 


du 

(a^ + u)^ (Jj^ + n)^ [c^ + ’ 


(4) 


where q is the greatest root of the equation 

// 

j_ — 4- — . — 

q Ir-^ q c? q 

(See Ex. 1, Art 52.) 

If the integral in (4) be denoted by I, and the corre- 
sponding integrals for the other two axes by and 3> it is 
easy to see that the potential F of an ellipsoid, uniformly 
magnetized in a direction inclined to the axes at angles 
whose cosines are X, /a, v, is given by the equation 

F = 2nlahc{kxlli + juygj + v)53)* (5) 

From (5) we can obtain the components of the magnetic 
force exerted by the ellipsoid at an external point. By 
differentiation we have 

-j— = 2 it lahc I X(|f + \x -y- + ug -7 — 1- vz ] ; 
dx \ dx dx dx J 

but 

^ 1 

dx (a* + (6’ + q)i (c* + q)^ dx 
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If we denote the semi-axes of the ellipsoid passing through 
the point xyz and conf ocal with the given ellipsoid by a', c\ 
we Lave 


^ 4 . - 0 

a'* \a'^ b'^ c'^J dx 


Hence, if y denote the central perpendicular on the tangent 
plane to the ellipsoid a'6V at the point xyz, we have 


da; <1 “ 


Hence . 

dx 

In like manner, 

^ 3 . 

dx dx (Pb'd'^* 


accordingly, by substitution, we obtain 


o r 7 ^ ^ 47rJr?Jc px 

— = 27rIabcA^ r.rr ”77 

dx abc a* 



pp'y vp'z\ 
b'* c'V’ 


If OT,, ^2, CT3, and 3^ denote the angles which the normal 
to the ellipsoid a'b'c' at the point xyz makes with the axes, 
and with the direction of magnetization, we have 


p^' n 

-7j = cos 7xji, &0. ; 

Cl 


and cos & = A cos wj + ju cos wj + v cos Wj. 

Hence Tve obtain 

n r 1 ^ If ^irlabc . 

-7- = 2vlabc\s TTF-r cos 3 ^ cos oti ; 

dx abc 
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and, since the components of magnetic force a, j3, y are 
expressed by & 0 ., we have 


a = cos & cos OTi - 27rIabcX^\ 1 

abc ' 

Q = cos S cos ZJ 2 - 27r/(76c?ju§)', 

' a oc 

y = coaS COSOT3 - 2TrInhcv^\ I 

‘abc J 


( 6 ) 


From equations (6) it appears that tlie force exerted by 
the magnetised ellipsoid abc at an external point P is the 
resultant of two forces of which one is in the direction of the 
normal at P to the ellipsoid a'b'c\ and is expressed by 


4:7rlabc 

“W7“ 


COS&; 


the other is the force due to a homogeneous solid ellipsoid, 
coinciding with abc, at the point Q in which a line drawn 
from the centre in the direction of the magnetization meets 
the surface of the ellipsoid a'b'c\ the density of the solid 


ellipsoid, supposed attractive, being where R denotes the 


distance of Q from the centre. 


185. Potential of IHagnetlsEed Body expressed as 
lium of Force Components. — Adopting the notation of 
Art. 182, we have, by (1), 


r drij r d^] r 


(7) 


Hence, if we suppose three bodies geometrically identical 
with the magnetized body, and having for densities A, P, 
and C, the magnetic potential is equal to the sum of the 
force components exercised by the first body parallel to tlie 
axis of Xy by the second parallel to the axis of y, and by the 
third parallel to the axis of s. 
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186. JMIagnetlc Force. — The differential coefficients of 
the potential with their signs changed are termed the com- 
ponents of magnetic force. Outside the magnetized body 
these are the actual components of the force which the body 
would exert on a north magnetic pole of unit intensity. 

Inside the body the actual force due^ to the body is 
indeterminate. In order to imagine that such a force should 
act, we must suppose a small cavity inside the body, and, in 
the case of a magnetized body, the force depends on the 
shape of this cavity. 

The components of the magnetic force are usually denoted 
by the letters a, /3, 7. 

It is easy to see that the normal component of the 
magnetic force as defined above is not continuous when we 
pass from the outside to the inside of the magnetized body. 
This follows from the consideration that the normal compo- 
nent of that part of the force due to the surface distribution 
I A + mB nC is diminished by 4:ir{lA + mB + nC). 

187. Magnetic Induction. — We can obtain a vector 
quantity whose components satisfy the solenoidal condition, 
and whose normal component at the boundary of the magnet 
is continuous, by adding to each component of magnetic force 
the corresponding component of magnetization multiplied by 
47r. This vector quantity is termed the magnetic induction^ 
and its components are usually denoted by the letters 6, c. 
We have, then, 

a ~ a + 4:7rA, Z/ = /3 + 47ri?, 0 = 7 + 47rC'. (8) 


Outside the magnet a = o, 6 = j3, c = 7 ; and 


da db 


dc 


= ~ V’-r=o. 


Inside the magnet 


da db dc 
d^ ^ dri ^ dZ 


= ~ V^F’+47r 



clB dC\ 
^ dij ^ dz)^ 


and, by (3), the right-hand member of this equation is zero ; 
accordingly, a, b, and c always fulfil the solenoidal condition. 
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At the siu:face bounding the magnet, in passing from a 
point outside to a point inside, la + + nj is diminished 

by 47r {lA + mB + nC ) ; but outside the surface la + mb + no 
is the same as la + mjS + wy, and inside the former exceeds 
the latter by 47r {I A + mB + nC). 

Hence, in passing through the surface, the value of 
la-^-mb + nc is unaltered. 

It is now easy to see that the surface integral 

I {la + mb + nc) dS 

taken over any closed surface is zero. 

If the surface be altogether outside or altogether inside 
the magnetic body, this follows from taking the volume 
integral of 

da db dc 

If tlie surface 8 be partly outside and partly inside the 
magnet, the enclosed volume is divided into two parts by 
the intercepted portion of the surface of the magnet. Tlirough 
each of these parts the integration may be effected, and in 
consequence of the continuity of the normal component of 
magnetic induction, the two surface integrals which are taken 
over the portion of the magnet surface are equal in magnitude 
and opposite in algebraical sign, and therefore the surface 
integral of induction over the closed surface 8 is zero. 

188. magnetic Force nnd magnetic Induction 
regarded as Forces. — If we imagine a small cylindrical 
cavity whose axis is in the direction of magnetization, and a 
north magnetic pole of unit intensity placed at the middle 
point of this axis, the actual force acting on this pole is the 
magnetic force when the cylinder is long and narrow, and 
the magnetic induction when the cylinder is sliort and broad. 

As the cavity is supposed to be small, the removal of the 
volume distribution with which it was occupied produces no 
sensible change in the force acting on the magnet-pole, and 
this force is therefore due to the volume distribution through- 
out the magnet, the surface distribution on its boundary, and 
to the surface distribution on the surface bounding the cavity. 



Energy due to Magnet. 


89 


In the case of a cylinder parallel to the magnetization axis, 
I A + mB ^-nC is zero except at the plane ends, where it is - / 
at the positive end, and + / at the negative. By (3), Art. (14), 
the force due to the surface-distribution is, tlierefore, 


47r/ 


( 


in the direction of magnetization, where c denotes the semi- 
axis of the cylinder, and a its radius. 

When c is large compared with a, this expression becomes 
zero ; and when a is large compared with c it becomes 47 r/. 
Hence, in the first case, the components of the total force 
acting on the magnet -pole are a, j 3 , 7 ; and in the second 
a + 47r^, )3 + 47rJ5, 7 + 

189. £nergy due to magnet. — When a magnet is 
placed in an independent field of force, if V denote the 
potential of the field at any point where there is a south 
pole of strength SD?, the energy due to the presence of this 
pole is - 5 KF, and that due to the corresponding north 
pole is 


■where dh is the axis of tlie particle whose poles are 3)1 and 
- 3)1. Hence the energy due to the particle is 


m%dh. 

dh 


If A, ;u, V be the direction-cosines of dh, we have 


dV dV dV dV^ 
dh~^ dx"^ dy dz’ 


also, = /c/@, and tlierefore 
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Consequently, if W denote the energy due to the presence 
of a magnet in an independent field of force, 


<*) 

where -4, jB, (7 denote the components of magnetization of the 
magnet at any point where the potential of the field is F*. 

190. Knergy of magnetic liystem. — When the field 
of force is due to the magnets which are present, it is plain 
that if the magnetization be everywhere increased in the 
same ratio, the potential is likewise increased in this ratio. 
Hence, by reasoning similar to that employed in Art. 50, we 
see that, if IF denote the energy of a magnetic system, and F 
its potential at any point, we have 



( 10 ) 


If we integrate by parts the expression for W given 
by (10), we get 

W= V{lA + mB+nC)dS-\^V\^-^+^-^ + '^yi<S> 






( 11 ) 


where the last two integrals are taken throughout the whole 
of space. 


191. Vector Potential of magnetic Induction. — 

We have seen. Art. 187, that &, c, the components of mag- 
voluikLinduction, fulfil the solenoidal condition throughout tlie 
sensible cl space, and that the surface integral of induction 
this force is \osed surface is zero. From hence it follows that 
out the magnet, t®fi:ral has the same value for any two surface 
to the surface distribuimu boundary. 
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Hence the integral of induction taken over a surface- 
slieet 8 must be expressible as a line integral taken round 
the curve 5 which is the boundary of 8. We have, therefore, 
an equation of the form 




The directed magnitude of which jP, (t, S are the components 
is called the vector potential of magnetic induction, 

192. litokes^is Theorem. — If v, w denote three func- 
tions of the coordinates, Stokes’s theorem is expressed by the 
equation 



where /8 is a surface-sheet, and s the curve which forms its 
boundary. 

To prove this, we observe that the terms in the surface 
integral which contain w may be written 



m dw\ 
I dx ) 


MS, 


or 



m dw' 
I dx , 


dy dz. 


In this double integral x is regarded as a function of y and z 
given by the equation of the surface 8, 

If ~ w; denote the differential coefficient of tv taken on 

dy 

this hypothesis, we have 


d dw dw dx 

dy dy dx dy ’ 


but the differential equation of the surface may be written 
Idx m dy •¥ ndz 0, 


and therefore 


dx m 

dy I 
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Hence 


d dw m dw 

dy dy I d,v * 


ii 

If the terms containing ti and v in the double integral be 
treated in a similar manner, we obtain the right-hand member 
of (12). If the axes be drawn in the usual manner so that 
counter-clockwise rotations round a?, y, and z bring y to s, 
z to iT, and x to y, equation (12) shows that the direction of 
integration round s is counter-clockwise as viewed from the 
positive end of the normal to S. 

If the surface-sheet 8 be contained between two curves^ 
the surface integral is equal to the difference between twa 
line integrals. 

193. Determiaation of ITeetor Potential. — It follows 
from Stokes’s theorem that as a consequence of (12) we may 
assume 


a = 

dH 

dG \ 

dn ' 

dK’ 

I. 

dF 

dH 

0 = 

dK ' 

■ rfE’l 


dG 

dF 

c = 

dl ■ 

' dn\ 


' dz , 
W -r* ds, 
I ds 


where ij, ^ denote the current coordinates. If jFi, Oiy Hi 
be three functions of 5, rj, ^ satisfying these equations, it ia 
plain that they will be satisfied also by 


F 4. llr 


dn 


Hi + 


dt' 


where ^ is any function of the coordinates. 

We see, then, that equations (13) are not sufficient to 
determine F, 6, and H, and we may assume 


dF dG dH 
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From (13) we have 

_ i'L - '19 - ('^9 

dn~ 'di\dt V </»/V 

= £('!E '9 

"rf,, " <y 


whence 

since 


V-H + 47r 




dfi da 
dl ~ Tii 


V^H-, 

(16) 


Equation (16) is similar in form to that for determining 
the potential of an attracting mass. Hence apparently we 
have 


II- 



(U\ (je 

dy) r ’ 


the integral being taken through the whole of the magnetized 
body. This integral is, however, indeterminate, as at the 
surface J?, and 0 are discontinuous, and their differential 
coefficients in the direction of the normal infinite. 

If we integrate by parts inside the boundary of the 
magnet, we get for H the expression 


We may therefore assume 


F = 
G = 

£ = 



C 


-A 


'>U. 

±\l 

dzj V 


-de 

r 


de 


B 


dx. 


r 






(16) 


provided these forms satisfy the differential equations (13) 
and (14). 
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It is easy to see that this is the case, for since 


di dv 


and A, B, C are not functions of ?/, we have 
dt! dZ, 





+ 

ay az 


-cm - V2 


J r 


( 17 ) 


but, by (1), the first term iu the right-hand member of (17) 
is a, and in space outside the magnet the remaining term is 
zero, and inside the magnet, when coincides with ir//s, its 
value is ^ttA, Hence we obtain 


(IH^ (m 

(iri dX, 


a + ^ttA - a. 


(18) 


We may write 


-jr in the form 


JL 

dKdri^ 



d^ 

dm. 



^ . dG 

and expressing — 

(Iri 

that 


and in a similar manner. 


^ (m m 

^ drt ^ dl 


we see 


vanishes identically. 
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194. Sector Potential of Magnetic Particle. — In 

the case of a magnetic particle equations (16) become 






where SW denotes the magnetic moment of the particle, and 
X, V the direction-cosines of its axis. 

If 01, 08, 03 denote the direction-cosines of r, we have 


whence 


^ 1 ^ Lli = ^ 

dz r r* r*** 




/i03 - V02 


&c. ; 


but iu 03 - V 02 = sin £ cos 

where denotes the angle which a perpendicular to r and 
the magnetic axis makes with the axis of Xy and e the angle 
between r and the magnetic axis. Hence 


F.S«4a-'oo.a„ 


m sin € 

6 r = COS ^25 

r 


^ u;c sin £ 

M ; cos 


Accordingly, the magnitude of the vector potential of a 

magnetic particle at any point is and its direction 

is perpendicular to the axis of the magnet and the line join- 
ing its centre to the point. 

If rotations from x to y, from y to z, and from z to x^ 
be counter-clockwise, the rotation from the magnetic axis to 
radius vector is counter-clockwise as viewed from the positive 
end of the vector potential. The vector potential of a magnet 
of finite size is the resultant of the vector potentials of the 
magnetic particles of which it is composed. 

195. Magnetic Moment and Axis of Magnet. — 
The potential energy of a magnet placed in a uniform field 
of force is determined from (9) by regarding o, j3, y, the 
components of force in the field, as constants ; we have then 


^ W ^ aiAd®-\- Bd® ^ y SCd@. 


(19) 
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If we assume 

J = Kl, J = Km^ J Cd^ = Kn^ - 1, 

then /, m, n are the direction-cosines of a line, and we have 


W = -KH cos 0, (20) 


where H denotes the resultant uniform force, and 0 the angle 
between its direction and that specified by /, m, n. This 
latter direction is fixed in the magnet, and tlie direction of 
is fixed ill space. Hence the magnet is acted on by a couple 

djjr 

expressed by - that is, - KS sin 0, which tends to 

diminish 0 and make the line /, m, n coincide with the 
direction of H. 

Accordingly, the magnetic moment of the body, Art. 181, 
is expressed by A, and /, w, n are the direction-cosines of 
the magnetic axis. 

If the potential of a magnet be expanded in a series of 
harmonics so that at an external point P, we have 


F = 



where r denotes the distance of P from the origin, the first 
Y 

term vanishes, since tlie total magnetic mass is zero, and 

Y 

in the second term, the spherical harmonic Yi is - iTcos 0. 

This is easily seen if we consider the expression for the 
potential energy W due to the presence of a mass at the 
point P. In this case W is given by the equation 


W=mV = '^+&o. 
/•* 


If we now suppose r to become infinite, but to be finite 

and equal to if, we have W = HY^^ but as the energy is 
that due to the presence of the magnet in a uniform field 
of force whose intensity is if, we have JF = - HK cos 0. 


Hence Fi = - A cos 0, 
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196. magnetic Shell. — A magnetic shell may be defined 
as a surface magnetized at each point in the direction of the normal. 

In this case, the expression for the magnetic moment of 
an element of the body is of the form IdvdS^ where dS 
denotes an element of the surface, and dv an element of its 
normal. The total magnetic moment of such a body is in 
general infinitely small ; but if we suppose Idv finite, this 
moment becomes finite. The quantity Idv is termed the 
strength of the magnetic shelly and may be defined as the ratio 
of the magnetic moment of an element of the surface to its 
area. If we put Idv = e/, then J denotes the strength of the 
magnetic shell. 

When the strength of a magnetic shell is the same at all 
its points, J is constant, and the shell is said to be uniform, 

197. Potential of IJniform magnetic Sbeli. — If r 

denote the distance of an external point P from any point Q 
of the shell, by (28), Art. 54, the potential at P of tlie 

IdS 

element of the shell at Q is cos c, where € denotes tlie 

angle between r and the normal at Q. 

But if do, denotes the solid angle which dS subtends 
at P, we have r^^dQ, = dS cos t. Hence 

JdS cos 6 ^ 

r = JdQ, : 

,.2 

and the potential V of the shell at P is given by the equation 

V = Ja, (21) 

where Q, denotes the solid angle subtended by the shell at P. 

This potential differs in character from those with which 
we have hitherto been concerned, as it is discontinuous at the 
surface of the shell. 

If we regard as positive the side of the shell at which the 
north poles of the elements are situated, or towards which 
they point, the potential at the positive side exceeds that at 
the negative by 47r«7. 

s 
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The solid angle subtended at P by the shell is in general 
the same as that subtended by its bounding curve, but the 
two solid angles differ in some important respects. 

The solid angle subtended by the curve is continuous 
except at the curve itself, and in a circuit embracing the 
cui*ve, by passing through its interior, is cyclic. Each time 
the circuit is completed the value of the solid angle is 
increased by 47r. 

These characteristics of the two solid angles we shall now 
consider. 

The solid angle subtended by the shell at P with its sign 
reversed is the same as Gauss’s integral of the normal com- 
ponent of force emanating from a unit mass at P. The sign 
is reversed, because in Gauss’s integral the positive direction 
of r is from P towards the surface ; but, in the present case, 
the positive direction is from the surface towards P. 

If P be on the positive side of the shell, the lines from P 
to the shell which fall inside the cone standing on the 
bounding curve meet the shell once externally, and possibly 
an even number of times afterwards. Those which fall out- 
side this cone meet the shell twice, or some other even number 
of times : first, externally, and then internally, and therefore 
contribute nothing to the integral representing the solid angle. 
Accordingly, the two solid angles are the same when P is on 
the positive side of the shell, and when P is infinitely near 
the shell, on this side, each may be denoted by Qi. 

When P moves across the surface of the shell from the 
positive to the negative side, the solid angle subtended by 
the bounding curve remains unaltered, but that subtended 
by the shell becomes - 47r. To see the truth of this we 
have only to suppose the closed surface completed of which 
the shell is part. Then, by Art. 26, the solid angle which 
the entire closed surface subtends at P is - 47r ; and it is plain 
that Hi denotes the absolute magnitude of that part of this 
angle which is subtended by the portion of this surface which 
has been added to the shell. Hence the solid angle subtended 
by the shell is - (47r - Hi). 

The solid angle subtended at P by the curve bounding 
the shell is everywhere continuous unless P be infinitely near 
the curve. As P moves about, the variations of the two solid 
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angles are the same except when P is passing through the 
surface of the shell. Hence we may take for the potential 
of the shell at P the expression e/iO, where Q denotes the 
solid angle subtended at P by the curve bounding the shell, 
with the proviso that when P passes through the shell from 
the positive to the negative side, 4e77r must be subtracted 
from the foregoing expression. 

If i/; be a function of the coordinates of a point, and if 

ds taken round a closed circuit be zero for every possible 




ds 


closed circuit, -dj is acyclic, but, if for some circuits , . 

}d8 

taken round the circuit be not zero, \p is cyclic. If a closed 
circuit 8 be such that we can draw a surface 8, of which s is 
the boundary, so tliat at every point of S the function \p has 
differential coefficients ?/, w which are finite and continuous, 
then by Stokes* Theorem, Art. 192, the function xp must be 
acyclic for the circuit s. Again, if a surface fulfilling the 
conditions stated above be bounded by two curves, and $ 2 f 
dxp 


the value of 


ds 


ds taken round the circuit is the same for 


one of these curves as for the other. It is now easy to see 
that 12, the solid angle subtended at P by the curve s 
bounding tlie shell, is acyclic for every circuit which does 
not embrace this curve, passing through its interior. For 
since the differential coefficients of 12 are finite and continuous 
for all positions of P not infinitely near the curve s, this 
follows immediately from what has been said above. 

If we suppose P to start from a point at an infinite 
distance on the positive side of the shell and to move in a 
straight line to a point at an infinite distance on the negative 
side, passing in its course through the interior of the curve «, 
the solid angle 12 passes from 0 to 47r. For if a unit sphere 
be described round P as centre, the edges of the cone having 
its vertex at P and standing on s initially converge to a 
point. As P approaches to s the cone opens out, and the 
area swept out on the sphere by the edges of the cone 
increases. After P passes through the interior of s, this 
area becomes greater than a hemisphere, and finally when P 
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reaches an infinite distance on the negative side of tlie shell, 
the edges of the cone again converge to a point on the sphere 
which is now opposite to that to which they originally con- 
verged. These edges have then swept out the entire sphere 
or 47r, 

We may now suppose P to return to its original position 
along a path on the outside of s, and such that all its points 
are infinitely distant from s. At all these points the diffe- 
rential coeflBcients of Q, are zero ; and hence the value of ^2 is 
47r, when P returns to its original position. It is now easy to 
see, from Stokes’s theorem, that for any circuit passing 
through the interior of the curve and embracing it once 

I ^ rfs must be 47r. 

Hence we conclude that the potential of a magnetic 
shell is expressed by a cyclic function, but that at the surface 
of the shell discontinuity occurs in the potential though not 
in the function. In consequence of this discontinuity the 
principle of the conservation of energy is maintained. 

In fact, if P start from a point 0 on the surface of the 
shell, at the negative side, and travel round the edge of the 
shell till it reaches the point O' on the positive side of the 
shell, opposite and infinitely near to 0, the function Qi in- 
creases by 47r, but in passing through the shell from O' to 0 
the potential of the shell is diminished by 47r. Hence the 
value of the potential at 0 is unchanged by the motion of P 
round the complete circuit, but the value of Q, is increased 
by 47r. 

198. CSnergy due to magnetic lihell. — The energy 
due to a magnetic shell placed in an independent magnetic 
field is given by (9). If m, n denote the direction-cosines 
of the normal to the element dS of the shell, and Jits strength, 
we have 

Ade - JldS, Bd® = JmdS, Cd® = JndS, 

and if a', j3', y denote the components of magnetic force 
due to the field (9) becomes 

W^-JJ(la'+m(i^ + ny')d8. 


( 22 ) 
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199. Knergy due to Two Magnetic lihells. — JK the 

magnetic field he due to a second shell 8' whose compon^ts 
of force are «' j3' y\ the energy W given by (22) repres^^s 
the result of the mutual action of the two shells. 

Since the one shell is outside the other, we may in (22) 
substitute tlie components of induction for those of force and 
for the components of induction we may put tlie expressions 
given by (13), Art. 193. Thus (22) becomes 



where s is the curve bounding the first shell. 

The values of F\ 6\ the components of the vector 
potential of the second shell, are given by (16), Art. 193. 
In this case 



but ^ = J'm'dS\ C'de' = JVd8\ 

and hano. F - xj j»,' iQ) - »' j ,IS'. 

In Stokeses theorem (12), Art, 192, if we make 

1 Cdr' 1 

= = we get r = J'\ ~-ds\ 

r J as r 

where s' denotes the curve bounding the second shell. In 
a similar manner we have 




Substituting in (23) the values obtained for F', Q', and H' 
we get 


W = -JJ' 


f dx dx dy dij 
\(/8 rfs' ^ rfs rfs' 


+ 


dzd£\ 
ds ds) 


- ds ds' 
r 


= _ //'IJ (25) 

where f denotes the angle between the curve elements ds andcfs'. 
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Section I. — Induced Magnetism, 

200. JHagnetlc Induction. — When a body is placed in 
a field of magnetic force, in general its magnetism is altered. 
The magnetism produced by the force is called induced 
magnetism. When the magnetizing force is small, the in- 
duced magnetization is, in general, proportional to and co- 
directional with the total magnetic force acting at the point, 
so that if denote a component of induced magnetization, 
and a the corresponding component of the total magnetic 
force, A 2 ~ ica, where ic is a coefficient depending on the 
nature of the body, and is called the coefficient of induced 
magnetization. 

It is easy to see that A^ and a are quantities of the same 
order, so that ic is a numerical, magnitude, which is positive 
in the case of paramagnetic bodies, and negative in the case 
of diamagnetic. 

If Ai denote the component of that part of the magneti- 
zation which is independent of induction, we have 

A = Ai + Kttf B = Bi + /cj3, C = Oi + Ky. (1) 

201. magnetism due altogether to Ind^tlon.-^ 

If there be no magnetism in the body independent of the 
induction due to the field of force, Ai- Bi- Ci = 0, and 

A = KOj B - ffj3, C = Kj. (2) 

In this case, by (3), Art. 183, we have 


dA dB dC (da dS dy\ , (dA dB dC\ 
dx^ dy^ d% \dx^ dy'^ dz) \dx ^ dy ^ dz / 


whence 


cU (W (W 
dx dy dz 


( 3 ) 


When the components of magnetization fulfil this con- 
dition, the distribution of magnetism is said to be solenoidal, 
and the potential corresponding may be regarded as due to a 
surface distribution of mass whose density is 


lA + mB + nC. 
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202. Component of Induction. — When a body has 
no magnetism independent of that induced by the acting 
force, the components of induction are given by the equations 

rt = (1 + 47r»c)a, 6 = ( I + 47ric) j3, 6* = (1 + 47ric)7. (4) 

If we put 1 + 47r/c = TSy the quantity ts is called by 
Maxwell the specific magnetic inductive capacity, and by 
Thomson the magnetic permeability, and in the case of a 
body magnetically isotropic, having no permanent magnetism 
independent of induction, we have, then, 

a = T^ja, b - 2sjj3, c = wy. (5) 

203. Distribution of Induced Dagnetlsm. — Let U 

denote tlie total magnetic potential, inside the body in which 
the distribution of induced magnetism is to be determined, 
the total potential in the external medium ; then, as the 
distribution of induced magnetism is solenoidal, and there is 
no other magnetism inside the field in which TJ and W are 
to be determined, we have 




also J7 *= ?7' at the surface bounding the magnetized body, 
and since. Art. 187, the normal component of induction is 
continuous, 


dU ,dU' 

rfi/ ® dv' 


= 0 , 


( 6 ) 


wliere ts and zs' denote the coefficients of permeability of the 
body and the external medium, and, v and v the normals 
drawn into them at the separating surface. 

If U' be assigned at the surface bounding the field exter- 
nally, U and JJ' can be determined in only one way so as to 
satisfy the given conditions. Let us suppose that the equations 
could be satisfied by two pairs of functions Z7i, U\ and ZTi, U' 2 , 
and let 
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then, if & be the surface hounding the field externally, we 
have 


f ^ = _ [ j + (^Y + 

av ^ jlv^J \dy) \(hj) 

I f + 1 dS' + j 0'vyrf@' = - j j J+ &o. j d<B'. 


If we multiply the first of these equations by w, the 
second by w', and add, observing that at the surface S' we 
must have = 0, and that ^ at /S, we get 



Since 



0 at 8, and = 0, V>' = 0, 


the left-hand member of this equation is zero. The coefficients 
-BT and w' are always essentially positive, even if k or k' be 
negative. Hence each member of the right-hand side, and 
each of the terms under the integral signs, must vanish 
separately, and therefore 0' = 0, ^ = 0. 

204. External IVIedluin not Magnetic. — If the ex- 
ternal medium is not capable of being magnetized, we have 
k' = 0, w = 1 ; also 

r+Q, F+a'; 

where V denotes the potential of the forces producing the 
induction, and £2 and O! the potentials, inside and outside the 
body, of the induced magnetism. In this case V is supposed 
to be given, and £2' is zero at infinity. 
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205. Anisotropic JHedlum. — When a magnetic medium 
is anisotropic or crystalline, the induced magnetism is not, in 
general, codirectional with the magnetic force ; but the com- 
ponents of induced magnetization are linear functions of the 
components of force, so that we have 

A 
B 
C 


= Kua + Kn(i + K137, \ 

= K21CL + If22^ + K2syn / (7) 

= IC;na + K32fi + »C337. ^ 


By (9), Art. (189), we see that, to increase by Sa the force 
acting in the element the work required is - ASa, and 
therefore we conclude tliat 


da " dfi " dy - 

Hence k-h = K 12 , &c., and (7) become 
A = »cu« + K 12/3 + ki 37, j 
B = fc'igfl + K 22/3 + fC23 7> ( 

C = Kisa + IC23/3 + IC337. ' 


( 8 ) 


( 9 ) 


200. Isotropic Kllipsold in Field of Vniform 
Force. — If an ellipsoid, free from magnetism and sur- 
rounded by a non-magnetio medium, be placed in a field 
of uniform force, the distribution of induced magnetism can 
readily be determined. In fact, we may suppose the ellipsoid 
to be uniformly magnetized in a direction to be determined ; 
and if the conditions of the question can thus be satisfied, 
we know by Art. 203 that we have reached the correct 
solution of the question. 

Let I denote the intensity of the induced magnetization, 
and X, p,y V its direction-cosines; then, by (5), Art. 184, the 
potential V of the induced magnetism is given by the 
equation 

r = I{\Lx + fxMy + vNz\ (10) 
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where, by (17), Art. 22, the constants L, M, N denote 




2Trfihc 


1^00 
J 0 


(hi 


+ ?^)i (6^ + ('o* + h )^ ’ 


two other integrals obtained by interchanging b and e 
with a. Hence, if the components of the uniform force due 
to the field be denoted by -F 2 , the total magnetic 

force a, parallel to the axis of is Fi - I\ L, and we have 

lX = Ka^K(F,-LI\). 

with two similar equations. Accordingly, we get 


(1 + kL)I\ = icF„\ 

(1 -f k3I)Iiul== kF,, j (11) 

(1 + kN)Iv = kF^. j 

The values of /, X, ju, and v obtained from these equations 
satisfy the conditions of the question*. 

207. Anisotroiiic Klll|isoid surrounded by UTon- 
Hagnetle IMIediuni in Uniform ifield of Force. — 

In this case, if we proceed in a manner similar to that of 
the last Article, we get 

/A = Kn {F, - IXL) k-,2 (F, - I^M) + (F, - IvN), 


and two similar equations ; whence we have 


(1 + KllZ/) /A + Ku MIjX-\-Ki^NIv = KiiFi + K 12 F 2 + IC 13 -F 3 , ^ 


K12LIX + ( 1 + KiiM) Ijx + Ki^NIv « KnFi + K21F2 + IC23-P3) r 


( 12 ) 


ki^LIX -f K 23 + ( 1 + 1 C 33 ) N1 V = ici3 F\ + 1 C 23 F 2 + 1 C 33 Fs. J 


Hence /, A, ju, and v are determined. 
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Skction III. — Terrestrial Magnetism. 

208. Karth’s Hagnetic Potential. — The components 
of tlie Earth’s magnetic force at any place can be deter- 
mined by observation. This can be done eitlier by finding 
the time of oscillation of a magnet, free to move in a 
horizontal or in a vertical plane, when disturbed from its 
position of equilibrium, or by arranging a position of 
equilibrium under the combined action of the Eartli and 
magnets whose strength and position are known. The 
investigation of the Earth’s magnetic potential was initiated 
by Gauss. In the British Islands some of the earliest observa- 
tions were carried out by Lloyd in the magnetic observatories 
of Trinity College, Dublin. 

When the Earth’s horizontal force has been determined 
at a sufficient number of places, the question of the existence 
of an acyclic magnetic potential can be investigated. 

If 8 denote any portion of a closed path on the Earth’s 
surface, H the iiorizontal component of magnetic force at 
any point, and 0 the angle which its direction makes with 
that of 8, on the hypothesis that a magnetic potential V 

dV 

exists, wo have Jff cos 0 = - — . Hence, if an acyclic 

as 

magnetic potential exists, J H cos 0 ds taken round the closed 
path is zero. By finding a sufficient number of values of 
J1 and d the numerical value of the integral can be deter- 
mined approximately. In fact, if Si and s^ correspond to 
stations not too far apart, we have 

f S2 

JETcosO ds = ^ {Hi cos 01 + H 2 cos (§2 - Si) 

J SI 

approximately. It is found in this way that JH cos ffds 
taken round a closed path is always zero. 

Hence we conclude that the magnetic action of the 
Earth can be represented by an acyclic potential, and con- 
sequently that electric currents passing from the outer 
atmosphere to the ground cannot be the cause of any part 
of this action. 
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209. liocality of tbe iiourcefit of the J!:artli’8 
magnetic Force. — If the Earth’s magnetic action be due 
to magnetism, or closed electric currents in its interior, the 
magnetic potential V at any point P outside its surface can 
be expanded in descending powers of r, the distance of P 
from the centre of the Earth. The difference between the 
numerical values of V at any two places can be determined 
from observations of the horizontal force. If the magnetic 
action be due to magnets or currents outside, the potential 
at any point nearer to the centre than the nearest of these 
sources of action can be expanded in ascending powers 
of r. 

Hence for a point P close to the Earth's surface at its 
exterior we liave 


r fr 


( 1 ) 


where a denotes the radius of the Earth and Ui and Yi 
spherical harmonics. At the surface of the Earth 


V=^X(U,+ r.), 

Of 


and if Oi denote a coefficient in Z7i, and hi the coefficient of the 
corresponding term in Ei, the coefficient of this term in V is 
(ti + hi. * By taking a sufficient number of numerical values 
of V at known places on the Earth^s surface we can deter- 
mine as many of these coefficients as we please so that 
we may regard ((% -h as known. 

If we now consider the vertical component Z, towards 
the centre, of the Earth's magnetic force, we have 


dV 


2 




( 2 ) 


At the surface (2) becomes 

Qi 


( 3 ) 
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Hence, from the observation of a sufficient number of 
values >f Z we can determine 


iai -(?:+!) hi. 


and consequently ai and hi are each known. 

It is found that ai is always zero, and accordingly we 
conclude that the Earth’s magnetic action is due altogether to 
sources inside its surface, and that V, the potential of the 
Earth’s magnetic action, is given by the equation 




SFi 


a* 


(4) 


210. JEarth’s JHasuetlc Poles. — A magnetic pole is a 
point at which the liorizontal force vanishes. At such a 
point this force changes sign so that at each side of the pole 
the same end of the needle points towards the pole. 

If there be two poles of the same kind on the Earth’s 
surface in going from one to the other along a magnetic 
meridian, the horizontal force must change sign and therefore 
vanish. Hence there must be a third pole between the two 
former. The end of the needle which pointed towards these 
poles points away from the intermediate one at each side. 

As a matter of fact there are only two magnetic poles on 
the Earth’s surface, and these two are of opposite kinds. The 
proximity of these poles to the extremities of the Earth’s 
axis of rotation appears to indicate a connexion between the 
Earth’s magnetism and the Earth’s rotation. From the 
properties of electric currents it is easy to see that sucli 
currents circulating round the Earth, and approximately 
parallel to the equator, would account for the magnetic 
phenomena exhibited. 
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CHAPTER X. 

KLECTRIC CURRENT’S. 

211. Introductory. — Not long after the discovery of 
current electricity it was observed by Oersted that a wire 
through which an electric current is passing exercises an 
attractive or repulsive force upon the pole of a magnet- 
needle. It was found also that wires along which electric 
currents are passing attract or repel one another. 

By a combination of experimental and mathematical 
investigations Ampere succeeded in arriving at the laws 
which regulate the attraction of currents on each other and 
on magnets. 

His original investigations must ever be regarded as 
worthy of the highest admiration, but some of his experi- 
ments, combined with the theoretical developments of other 
phjsioists, enable us to arrive at his results by methods 
shorter and simpler than those employed by him. 

212. Electric Currents. — An electric current may be 
produced in various ways ; but in all cases the maintenance 
of an electric current requires an expenditure of energy 
supplied by an external source. 

The source of energy may be chemical, as when two 
substances unite chemically, or mechanical, such as the 
action of a steam-engine or water-mill, but in all cases there 
must be a source of energy outside the current itself on 
which its continuance depends. The force due to an electric 
current is not therefore a permanent natural force, and pro- 
positions depending on the principle of energy cannot be 
applied to it in the same manner as to gravitation, or to the 
attraction of static electricity. 

.The currents whose attraction we are about to consider 
are those transmitted along a wire of small section. 
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The quantity of electricity which passes through an 
orthogonal section of the wire in the unit of time is called 
the strength of the current. The quantity which passes 
through the unit of area is called its intensity. When a 
steady current is established, the strength is uniform through- 
out the wire. The force which causes and keeps up the 
current is the electric force. When there is a potential 
corresponding to this force, the force is the rate of diminution 
dV 

of the potential, or — where s denotes an element of 

length along the wire. As the current is supposed to be 
constant, this force must be equilibrated by another of equal 
magnitude. 

The current is thus analogous to the uniform motion of a 
body sliding on a rough surface. 

The retarding force on a unit of electricity is found to be 
proportional to the intensity of the current, that is, its 
strength per unit of area. 

Thus we have 

ds a 

where k is a coefficient depending on the material of the 
wire, <r denotes the area of its section, and i the strength of 
the current. If we integrate the equation above, we get 

Vx- 

If k and a be constant, this becomes 


Vx- = lit, 

a 


where I denotes the length of the wire, and E - 


quantity E is termed the resistance of the wire. 

^ If Fi - Fj, the difference between the circuit 

potehH8l . at the extremities of the wire, be valent to a 
equation (I") -may be written to d8 and 

' ■ Pi strength of the 
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This expresses what is called Ohm’s Law. E is termed 
the electromotive force^ and may be defined as the difference in 
potential between the extremities of the wire, or, more generally, 
as the line integral taken along the wire of the electromotive 
intensity. 

The term ‘ electromotive force ’ applied to this integral 
seems highly objectionable, but is sanctioned by long usage. 


213i Solenoids. — If a wire be bent into the form of a 
circle, not quite closed, be carried on for a short distance at 
right angles to the plane of the circle, bent into another circle 
equal and parallel to the first, carried on again, and so on, 
and finally brought back in a straight line perpendicular to 
the planes of the circles and close to the connecting portions 
of wire between them ; and if an electric current be sent 
through the wire, we obtain what is termed a solenoid. As 
the portion of the current which is perpendicular to the 
planes of the circles consists of two parallel parts close 
together and flowing in opposite directions, it produces no 
attraction on a magnet-pole, and the solenoid may bo 
regarded as being composed of a number of equal circles 
whose planes are perpendicular to a straight line passing 
through their centres. 

It is found that at distances which are large compared 
with the diameter of one of the circles, the solenoid exercises 
the same action as a linear magnet. 

If or denote the area of one of the circles, S the perpen- 
dicular distance between two of them, I the length of the 
solenoid, and i the strength of the current, it is found that 
the magnetic moment of the solenoid is expressed by 



muTJie magnetic moment of a linear magnet of equal length, 
which of small magnets having each a magnetic moment p 
current is . . , u , , 

positions de^®“8t“ ^7 ^4 The axial length A 

applied to it i^etween the centres of two of the sma" moo-nets 
attraction of s. ia compoj.«^,„ If t; . 

.The currents v. tV'- ''** ?»?>?• moment oi 

are those transmitted aiv,. *[^i^gnet, and if we suppose 
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As the equivalence of the solenoid to the magnet holds 
good, whatever be the number of circles in the solenoid, we 
conclude that, 

A small circular current is equivalent to a small magnet 
whose centre coincides with that of the circle, whose axis is 
perpendicular to the plane of the circle, and whose moment 
is equal to the area of the circle multiplied by the strength of 
the current. 

The equivalence of a solenoid to a linear magnet holds 
good equally well if another plane curve be substituted for a 
circle, and becomes more rigorously true according as the 
diameter of the curve is diminished, compared with the 
distance of the magnet on which the solenoid acts. Hence 
we conclude that. 

The magnetic action of an infinitely small electric circuit 
is equivalent to that of a magnetic particle whose axis is 
surrounded by the circuit and is perpendicular to its plane, 
and whose magnetic moment is equal to the area of the 
circuit multiplied by the strength of the current. 

214. Equivalence of Electric Circuit to Magnetic 
illiell. — If a single-sheeted surface be described of which an 
electric circuit is the boundary, and a network of lines be 
drawn on this surface dividing it into a number of small 
elements, the electric current is equivalent to a current of 
equal strength circulating in its direction (clockwise or 
counter-clockwise) round each of these elements. This is 
obvious if we remember that along the boundary line 
between two adjacent elements there are two currents in 
opposite directions, one for each element. As these currents 
are equal, they neutralize each other ; and the only current 
which remains uncompensated is that in the outer boundary. 
By increasing the number of lines in the network, the size 
of each element can be diminished without limit. 

From Art. 213 it appears that the electric circuit 
embracing the element dS of the surface is equivalent to a 
magnetic particle whose axis is perpendicular to dS and 
whose moment is idS, where i denotes the strength of the 
current. 

I 
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Hence, the total electric circuit is equivalent to the 
assemblage of small magnets, normal to the surface 8, 
whose moments are tlie areas of the elements surrounding 
the magnets multiplied by the strength of the current, that 
is, to the magnetic shell whose surface is S and whose 
strength is i. 

215. magnetic Potential of Klectric Circuit. — 

Since the magnetic action of an electric circuit is the same 
as that of a magnetic shell bounded by the circuit, the 
magnetic potential of an electric circuit whose strength 
is i at a point P is expressed by iO, where 12 denotes the 
solid angle subtended by the circuit at P. This potential is 
continuous everywhere except at the circuit itself. 

For any closed curve not passing through the space 
surrounded by the circuit the potential is acyclic. 

For a curve passing through this space and embracing 
the circuit the potential is cyclic, and the value of the cyclic 
constant is 47ri. 

These characteristics of the potential show that in moving 
a magnet-pole round a closed curve which does not embrace 
the circuit no work is done, but that in moving the unit 
pole round a curve embracing the circuit and passing 
through its interior, if tlie direction of motion be opposed 
to the force, work is done represented by ^iri. 

If we imagine a person to stand on the positive side of 
a. shell equivalent to the current, that is, on the side towards 
which the north poles point, the current as seen by him will 
circulate counter-clockwise, and if a person is placed lying 
along the current which enters at his feet and goes out at 
Ids head, the motion of a nortli magnetic pole moved by 
the current round his body will as seen by l)im be counter- 
clockwise. 

The first of these statements is deduoible from the 
experiments made on solenoids; the second follows from 
the equivalence of the current to the magnetic shell. 

216. magnetic Force of Current®. — Since an electric 
circuit is equivalent to a magnetic shell, the components of 
force due to the current are in space outside the shell the same 
as a, j3, y, the components of magnetic force due to the shell. 
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Outside the shell a, j3, 7 are the same as a, b, c, the com- 
ponents of induction due to the shell. At tlie sliell a, j3, 7 
are discontinuous, Art. 186; but since the magnetization of 
the shell is normal to its surface, a, A, c are continuous, Art. 
187. The force-coraponents of the current are everywhere 
cvontinuous except at the current itself. Hence we conclude 
that for all space outside the current, the components of its 
magnetic force are expressed by r/, b, c, being tlie same as the 
components of induction of the equivalent magnetic shell. 

217. JEiiergy due to presence of Klectrie Current 
in Independent Magnetic Field. — Let a\ y' denote 
the components of magnetic force ; r/', 6', c' those of induction, 
due to a magnetic shell 8 equivalent to the current; a, |8, 7 
the components of magnetic force; r/, 6, c those of induction 
due to the field @, and A, jB, C the components of its mag- 
netization. Let U denote the energy due to the presence 
of the shell in the field, and W that due to tJie presence of 
the current. 

By Art. 216 and (9), Art. 189, we have 

W^-j(n'A^b'B + c'C)d<B, 

Except at the surface of the shell, a' = a, b' = /3', c' = 7' ; 
but at 8 we have ad® = ad® + 47rild8y where I denotes the 
direction-cosine of the normal to 8, with similar equations 
for b' and c\ 

Hence U - 47ri j {I A + mB + wC) dS. 

Again, by (22), Art. 198, 

CT = - ^ J (/a + m(i + ^^7) dSy 

and therefore 

W ^ — i j {/o + wjS + ny + 4w {lA + mB + nC)}d8 
= - i j (la + mb + nc) d8. 


12 


( 3 ) 
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218. Force-Components of Current expressed as 
Integrals. — If c denote the components of magnetic 

induction due to a shell equivalent to the current, by Art 193, 
and (24), Art. 199, we have 

dH (IG ^ d r d i idy' 
dy dz dy] X (/gj r ’ 

where y\ z' denote the coordinates of a point on the 
circuit, and r the distance between this point and the point 
ir, y, 2 , and the integrals are taken round the entire circuit. 

Q. 1 1 y-y' 

dy r r^ r 

we get 2-2' dz' y - y'\ ds' 

* J \c/6*' r ds' r / r'”^ ’ 

with similar expressions for b and c; and if jF,, F*, denote 
the components of force exercised by a circuit of strength i 
on a magnet-pole of strength situated at the point Xy y, «, 
we have 


fdz' x-x 


ds' r 7 r* 


. [(dx' y-y dij x--x'\ds' 

\dH r ds T J ^ 

219. Force exerted by Flement of Current on 
JNtagnet-Pole. — The components of force given by (4) are 
the sums of tlie components of force contributed by the 
various elements of the circuit. 

Hence, the circuit acts as if the force- components due to 
a single element ds' of a current whose strength is », acting 
on a magnet-pole of strength m, were expressed by 
nnds' fdy' %-z' dz' y- y'\ 

r* r ds' T )' 

imds' fdz' x-x' dx' z-z'\ 

\d8' T" - ds ~ r 

imds' ^dx' y - y' dy' x-cif\^ 
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[I’hat these are tlie actual force-components due to a current 
element is shown at the end of this Article. 

In the above equations, 

dx^ dif ^ dz' 
rfs' * ds'' rfs' 


are the direction-cosines of the current element ^ 4 *', and 

X -X // - y , s - s' 

• — 1- and 


those of r. Hence, if 0 denote the angle between ds' and r, 
and 3^1, ^2, &3 the direction-angles of a perpendicular to their 
plane, the force-components due to the current element are 
expressed by 

im sin Odd ^ Ini sin Odd ^ im sin Odd 

cos Sfi, T cos ^2, 5 cos ^^3. 

y3 


Hence the force which a current element dd of strength / 
exerts on a magnet-pole of strengtii m is perpendicular to 
the plane containing the pole and the current element, and 
tends to make the polo move in a counter-clockwise direction 
round the current element, along which the observer is sup- 
posed to be situated with the current entering at his feet and 
going out at his head. The magnitude of the force is 

im sin Odd 
r* 

This result can be proved directly from the expression for 
the magnetic potential of the circuit. 

If we suppose an element dd of the circuit to be free to 
undergo a displacement under the action of a magnetic pole w, 
the work done by the force in this displacement will be equal 
to the loss of potential energy of the system. 

The potential energy W of the system is denoted by imQ, 
where Q is the solid angle subtended at m by a surface 8 
bounded by the circuit. 

Let dd receive three displacements ; one, 8^, along dd 
itself ; one, 8i|, perpendicular to dd in the plane of dd and r ; 
and one, dZ, perpendicular to the two former. 8? does not 
alter the surface 8 , The displacement 8i| by the motion of dd 
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generates an increment of the surface /S, but the element of 
surface so generated is in a plane containing r ; and, as its 
normal is perpendicular to r, it subtends no solid angle at m. 
The displacement 8^ alters Shy the amount and the 

normal to this element of surface lies in the plane of r and ds', 
and is perpendicular to the latter. Hence, if 0 denote the 
angle between ds' and r, the angle between r and the normal 

is ~ Accordingly, tlie element of solid angle subtended 

at m by the element of surface is 

^u^eds'8l 
r- ’ 

and therefore 8W = m8Q, = 

r 

Hence the force between m and ds' is in the direction of 
the displacement 8?, and is expressed by 

im sin 0ds' 


The direction in which the force exerted by dd on m tends 
to move the latter is in the direction in which the solid angle Q^ 
at m is diminishing. Thus, we arrive at the results already 
stated. 

220. Energy due to mutual action of two Electric 
Circuits. — Since the action of each circuit in space outside 
itself is the same as that of a magnetic shell, if W denote the 
energy due to the mutual action, by (25), Art. 199, we have 

Jr=-n'|^rf6rfs'. (5) 

It is here assumed that the strength of each current is 
maintained constant. 

221. Forces between two Electric Clreuito.~ 

If Jf, Yf Z denote the components of the force acting on a 
current element in consequence of the mutual action between 
the circuits, for any system of small displacements we have 

^{X8x + Y8y ^ Z8z)^- 8W. 
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In order to determine the variation of W we must express 
cos c in terms of r and its differential coefficients. 

If 07, z denote the coordinates of an element, ds of one 
current, and o;-, y\ s' tliose of an element, ds' of the other, 
find r the distance between these elements, remembering that 
O’, 7jy Z are functions of s, and x\ y\ s' of s', and that s and s' 
are independent of each other, we have 


= (;r - x'Y + (// - !/J {z - 

dV /X / /X / /X 


dv tk 
<h (k' ■*" 


(Pv fdx dx' dj/ dy dz dz' 

dsds v/s ds' ds ~ds ds ds 


== - cos f. 


Substituting for cos £ in (5) we get 

f/ dh 


W 


t i 


1 dv dx 
\ds ds' ^ V ds ds' 


ds ds'. 


Tile first term under the integral sign can be integrated 
round either circuit, and, as the circuits are closed, it 
vanishes. 

Hence 

[ 1 / dv dSv dv dSv \ dv dv Sv ) 


8jr = 


v\ds ds' ds' 


T aor \ dv dv or ) , , , 


we iutegrate by parts tbe first two terms, since the 
circuits are closed, the single integrals vanish, and we have 

a /I <Jr\) 

'6-'\r cisjj 




+ 


( 2 dh' 


1 dr dr 
r dsdd r^ ds ds' 


ds' 

Sv ds ds' 


Sv ds ds' 


-■•11 

= i II ^ I ^ I ® I 


( 6 ) 


where 0 and 6' denote the angles which r makes with s and 



120 


Electric Currents. 


Hence 


S(X8it! + 7hj + XSz) 
2 


= - 


•'2 


cos a + ^ COS B cos 0' 


dsds'Br. 


(7) 


Accordingly, tlie forces due to the mutual action of the 
circuits are equivalent to a system of forces acting in the 
lines joining ilie elements of one circuit to those of the otlier. 
If JB denote the magnitude of the force acting in the line 
joining the elements ds and (Is\ by (7), we have 




2ir 


cos e -f - cos 0 cos 6 ' ) ds ds\ 


( 8 ) 


The negative sign sliows that the force between the 
elements is attractive when the currents are both approaching 
the sliortest distance between their lines of direction. 

The magnitude of It was discovered by Ampere. He 
assumed that the direction of the force between two current 
elements is the line joining them. 

In the investigation above, nothing has been assumed ; 
but it has been shown that two closed currents act on eacli 
other as if there were a force It along eacli line joining an 
element of one current to an element of the otlier. 

So far as this investigation goes there may be other forces 
acting between each pair of elements, but these forces must 
be such as to produce no effect on the total action between 
two closed currents. 

If Uu U 2 J and Uz denote three functions of s and s', in 
addition to B acting along r, there might be three forces: 


-77- a.s‘ 


parallel to the axis of ir, 


ds ds' parallel to that of and 
ds ds' parallel to that of z, 


due to the action of ds' on ds. 
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In this c, ^ , 
on dSy resulti. 


taken round the closed circuit d j M 1 XX 14 . 

As the expression for the force between two eie**. 
must be symmetrical with respect to these elements, ti 
force exercised by ds on dn' parallel to the axis of x would, 
ill this case, be 


dU, 

ds 


ds dfi ' ; 


and as this must be equal and opposite to the force exercised 
by dd on rf.§, we have 

dl/, dU, 

ds ds' 


Again, as is a function of s and s', 


dU, dlh ,, 
dU\ ^ -—r- —rr 

ds ds 


{ds - ds'). 
ds 


dJJ 

Hence is a function of s - s', and therefore 
ds 

i/. =/,(«-. s')- 

111 like manner, 

IT, = f,(s - s'), U,=Ms-s'). 


222. Force on Current Flement In Magnetic 
Field,— If A, /i, V denote the direction-cosines of a current 
element ds, we have seen, Art. 219, that the components of 
the force which a magnet-pole exerts on ds are 

i 

{fiy - v]3) ids, (va - Ay) ids, and (A/3 — /la) idSj |h 

where a, 7 denote the components of the magnetic for^ 
due to the magnet-pole. ^ 
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225. Kleciric Ilis|ilaceiiieiit or I*4>larlzatt«»n. — 

When a conductor is electrically excited the conductors in 
the vicinity become electrically excited also, and a change 
is produced in the intervening medium or dielectric whereby 
at each point a directed or vector quantity is brought into 
existence in the medium. 

This directed quantity is called by Maxwell the electric 
dispheement^ and by Professor J. J. ^.I'homson the electric 
polarization. The latter term is no doubt scientifically the 
more correct ; but tlie word ' polarization' is used so frequently, 
especially in the theory of light, that MaxwelFs term is in 
practice the more convenient. 

In order to bring about this change in the dielectric the 
expenditure of work is required. If the electric displacement 
per unit of volume be denoted by i>, and its components 
hy /j .<5^, the expression for tlie total work 8Z7 per unit 
of volume, required to increase D by SD, is of the form 

X§/-i- F8r/ + ZU. 

The quantities by which 8/, S//, and are multiplied in 
tliis expression are called the components of the electromotive 
intensity 72. 

It will be shown that the vector quantity thus defined has 
proi)erties for the most part the same as those whicli belong 
to the electromotive intensity in the theory of action at a 
distance. 

Xhfdxdydz represents an element of work, X^fdxcly 
is of the nature of a mechanical force. Hence, if X be 
regarded as of the same nature as the force acting on the 
unit of electricity, fdxdy may be regarded as a quantity of 
electricity, and./ as a surface-density. 

In an isotropic dielectric wliose properties are the same 
in every direction, the electromotive intensity is co-dircctionai 
withy and proportional tOy the electric displacement. Hence foi 
sucli a dielectric we may write 

47r/-/iX, 4:ir(j^kYy 4wh ^ kZ. (1) 

The constant k depends on the nature of the dielectric: 
Mjd is called its specific inductive capacity. 

^ Since / is of the nature of an electric surface-density, 
)y (6), Art. 29, k must be a numerical quantity. 



Energy due to Electric Displacement, 


125 


226. Unergy due to Kleetric Dlis^placeiiieiit. — If U 

denote the energy per unit of volume due to an electric 
displacement, by Art. 224, we liave 

8?7=-XS/‘+ rSfj-hZSlL 

Substituting for A, F, Z, from (1) we get, by integration, 

ir. 

= ^(Xf+ Yg-^Z/i). (2) 

Hence the total energy IF, stored up in au isotropic 
dielectric © in consequence of an electric displacement, is 
given by the equations 

JF := I R-d® RDd®, (3) 

The second of the expressions for W given by (3) differs 
from that in Art. 77 only by containing the factor k, 

227. Conductors and Currents. — A permanent electric 
displacement cannot be set up in a conductor, but passes away 
immediately if not renewed. A displacement which is con- 
tinually passing on and being continually renewed constitutes 
an electric current. The intensity of a current is the rate of 
change of the corref<ponding displacement. When a conductor 
in electric equilibrium is situated in a dielectric in which 
there is a displacement, it constitutes a boundary to the 
dielectric ; and the surface integral of the normal component 
of the displacement taken over the conductor constitutes 
what is called the charge on the conductor. 

228. iioleiioidal Distribution of Dispiacement. — 

If a closed curve be drawn in a dielectric, and through each 
of its points a line be drawn in the direction of the electric 
displacement, we have what is called a tube of induction, or, 
in the language of Professor J. J. Thomson, a Faraday tube. 
Such a tube terminates at each end on a conductor, and, 
whatever be the electric charge at one end, an equal and 
opposite charge is found at the other. In an isotropic medium 
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tubes of induction are in the same direction as tubes of force, 
and are therefore at right angles to the surface of a conductor 
in electric equilibrium. Hence, if the tube be small, the 
positive displacement over the normal section directed into 
the tube at one end is equal in magnitude to the negative 
displacement directed into the tube over the normal section 
at the other end. Hence if Si and Sj denote the two normal 
sections, and i), and Zlj the two displacements in the positive 
direction of the line of induction, we have ZliSi = -DjSj. 

We conclude that, for any small tube of induction drawn 
in the dielectric, the product of the displacement and the 
normal section is constant. 

Prom this it follows that, if any closed surface 8 be 
drawn whose interior is occupied continuously by the dielec- 
tric, and if I, m, n deiicte the direction-cosines of the normal, 
we have 

f (If + my + «//) dS = 0. 

For, if ^ be the angle which a line of induction makes with 
the normal to the surface at any point. 


])% = D cos ■ipdS = (If t- mg + n/t) dS ; 

and, as every tube of induction is cut twice, or some other 
even number of times by the closed surface, 

J' (If + mg + n/i) d,S = cos ipdS = 0. (4) 


If the volume enclosed by S be the element dx dy r/s, we 
ebtain 


fdyd% 


-{f 


dy dz + y dz dx 


[ 

V ^(1/ ' 


dy ) dz dx 


f dh ^ 

h dx dy ^ r ^ ^ J ^ > 


that IB, 


(if dy dh 
dx dy (Iz 


(«) 


This equation expresses a fundamental property of the 
eleotrio displaoement, and is analogous to the condition 
fulfilled by the components of velocity in an incompressible 
fluid- 
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In tlie case of a conductor, /, h cannot exist except; in 
the form 


di 


df, &c. ; 


but the solenoidal condition is still fulfilled, so that for a 
conductor we have 


dt \dx dtj (h ) 


( 6 ) 


229 . Constancy of Chargee on Insulated Con- 
ductor. — If a conductor be insulated, its bounding surface, 
or surfaces, remains unchanged, and throughout the conductor 
hy (b) we have 

d df d d(j d dh . 

dx dt dy dt dz dt 


Multiplying by dxdydz^ and integrating, throughout the 
conductor we get 



df d(f dh 

dt dt dt 


dS^O; 


that is, 


d 

di 


If + mg + nhj dS = 


0 . 


Hence / (^ + nig + nh) dS^ taken over the surface or 
surfaces of the conductor, is constant. 

When a conductor is touched by anotlier conductor, the 
bounding surface of the space through which the integration 
is effected is altered, and there is no longer any ground for 
asserting the constancy of the charge. 

230 . Uisplacenieiit due to CSleetrlfied Upliere. — 

If a conducting sphere, of radius a, placed in an isotropic 
medium, be uniformly electrified, the lines of force and of 
induction are perpendicular to its surface and pass through 
its centre, since there is perfect symmetry round this point. 
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Hence the spuere is in electric equilibrium, and over any 
concentric sphere of radius r the displacement D is uniformly 
distributed ; and if Da denote the displacement at the surface 
of the sphere of radius a, we have 47 rr® 2 ) = iira^Do = e, where 
e denotes the total charge on the electrified sphere. Hence 


If a be suflScieutly small, we may regard the sphere as 
an electrified particle. 

The electromotive intensity R is given by the equation 


R 





( 7 ) 


and we have the result, that in an isotropic medium the 
force due to an electrified particle varies directly as the 
charge on the particle and inversely as the square of 
the distance. 


231. Energy due to two Klectrified 

lipheres. — Let the radii of the spheres be denoted by a and 
j3, and the spheres themselves by A and B. The electro- 
motive intensity due to the sphere A, on which there is a 
charge is by Art. 229, on the hypothesis that the oliarge is 


uniformly distributed, 7 “.;, where denotes the distance 
from the centre of the sphere. The electromotive intensity 


due to the sphere B is in like maimer 

It is plain that the resultant force may be derived from 
a potential function F, where 


AF= -+ - = + r.). 

Vx 

If W be the energy due to the spheres, we have, then, 
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The surface integral is to be taken over a sphere of infinite 
radius and over the spheres whose radii are a and (5. 

At the surface Si of the sphere A, if a be sufficiently 

small, , rr ^2 


where e denotes the distance between the centimes of A and J5, 


ei j dV2 

- + A* • 

fr f/v 


By Art. 26, 


dv 


and, as V is constant at Su we have 


Hence 


r^d8t = 0; 

dv 




a c ix 


In like manner, 


V"l 


dV ete.i 

— d Si = - 77 + — 

dv \p c , 


The integral over the sphere of infinite radius is zero^ 
also V’F'=0 througliout the field. Hence 


47r ej® 2^1 fj' 

T (,7 ^ ^ T" 


I fei\ ^ 2eiei\ 


If the sphere A were alone in the field, the expression 


above would become 


Similarly, if B were alone, it 


would be • Hence the energy due to the mutual action 
of the two spheres is ^ • 
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232. Force between Electrified Particles. — If W 

denote the energy due to the mutual action of two electrified 
particles, by Art. 231 we liave 


W = 


^ 1^2 
lev * 


where r denotes the distance between them. Hence, if be 
the mutual force which they exercise on one another, 


dW _ 1 ^1^2 
dr k 


( 9 ) 


Accordingly, the force between two electric particles acts 
in the line between them, and varies directly as the product 
of the quantities of electricity and inversely as the square of 
the distance. 

Also, by (7), Art. 230, the electromotive intensity due to 
an electric particle is equal to the force which it exercises on 
the unit of electricity. 

233. Irrotatlonal Distribution of Eleelromotive 
Intensity. — The components Z, F, Z of tlie electromotive 
intensity, due to a permanent statical distribution of elec- 
tricity, must be the differential coeflScients of an acyclic 
function of the coordinates. 

For, if we draw any closed circuit and suppose it occupied 
by a conducting wire, 



as as 




taken round the circuit, must be zero, as otherwise a perma- 
nent electric current would be set up in the wire without any 
expenditure of energy, which is impossible. 

Hence / (Xdx + Ydi/ + Zdz) between two points must 
be independent of the path, and therefore 


Xdx + Yd ff + Zdz = - dV, 
and V must be acyclic. 
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234. Distribution of DIectricity on Conductors. — 

If a conductor be in electric equilibrium, there can be no 
electromotive force acting in it, and therefore the potential 
is constant throughout. In the surrounding dielectric, 


4,/. -if, 4c., 


and, accordingly, from (5) we liave V = 0. The poten- 
tial V is therefore determined in the same manner as on 
tlie hypothesis of action at a distance. 

The charge on a conductor is J {If + hkj + nh) dS ; and, 
by (1), Art. 225, this is equal to 


A 

47r 


If-. 


where D and v are both drawn into the dielectric surrounding 
tlie conductor. 

Hence, if the total charge be given, so also is 

235. Conditions at Boundary between two Dielec- 
trics. — If two dielectrics, wliose specific inductive capacities 
are h and /ca, be in contact, at the boundary between them 
in passing from one to the other, V is continuous, as otherwise 
the electromotive intensity perpendicular to the boundary 
would be infinite. 

Again, the normal component of the displacement must 
be the same in one medium as in the other. To prove this, 
let us suppose two small tubes of induction resting on the 
same element of the boundary surface and drawn one in each 
medium. Let A and A denote the displacements, 2i and Sa 
orthogonal sections of the tubes drawn close to the boundary 
surface /S, and and ^2 the angles between the lines of 
displacement and the normal to 8. Then, by Art. 228, we 
have A Si = AS3 ; but 2i = dS oosipi, S2 = dS cos ^2, and 
therefore Di cos ^i = A cos \p 2 . 

The conditions stated above give the equations 



= Fa, 




dVr 

dvi 


+ 1(2 


dV. 

dv2 


0 . 


( 10 ) 
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If tlie positive direction be that of the normal drawn 
into tho medium whose inductive capacity is ku the second 
equation above may be written hNi = 1 hN 2 ; whence 

iv, - iV", = iV,. (11) 

k\ 


If we suppose k% to be greater than ki^ we see that a 
dielectric of greater inductive capacity, relatively to one of 
less, behaves like a conductor on which there is a charge of 


density 


Ih *“ ki 
^nki 


N 2 , In the case of a conductor, we must 


suppose kz infinite, then from (10) N^, = 0. In what precedes, 
Ni and denote the components of electromotive intensity 
normal to the boundary. 


236. Attraction on nieleetrle in Field of Force. — 

If a body composed of dielectric material be placed in a 
medium of different specific inductive capacity, the body in 
general behaves like a conductor in tending to move. 

To see the reason of this we must remember that, in 
general, if a conductor or a dielectric of different inductive 
capacity be introduced into a medium occupying a field of 
force, the total energy of the field is altered ; and, unless the 
field be uniform, the alteration is different according to tlie 
part of the field into which the conductor or dielectric is 
introduced. 

If a small change in the position of the conductor 
diminishes the total energy of the field, the conductor will 
have a tendency to move in the direction, producing a 
change of position whereby the total energy of the field is 
diminished. 

The same thing liolds good in the case of a body com- 
posed of dielectric material differing in inductive capacity 
from the medium by which the field of force is occupied. 


237. Cryiitalline Dielectric. — In a crystalline, or 
anisotropic, dielectric different directions differ in their 
electric properties, and the electromotive force is not 
necessarily co-directional with the displacement. 
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In this case, the components of the one are linear 
functions of those of the other, so that we have 


^Trf = 1(\\X. + h\iir + 

= A*2i-X^ + y " 

4i7rh = kz\jS^ + A*32 Y + k^Zm ^ 


( 12 ) 


If there be a function U of the components of force, 
representing per unit of volume the energy due to the 
displacement, we have 

S?7=X8/+ Yhj-^ZhU. (13) 


Substituting from (12) for S/, &c., in (13), and arranging, 
we get 

47r8C^ - (A’„X + A‘2 iF + kz,Z)lX + {knX + kzzY^l\.Z)hY 


but 

ilX. cljf (Ik> 


and therefore 

kzi — A’la, 

kz2 — /‘’23, A*31 ~ kl3 y 


and 





kiiX” -f- A ’22 

F* + hzZ^ + 2knXY 


also, 


+ ^kzzYZ + 2/»’i3 XZ ^ 

(14) 

II 

II 

11 

(15) 


By transformation of coordinates, 9>irU can be reduced 
to the form 

k,X^ ^kzY^^ hZ\ 

When XJ is reduced to this form, the coordinate axes are 
the principal axes of electric displacement, and Ai, /fa, denote 
the principal inductive capacities of the dielectric. 
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For the components of electric displacement we have, 
then, the equations 

‘^irf = kiXy Airg = IhYy ^wh = ( 16 ) 

If we take any point P of the dielectric as origin and draw 
the ellipsoid whose equation referred to the principal axes is 

kiiv^ + ka// + A-gS* = £®, 

it is plain that if we draw a line through P in the direction 
of the electromotive intensity, and draw a tangent plane to 
the ellipsoid at the point in which it is met by this line, the 
perpendicular on this tangent plane is in the direction of 
the electric displacement. 


238. Differential Dquation for Potential in Cry»» 
talline Medium. — If we express the principal components 
of displacement in terms of the electi’omotive intensity by (16), 
equation (5) becomes 


and therefore 


r/X 




0, 


/•i 




. . c/^F 


0 . 


(17) 


239. Distribution of DIectrieity on Conductors. — 

As there is no electromotive intensity in the substance of a 
conductor in electric equilibrium, F must be constant at the 
surface. 

In the dielecti’ic outside the surface, F must satisfy the 
equation 


. cPF . (PF , iPF ^ 


Also, the product of V and its differential coefficient 
integrated over a sphere of infinite radius must vanish. 

This appears from the consideration that 1/ + mg ■¥ nk 
integrated over a sphere of infinite radius is finite. Hence, 
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if E denote the radius of the sphere, / must be of the order 
but /, h are of the same order as the differential 

coefficients of V, Accordingly, V is of the order and 

dV 1 ^ 

V &C .5 are of the order 
dx ic** 

Finally, if the charge on each conductor be assigned, 


kj 


dj^ 

dx 


+ kim 


dV 

dy 


-f ki U 


(iV 
dz J 


]dS 


is given for each conductor. 

There is only one function Fwliich satisfies these conditions. 
If there were two, let ^ be the difference between them. 
Take the expression 



and integrate it by parts throughout the whole of space — 
the first term with respect to the second with repect to f/, 
and the third with respect to z ; then we get 



From wliat precedes, it is plain that eacli term on the right- 
hand side of this equation is zero ; and, as /c,, /t 2 , and are 
always positive, we have 

dx dy dz * 

and therefore ^ is constant, and consequently zero for the 
whole of space. 
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240. Energy expressed as liurface Integral. — If V 

denote the potential, and W the energy due to the electric 
displacement, by an integration similar to that employed in 
the last Article, by (14), Art. 236, we obtain 



Hence, by (16) and (17), we have 

2 IF = 2 J f^(^ + + nh) dS. (18) 

Since V is constant at the surface of each conductor, and 
since J {If 4- mg + nh) dS denotes the total charge on tlie 
conductor, equation (18) may be written 

(19) 

241. Energy due to Electrified Particle In Electric 
Field. — Let us suppose the field to be due to a single con- 
ductor, whose surface may be denoted by Si, on which there 
is a charge /»i. Let a conductor, whose surface may be denoted 
by S 2 , on which there is a charge ^ 2 * be introduced into the 
field. Let V denote the potential at any part of the field 
before the introduction of S 2 , and V + v the potential after- 
wards ; also, let fT and IF -i- w denote the total energy of 
the field before and after the introduction of S 2 . 'rhen we 
have 

2 JF== Viei, 2 (TF + tv) = (Fi + n) ei + (F 2 + Vz) 02 , 

If we now suppose 82 and <?2 to be infinitely small, so also 
is f?, and the term «? 2^2 is of the second order, and therefore 
negligible. Hence we have 

2w = F2^2 + ViVi. 
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On tlie hypothesis tliat 8 ^ and et are infinitely small, we 
have ViCz = ; for, if we integrate the expression 




, dV dv , dV dv , dV dv\ 

ki — -- + h T -T T" T 

ax ax ay ay az ilzj 


since V and v each satisfy equation (17), we get 


j’^( 


Iki ^ + mki '-y + nJcz '-y ) ( dSi + dS^ 
dx (ly dz 


) 




At the surface 8 \ tlie potential V is constant, and 



do 

dx 


+ ml'2 7 - + nh 
dy 



d8\ 


is zero, since the introduction of 8 ^ does not alter the total 
charge on 81 , Again, before the introduction of the con- 
ductor 8 i the total charge on the space surface 8 % was zero, 
and therefore 

f ( Ik, -y + mk. *‘y + nk, dS, = - 4ne,. 

J\ ^dx dy dzj 

Hence, as S 2 is infinitely small, 

\v{lk,^ + mh ^ + iih dS2 
j V dx dy dz J 

cannot differ from ~ 47rF2^2 by more than an infinitely small 
quantity of the second order. Accordingly the left-hand side 
of (20) is equal to - 47rF2^3. 

Again, as V o and V are eacli constant at 8 i, so also is t?, 
and therefore 
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Also, at the surface & we have 




iV , dV , dV\^^ _ 

— + tnki + ulh -T“ = 0 ; 

dx dy dz } 


and, since V v is constant, and V can vary only by an infi- 
nitely small quantity, tlie variation of v must be infinitely 
small, and 



7 

+ mlh ~r“ 
dy 




taken over the infinitely small surface /Sf2, cannot differ from 
zero by more than an infinitely small quantity of the second 
order. Accordingly the right-hand side of (20) is equal to 

TJTe have, then, V^e^ = viCi ; and therefore 2«r = 2 ^2^2. 
Hence, by bringing an electric particle e% to* a point where 
the potential of the field is F2, the energy produced is V^e%, 

It is obvious that the result arrived at above can be 
extended to an electric field due to any number of conduc- 
tors, so that in general Ve denotes the energy produced by 
bringing a small body having a charge of electricity e to a 
point where the potential is F. 

If, instead of supposing a small charged conductor intro- 
duced into the field, we suppose the charge on one of the 
conductors S\y already in the field, increased by the amount 
we can show in a manner similar to that employed above that 
8 TF, the increase of energy, is given by the equation 

In fact, 

2 fF = T^iCi + VzCi + Fs^a + &c., 
and 28 fF = Fi -f CiS Fi + ^28 F2 + &c. i 


but Fi 86-1 = eiSVi + C2S F2 + ^38 F3 + &c.. 


as may be shown in the following manner. 
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Let the original potential at any point of the field be 
denoted by F", and the increase of potential due to the intro- 
duction of Sci by V ; tlien, by an integration similar to that 
already employed, we have 

I ^ ^ £) 

r /, ,dv . dv . dv\f \ 

= j ^ ^ H rfoi + dSi + &c. j • 


At each of tlie conductors V is constant, and also V+Vy and 
therefore «?. 

Again, at each conductor, except the first, 




kd ~ + k^m ~ -f k^n ^^dS 
d,T dy dz ) 


is zero, and at the first this integral is 
Also, at each conductor, 


47r8(^i, 



dV , 

-p- + kirn 
dx 


(/V 

dy 


H k^n 



hence the equation above becomes 

47rFi8^’i = iir (ciVi + e-iVz + &c.), 

that is, 

= e,SVi + ezSV, -i + &c. (21) 

Hence we obtain 

SJr^ r,Se,y ( 22 ) 


and therefore we conclude that under any circumstances the 
energy produced by bringing a small quantity of electricity e 
to a point where the potential is V is denoted by Ve, 

242. i^ystem of €liarf(ed Conductora. — It is now 

easy to see that, if there be a system of charged conductors 
in a crystalline dielectric, the equations which hold good 
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between the charges and potentials are of the same form as 
those belonging to an isotropic medium, that is, 

Vi = pne, + piiCi ■\ pizCi + &c., 

Vi = J)\iC\ + P22^2 + pizCi + &C., 

V^ = Pv^Bl + /^23^2 + Jl?33^3 + &C., 

&c., &c. 

In fact, every step in the process by which these equations 
are proved in Art. 128 holds good here. 

For, from equation (16), it appears that if each component 
of displacement be altered in the same ratio, so also are the 
dijSerential coefficients of the potential and the total dis- 
placement. Accordingly, if the mode of distribution of the 
displacement be assigned, the potential at any one point 
varies as the displacement to which it is due. Also two 
systems of displacement which are each in equilibrium may 
be superposed without disturbing the equilibrium. 


248. Force on Flectric Farticle . in JBleetric 
Field. — We have seen that the energy due to an electric 
particle e in an electric field is Ve, 

If the particle receive a displacement whose components 
are Sir, S//, and Sz, the energy of the field is increased by 


fdV. dV. dV. 


This is the work done against the forces of the field 
which must therefore be 

dV dV , dV 

-’Hi' -‘Hi' -’Tz 

Accordingly, the force acting on an electric particle per unit 
of mass is the same as the electromotive force of the field. 

244. Potential due to Spherical Conductor. — If a 
charged spherical conductor whose radius is a be alone in the 
field, the potential V is constant at the surface of the sphere, 
and in the space outside satisfies equation (17). If we assume 

y/hn-y^ aud ^/k%Z,^Zy 
(17) becomes d^V (PV iPV 
^ rfi,’ dl^ 


( 23 ) 
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and, at the surface of the sphere, 

= al (24) 

We have therefore to find a function of tj, and ^ whiclk 
satisfies (2d), and which is constant when ly, ^ satisfy (24). 
This is the same problem as to find the potential of an ellip- 
soidal charged conductor. Hence the form of F is given by 
Ex. 3, Art. 76. 

If ki > ky > ksy we have 


rr = 


a'’ 


JO .. a® 

*- = r. 


h’ “ /n’ '' kx 

In-h 


Ji-i = (,a F=<?r — 

hh ’ kxh ’ J.y{(A*-A»)(X»-A»)l’ 

(2ft) 


where X is the greatest root of the equation 




= 1 . 


X» X*-A* X‘-k^ 

The constant C is determined from tlie equation 
dV 


(26) 


K 


„ dV - dV , dV 

Ikx -r- + mkt --r- + ««■ 
dx dy 


-"jdS^-ixre, 


where e denotes the charge on tlie surface 8 of the 
conductor. 

245. Force due to Spherical Conductor. — Diffe- 
rentiating (26) we obtain 

2Kd^ _ 01 A ( I! . 

X» - k^ ] (V - ky (A» - uy A‘ 

The quantity inside the bracket on tlie right-hand side 
may be denoted by and we get 

dX p^l 
dl X(X2-A*)' 

In like manner we have 


dX 


P^n 


d\ p^Z 
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Again, 


"'I = _i_ '!l 

dm dx 


d^ A ^ 
^.0,4o., 


accordingly, if X, y, Z denote tlie components of force, 
we have 

Y dV C /- f 1 

dx \‘-k^ y\(\^-h^){X^-k)\’ 

y^_dV_ C », 1 

df/ '"yh X x‘-/i^ y{ix-h%x‘-krr 

* X y{(X-^- /r)(A- 4-^)1 • 


246. Force due to diilierical Particle. — In the 

ease of a particle, a becomes infinitely small, and so also 
do Ic and h ; tlien 

A* = A' - = A- - Ir, 

and j)‘‘ = A* = 5“ + n® + 


Accordingly, we have 


r=£ _1_ v=9 L 

^“A*//a’ ■A»y/7.A' 

If we substitute for 5, »», and ? in terms of x, y, and s, 
we obtain 


where 


- V - y _ 

~kj?' ^~ky?' 


A* 


a? y* 
k^h 


+ 



(27) 


Hence we conclude that the force exercised by a spherical 
particle at a point P is not in the direction of the line joining P 
to the centre of the particle, and does not vary inversely as 
the square of the distance of P from the centre of the sphere. 
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From (16) and (27) we have 

Cr X _ Or y Cr z 

' ~ 4^X> r’ •' ~ 4^ r’ " ~ r’ 


(28) 


wliere r denotes the distance of P from the origin. 

Hence, the direction at any point of the displacement 
due to a spherical particle passes through the centre of the 
particle, but the magnitude of the displacement does not 
vary inversely as the square of the distance. 
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CHAPTER XII. 

ELECTROMAGNETIC THEORY OF LIGHT. 


247. lolroduetory. — The electromagnetic theory of light 
cannot be considered part of the theory of Attraction ; but it 
is so intimately connected with the properties of dielectrics, 
and with those of electricity and magnetism which have been 
explained in the foregoing chapters, that some account of 
Maxweirs great investigations does not seem out of place 
here. 

248. £nergy of Current in magnetic Field. — From 
the identity of the action of an electric current with that of 
a magnetic shell, in Art. 217 it was concluded that the 
potential energy JF oi a current in an independent magnetic 
field is given by the equation 

W ^ - t j (la 4 inb + nc) (18, ( 1 ) 

If i assume the infinitely small value Si, equation (1) 
becomes 

g IT = - 8/ J {la + mb + nc) dS. (2) 

This equation holds good whether the field be independent 
or not, as a change in the integral due to an infinitely small 
value of i must be infinitely small, and when multiplied by Si 
becomes evanescent. We cannot, however, regard the energy 
due to the presence of an electric current as potential, because 
the current is not a permanent natural agent whose action 
varies merely with its position. The current may cease, and, 
if so, its energy disappears. 

We must therefore consider the energy due to an electric 
current as kinetic. 
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In a dj'^namical system, if work be done against the 
natural motion of the system, the energy, if potential, is 
increased, but the energy, if kinetic, is diminished. 

Hence we conclude that, if the energy of an electrio 
current be kinetic, the expression for tlie variation of energy 
due to a variation of i must have the opposite sign from 
that which it would have if the energy were potential* 
Therefore, if T denote the kinetic energy due to the presence 
of an electrio current in a magnetic field, we have 

J {la + mb + nc) dS. (3) 

dT 

Since — S/, we obtain 

di 

dT 

— =/(/« + mb + nc) dS, (4) 


249. £fier|^y and Kleetroinoilve Force. — The con- 
nexion between variation of energy and force is given by 
Lagrange’s Equations, Dynamics, Art. 207. 

In the present case of the dynamical system consisting 
of electric currents in a magnetic field, the position of the 
system is specified by the geometrical coordinates of the 
various magnets and electrio circuits, and in the case of 
each current, by the distances along the circuit which the 
electrio molecules have travelled at any time since a definite 
epoch. If s denote tlie distance along the circuit which a 

ds 

molecule of electricity has travelled, its velocity is 

at 

Again, if p denote the density of the electricity, and a the 
area of a section of the circuit, the quantity of electricity 

ds 

which passes the section in the unit of time is but 

at 

this. Art. 212, is i the strength of the current. Hence, 

i i ds f ^ 

idt ~ pcf ds (s) - xp (so), 

] Iq jIq dt J 8q 

since p is constant, and a a function of «. Hence, if idt 

J ^0 
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and So be assigned so also is s. Accordingly, instead of 
specifying the position of a molecule of electricity by s we 

may do so by the coordinate tj, where r\ = idt. Again, 

J ^0 

since i is uniform throughout the circuit, r\ is the same for 
all the molecules of electricity. 

If now X\ Y\ Z' denote the components of the total 
electric force at any point of the circuit, Lagrange’s 
equation of motion, corresponding to the coordinate ly, is 

at (iTj (hi \ (iti (hi dril 

and if X, y, Z denote the components of electromotive 
intensity, the corresponding forces X", F", are given by 
the equations 

X" = /xrXfs, F" = p(T Yih, 2r" = po-Zr/s, 
also dri = 'i)dt = par/v. 

Hence 



From Art. 212 it appears that when a current is passing 
the electromotive force is opposed by the resistance of the 
circuit, so that the generalized component of force tending 
to increase rj is not E but E - RL 

Again, the kinetic energy T does not depend on ri but 
dT 

on ri or i. Hence — is zero always, and Lagrange’s 

equation of motion corresponding to the generalized co- 
ordinate 1 } becomes 


(it di 


E - Bu 


(5) 
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This equation may be written 


^ dt di ^ 


( 6 ) 


If T remain unchanged, (6) becomes Ohm’s equation (2), 
Art. 212. 

If T vary in consequence of a change in the electro- 
magnetic field, the electromotive force keeping up the 

current is diminished by — — . If this be negative, the 

electromotive force is increased. 

This property of currents is abundantly confirmed by 
experiment. It is indeed on this property that almost all 
the modern applications of electricity depend. It was 
originally discovered by observation; but its exact mathe- 
matical expression as given above is due to Maxwell. 

A simple dase of this phenomenon is exhibited if two 
currents which repel one another be made to approach. 
An additional electromotive force is then developed in each 
circuit tending to increase the current. 

This still holds good if E be originally zero in one circuit. 
A current is then produced tending to oppose the motion. 
Such currents are called ‘induction currents.’ It is on their 
existence that the whole theory of light as an electromagnetic 
phenomenon depends. 

The general principle exemplified in the production of 
induction currents may be expressed by the statement 

Li any circuit contained in an electromagnetic field every 
variation in the strength of the field produces an electromotive 
force which tends to diminish the variation, 

250. maxwell’s Theory of liight. — Maxwell supposes 
the entire universe to be filled with a dielectric called the 
luminiferous ether. 

If there be a variable electric displacement in any part of 
this dielectric, the variation of the displacement constitutes an 
electric current which produces an electromagnetic field. The 
variation of the current produces an electromotive force in all 
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the surrounding circuits. These electromotive forces produce 
currents which again give rise to otlier electromotive forces 
and currents, and so the original variable displacement is 
propagated through space. 

In the case of light, tlie original displacement is vibratory ; 
tliat is, it begins in a certain direction, increases in that 
direction up to a certain amount, and afterwards takes place 
in tlie opposite direction till it reaches the same amount as 
before, only in the opposite direction, when it ^is again 
reversed ; and this process is repeated so long as the light 
remains steady. 

The displacement is therefore quantitatively the same as 
the distance moved through by a vibrating particle, and may 
be represented by an expression of the form 

. 27r 

a sin — /. 
r 

The whole plienomenon may therefore be termed an 
electric vibration ; and, when propagated tlirough space, 
may be called an electric wave. 

From the results already arrived at, the laws which 
govern this propagation may be deduced, as will be shown 
in the following Articles. 

In the study of an electric vibration we have to do with 
five vector quantities : the displacement, the electromotive 
intensity, the current intensity, the magnetic force, and 
the magnetic induction. 

Let 

/, g, h denote the components of electric displacement ; 

Xy Y, Z those of electromotive intensity ; 

w those of current intensity ; 
a, fi, 7 those of magnetic force ; 
a, by c those of magnetic induction. 

We seek to determine difierential equations for the com- 
ponents of one of the vectors which will enable us to arrive 
at the laws of its propagation. 
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201. Magiietlc Inductloa and Electromotive In- 
tensity. — We have seen, Art. 249, that for any oirouit s, 
if X, F, Z denote the components of electromotive intensity 
due to current induction^ 


d dT [f^-dx ,j.du r^dz\ , 


m 


If we imagine a surface-sheet S filling up tlie circuit s, 
by Stokes’s theorem, Art. 192, the right-hand side of (7) is 
equal to 



and, by (4), 


Hence, 


d^ 

Ti 


{la + mb + nc) dS, 


^ (fo + mb + nc) dS • 



dX 

fdX 

__ 

dz 

\d!, 



dS. 


In the case of an electric disturbance in a continuous 
medium, this equation holds good for every circuit which 
can be drawn ; and therefore we have 

da ^ dY dZ dh dZ dX dc dX dY 

dt dz dy ’ dt dx az ’ dt dy dx 


252. Current Intensity and magnetic Force. — If 

we suppose a surface-sheet 8 drawn in the dielectric, the 
total current passing across it is denoted by 

J {lu + mv + mo) d8. 

The line integral of the magnetic force, taken round a 
circuit s,. bounding the surface /S', is due altogether to the 
current passing across 8^ since for magnetic forces due to 
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currents not embraced by s this line integral is zero. Hence, 
by Art. 215, we have 

4jr /(/« + mv + nw)dS = J + /3 + 7 ds 


me) dS = J j 


r% “ y * 


dy d(i 
dy dz 


[ da dy 


and, since this equation holds good for every circuit and 
corresponding surface whicli can be drawn in the medium, 
we have 

(9) 

dy dz dz dx dx dy 

253. Relation between Magnetic Force and 
Induction. — We have seen, Art. 201, that in a body 
magnetically isotropic, in wliicli tliere is no permanent 
magnetism, the components of magnetic induction are in a 
constant ratio to those of magnetic force, so' that 

a * wa, h = = wy. ( 10 ) 

In what follows, we shall always suppose the medium to 
be magnetically isotropic. 

254. Rquatlona of the Flectromagnetic Field and 
of Propagation of Rlsturbance. — In tlie general case of 
a dielectric electrically crystalline, collecting the results given 
by (16), Art. 237, by Art. 227, and by (9), (10), and ( 8 ) of the 
present Chapter, we have the following group of equations 
holding good in the electromagnetic field : — 



inf = 

KrX, 

ing -- 


4tTrh = KzZ. 

( 11 ) 



df 

di 

dg 
' dt 

= «?, 

dh 

dt 

= w\ 


( 12 ) 

dy 


inv - 

da 

dy 


da 

(13) 

■ dy 

dz ’ 

dz 

dx* 


dx 

dy 



a = nstt. 



c = 

■ zsy. 


(14) 

dY 

dY 

dh 

dZ 

dX 

dc 

dX 

dY 

/I f\\ 

dz 

dy' 

dt “ 

dx 

dz' 

dt 

dy 

dx 

(ID) 



Solution of Equation of Pvopagaiion. 151 

By (liffereutiatioii from (12) and by (13), &c., we have 


rlt^ 


du d f dy d(5 


^ ^ dt dt \dy dz ) w \ dg dt dz dt 


Hi. 


d dc d db 


w ( rfy \ df/ diV ) dz \dir dz J) * 


that is, 

^ d^f d fdX dV dZ 

47r® -i - ( -T + — + — 


dl} 


dx \ dx dy dz 

4,il 

l/f, ■ dx\K, dx Kt dy K, ds js 


If we assume 


we get 


/i*? - ^ ri‘i L- 


In the case of an isotropic mediiim, 
= 2?“ = C« = V\ 


(16) 


and we liave 


F^'vy. 


(17) 


255. Solution of CSquation of Propagation. — 

Equations (16) and (17) are very general in their character; 
and to obtain a solution suitable for the present investigation 
we must consider some of the characteristics of a ray of light. 

When light emanating from a point passes through a 
lens whose focus is at the luminous point, a cylindrical 
beam is obtained whose parallel sections are planes having 
similar characteristics in reference to the beam. We may 
assume therefore that one of the vibrations which constitutes 
the light is propagated so that its direction remains parallel 
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to a line fixed in gpaoe, and that at all points of a section of 
tlie beam parallel to a certain definite direction the vibrations 
are in parallel directions, and in a similar state or phase. 
Consequently, if D denote one of the displacements whose 
vibrations constitute the light, the direction of D is constant, 
and the direction in which D is propagated through space is 
also constant. 

If X, fXy V denote the direction-cosines of D, we may there- 
fore assume that X, /x, v are constant for all positions of 2>, 
and we have 

/ = X2), g = fxDy h = vD. 


Equation (17) assumes its simplest form when /is a function 
of one coordinate; and, as a particular case of (17), we may 
write 


dt* dz^ 


(18) 


By Art. 63, the solution of (18) is 
/ = (t>(rt- z). 

This expression for / indicates a variable quantity whose 
magnitude at a given point is continually altering and 
wliose every state or phase advances through space in the 
direction of z with a velocity F*. 

This is obvious, because 


0 { F (^ + ^') “■ (2 + «')} = 0 ( ~ 2)> provided Vt' *= / ; 


and, accordingly, the value of / at the point z at the time 
t is the same as the value of f at the point z + z' at the 
time t + t\ 

If 0 be a periodic function, the disturbance in the medium 
is called a wave. 

The distance between two points on the line of propaga- 
gation at which the disturbance is in the same state is called 
the wave-length. 

If T denote the period of the disturbance, that is, the 
length of time in which the disturbance at a fixed point P 
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passes through all its phases and returns to its original state, 
the wave-length is equal to Vr. For, during the period t, 
the original disturbance reaches a point Q whose distance 
from P is Fr, and the disturbance at P has during the same 
time returned to its original state. Hence, at the end of the 
period r the disturbance at Q is in the same state or phase as 
that in which it is at P, and consequently PQ is a wave-length. 

When a wave is passing through a niediun), the locus of 
the points at which the disturbance is in tlie same phase is 
called the wave-front. 

If the wave-front be a plane parallel to a plane fixed in 
space, the wave is called a plane wave. 

In the case of a plane wave, the direction of propagation 
is the normal to the wave-front, and the direction of vibration 
is parallel to a line fixed in space. 

We can now generalize the solution of (18) so as to satisfy 
(17), and to represent the propagation of a plane wave of 
electric displac^^meiit through the dielectric. 

We may assume ^ 

/=XZ>, g-^D, h = vDy D - ^[Vt-{lx'\-my -v uz)]y (19) 


where /, m, n denote the direction -cosines of a line fixed in 
space. 

Then V*P = (/* + and 




FV'; 


and, accordingly, (17) and the corresponding equations for 
g and h are satisfied, also D represents the displacement in a 
plane wave whose line of propagation is in the direction /, w, n, 

256. Direction of Displacement in Isotropic 

medium. — The expressions for — , &c., given by (12) 

at 

and (13), show that 
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If there be an electric displacement in the medium before 
the disturbance takes place, by (5), Art. 228, 


dx dy 



( 20 ) 


Hence this equation always holds good ; but 
f =^X<^[Vt - {lx + my + ms)) , 

with corresponding equations for g and A, and therefore 
by (20), 

{\l + ym + = 0 , 

and accordingly 

XI + /i/M + vu = 0, 

and we learn that in a plane wave the disturbance is per- 
pendicular to the wave-normal, and is therefore in the 
wave-front. 

This is often expressed by saying that the disturbance 
is in the plane of the wave. 

257. JHagnetlc Foree In lsotro|ile medium. — By 

(14), &o., we have 

dt Z3 dt 7ff\dz dy J 7^K\dz dy) 

= 47r {mv - Wju) 

Integrating with respect to we obtain 

a = in V(niv - ny)!) + constant. 


As we are concerned only with the magnetic force due to 
the disturbance, the constant may be neglected, and we have 


a = 47r VJD (mv - Mfe), ^ 
j3= iirVDinX - M, . 
7 = iTrVD{lfi - mX). 


( 21 ) 
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Hence the magnetic force is in the plane of the wave and 
perpendicular to the displacement, and its magnitude H ia 
given by the equation 

H = ^irVD. (22) 

258. Crystalline medium. — The solution found, Art. 
265, for (17) holds good for (16) with some modifications. 

In fact, if we assume equations (19) and substitute in (16) 
and the two correspoiiuing equations, we get 

r*A = A^X - 1{AHX mnji + Chiv), 

FV = -BV - m {AHX + 

V^v = C^v - n (AyX + C^riv). 

In the solution of (17) F is given, and we find tha 
are indeterminate, provided they fulfil the condition lugate 

IX + nin + v« = 0. ipsoid 

In the present case, when /, n are given, equations ^ 
determine F* and A, fi, v. • If we eliminate A, fi, v from (23 
we get a cubic equation to determine F*. The absol\jtft 
term of this equation is f. 

^*(/*-l) BHm CHn 

AHm jB*(m*-'l) Chnn 

A^ln B^mn (7®(w*-l) 



If we call this determinant Q, we have 
1 

— -1 — m n 

V 

Q-^A^B^CHmn i ^ z}. „ 

m 

, - 1 

I m 

n 




p-l 




0 m* w* 

p 

nP-1 n“ 

= A*B‘C^ 

0 m®~l 

p 


n*-l 


0 w* 


= 0 . 
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Hence one value of F* is zero. The oorrespondingp values 
of A, /u, V are proportional to 

I m n 
'^25 

but they have no physical import, as the displacement to 
which they belong is not propagated through tho dielectric. 

For each of the values of F* w’hicli are not zero there is 
^ corresponding set of vabj^^a ol indicating two possible 

witL?^^^*^® displace^eiit a given wave-front, 
by equations (23) by /, the second 

’ .a the third by w, and add, we get 

and accc ^ ^^2/x + E^myi + C^vn) (1 - - n^) = 0. 


and V /\ + my. + nv =- Oy (24) 

pendi 

we infer that the two directions of displacement corre- 
sponding to a given plane wave-front lie in the plane of the 
.^yave. 

If we multiply the first of equations (23) by X, the second 
by y, the third by i/, and add, we get 

F*(A’ + ia*+ V*) 

= + By + - (AVX + B^my + C^nv)(lX + my + nv)y 

and therefore, by (24), we have 

F* = ^*A’ + By + CV. (25) 

If Ai, yiy vi ; X 2 , ft 2 , vz denote the direction-cosines of 
the two displacements perpendicular to /, m, w, and Fi 
and Vz the corresponding velocities of propagation, we have 

{A^ - Fi*) Ai = 1{AHX, + B^my, + Chivi)y 

with two corresponding equations. 

Multiplying the first by X 2 , the second by yzy the third 
by vz and adding, since /A 2 + myz + nvz = 0, we get 

Fi* (A 1 A 2 + y]y2 >^ 1 ^ 2 ) = ^^AiA’i + B^yiyz + C^vivz* 
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111 like manner, we have 

F2*(\iA2 + JULlfUli + VWi) = -4*XiX2 + + C^ViV^- 


Consequentl 3 % unless Vi = F*, we obtain 


A 1 X 2 + + V 1 V 2 ^ Oy 

Hence we ^ 
belonging to t, 
other in the pla^s. 
ellipsoid whose equation is 

^ V + jBV + CV = constant 


(26) 

I . I - 


ment 
t each 

auu are also conjugate in the 



Since these two directions are perpendicular and conjugate 
to each otlier, they are axes of the section of this ellipsoid 
made by the wave- plane. 


259. Wave - Surface. — If a vibratory disturbance 
emanate from a point 0 and spread in all directions 
through a medium surrounding 0, the locus of points at 
which at any time the disturbance is in the same state or 
phase is called the wave-mrface. 

If the medium surrounding 0 be isotropic, the disturbance 
is propagated with equal velocities in all directions, and the 
wave- surface is a sphere having 0 as centre. 

If the medium be not isotropic, we may suppose a number 
of small plane waves to start simultaneously from 0 in all 
possible directions. Each of these is propagated with a 
velocity corresponding to the direction of its normal. The 
envelope at any time of all these plane wave-fronts is tlie 
wave-surface corresponding to the medium. 


260. Construction for Wave-liurface of Crystalline 
Medium. — When an electric disturbance takes place in a 
crystalline medium, the equations of Art. 258 enable us to 
give a construction by which the wave-surface may be 
obtained. 
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If we take any period of time and assume 
a = Atly b = Btly c = Ctly 
the ellipsoid, whose equation is 

tf s* 

I^i.ivyug. J 

For each of the values of F* which are n©lft8ticity, and 
corresponding set of values of X, ju, v, indicati* 

^jrections of displacement with a given wave-f ® centre 0 of 
auu ivi 2 . ''-.^4- «nnn.fioTio/'98^ h\Q eleotric dis- 
^'^'aent in a plane wave due to the disturbance. Draw a 
nt plane to Fresners ellipsoid perpendicular to OP; 
and d itg point of contact, and draw OY perpendicular to 
plane POQ. 

Then OT and OQ are conjugate; and, being also at 
®’*ght angles to eacli other, are the axes of the section of 
I^VesnePs ellipsoid. 

Let the direction-cosines of OP, OF, and OQ be denoted 
by Ml, vi ; A 2 , M*, ; A', v ; tlien Ai, jui, m are proportional 

to 

A / / ^ 

M 1 V 

¥• 


and therefore, since X'Aa + fxfULz + v'v^ = 0, we liave 
a^AiA2 + b^fi\ii2 + c\iVi = 0 ; 
that is, -^’*AiA 2 + + C^viv^ ~ 

Also, A 1 A 2 + Viv^ = 0. 

Hence, by Art. 258, OF must be the second possible 
direction of displacement in the wave-plane corresponding 
to OP, and this wave-plane must be POF. 

Draw 08 in the plane POQ perpendicular and equal to 
OP; then OS is the wave-normal, and its length is the 
distance through which the wave-front has advanced in the 
time ty. If OP be drawn in the plane QOP perpendicular 
and equal to OQ, the locus of Pfor all possible positions of OP 
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is a surface wliich touches at T the wave-front perpendicular 
to 08. 



To prove this, take on, the tangent to Fresners ellipsoid, 
QP, a point Q' infinitely near Q, and in the plane QOF 
draw OT' perpendicular to OQ' ; then 02" = 0Q\ and if a 
plane be drawn perpendicular to 0T\ it passes through 0Q\ 
and one axis of the section of Fresners ellipsoid made by 
this plane is infinitely near OQ' and, being an axis, is there- 
fore equal to OQ' and consequently to 02". Accordingly T' 
must be a point on the locus surface, and TT' a tangent to 
this surface. 

Again, draw TT' parallel to OF, and take on it 2"' 
infinitely near T. Then, since T2’" is perpendicular to 02’, 
we have 02’" equal to 02’ Again, since OQ is perpen- 
dicular to tlie plane TOT'\ the plane perpendicular to OT" 
passes through OQ, and the axis of the section of Fresners 
ellipsoid made by this plane, being infinitely near to OQ, is 
equal to it, and therefore to OT and OT". Hence T" is a 
point on the locus surface, and TT' a tangent to this surface. 

Accordingly, the plane STT" is a tangent-plane to the 
locus surface ; but this plane is the position of the wave-front 
at the time Hence the locus-surface is the envelope of all 
possible wave-fronts at the time ^i, and is therefore the wave- 
surface. 
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261. dquatioii of Waire-flurface. — It is now easy to 
find, in the manner of MaoGullagh, the equation of the wave- 
surface. 

If r denote the length of any radius- vector of Fresners 
ellipsoid, a sphere, having 0 as centre and r as radius, meets 
the ellipsoid in the cone whose equation is 



A tangent plane to this cone meets the ellipsoid in a 
section in which two consecutive radii vectores are equal to r. 
Hence the line of contact is an axis of this section, and 
the extremity of an intercept equal to r on the perpendicular 
to the tangent-plane to the cone is a point on the wave- 
surface. If r be regarded as constant, the equation of the 
cone reciprocal to the cone of intersection of the sphere and 


ellipsoid is 


■ Qfi* 4 /fj* 4 - 

d* - r* - r* ^ c’ - 



(28) 


The coordinates of a point on the wave-surface whose 
distance from 0 is r satisfy this equation. Hence, if 

r* = + y^ + 

equation (28) becomes tlie equation of the wave-surface. 

Eejecting tlie factor and getting rid of fractions, we 
have 

{fr" - 7’*) {& - 7’*) + b^ (c?® -- r®) (r;® - 7’®) 

+ c®s® (a® - 7’®) (i® - r) = 0. 

Arranging in powers of r, and dividing by 7’®, we get, 
finally, 

(a* iu® + ^® ]f + c® 3®) 7’® - a® (6® + c®);*?® - ft® (c® + rt®) 7/® 

-c®(a»+ft®)3® + a®ft®c® = 0. (29) 

The surface whose equation we have obtained was 
discovered by Fresnel, and is known as FresneFs wave- 
surface. 
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262. magnetic Force. — From Art. 254, we have 

da _ 1 da 1 (tjY (j^\ _ ^ / I 'IlL - 

dl CT <(i w ill/ ) w \^3 dz ICz dyj 

= 47r ~ - <7’ = 4tir{C"-vzin - 


Integrating with respect to t, we get 
« = ^(wC'-vi - 
Hence we have 


a = (mC'Sm - nB-fix), 

)3 = ~ (nA^-Xx - . 

' .* 

.y = - mA%). 


(■•30) 


From (30), we see that 

la + H- ny =* 0 ; 

also, 

Ai« + /ui/3 + vij 

- - riwO + ^ 

= (-4®AiX2 + C^ViVt) - 0. 

Accordingly, the magnetic force is in the wave-plane, 
and perpendicular to the displacement; that is, its direction 
coincides with the second possible direction of displacement 
in the wave-plane. 


M 
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Hence if H denote the magnetic force, we have 

ivitli two other corresponding equations. Multiplying first 
by As, second by /ij, third by vs, and adding, we get 

47r2> 


n = {A^X\(fi2n - vitn) + B^iJii(vil~Xtn) + C’viQ^im- 

A.» + + CW) = ^ F* = 47r FD. (31) 


V 


263. Electromotive Intensity. — If F denote the 
resultant electromotive intensity, and 9i, 02 , and 03 its 
:lirection-angles, we have 

F cos 0, = Z = ^ = AttzsDA^Xi, 

Ki 

with two similar equations ; then F is in the direction of OQ, 
Bg. 1 ; and if x denote the angle between the displacement 
md the resultant electromotive intensity, we have 

F cos X = F (Xi cos 01 + fXi cos 6% + vi cos 03 ) 

. = 47r^J9(^»Ai^ 4- 

iiid 

F = 4:ttz5 V^D sec x- (32) 

264. Conditions at a Boundary. — When a disturbance 
passes from one medium into another, six conditions must be 
(ul filled at the boundary ; but of these six, only four are 
independent. 

By Art. 228, the normal component of electric displace- 
ment must be continuous. Hence, if /, m, n denote the 
direction-cosines of the normal to the boundary, /, g, h the 
components of displacement on one side of tlie boundary- 
surface, and /', g\ // those on the other, we have 

I if “/') + ^ (/7 - g') + « (h - //) « 0. (33) 
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Again, the tangential components of electromotive 
intensity are continuous. In fact, each component of 
electromotive intensity must be continuous in a direction 
perpendicular to its own, as otherwise, by (15), tliere 
would be an infinite rate of change in the magnetic 
induction. 

Accordingly, if X, Y, Z, and X', F', Z' denote the 
components of electromotive intensity at the two sides of 
the boundary-surface, and Xi, jui, vi ; Xs, lun, vg the direction- 
cosines of two mutually perpendicular tangents to the 
surface, we have 

Xi(X-X') -f (F- F') + Z..(2r-F') = 0, \ 

X2 (X - XO + /ig ( F- F') + 1^2 if') = 0. j 

As the magnetic induction fulfils tlie solenoidal condition, 
eacli of its components must be continuous in a direction 
coinciding witfi its own, and therefore the component normal 
to the surface must be continuous. Hence we have 

I {a- a') + m (b - V) + n{c - c) - 0. (35) 

Also, by (13), the components of magnetic force tangential 
to the surface are continuous, and therefore 

Xi (a - a') + (/3 - /3') + VI (7 - 7') = 0, 1 

Xs (a - a) + jU3 (^ - j3') + V 2 (7 - 7') = 0, j 


Equation (35) follows from equations (34), as may be 
shown in the following manner; — 

By (15) we have 
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From equations (34) it appears that X-X', Y- Y\ 
and Z-Z' are proportional to (juiy2-- V1/X2), &c., that is, to 
/, w, n ; or, if U = 0 be the equation of the boundary-surface 

to 

dU dU ^ dU 
1 > j i and , • 
dx dy dz 

Hence, if A denote an undetermined function of the 
coordinates, and Q be put for 


1 



we obtain 

^ dU d ^ dU\ (d ^ dU d ^ dU 

dy dy dz)^ \dx dz dz dx 

fd dir d dU\ 


'dir dA 
jiy dz 


<[AdU\ f(m(lA ^(a(W\ 
dy dz / ^ \ dx dz dx ) 

(djr(lA^dAdJI\ 
^ ^ \ dx dy dx dy ) 


\ dx \ dy dz dy dz J dy \ dz dx dz dx / 
dz \dx dy dx dy J) 


In a similar manner, from equations (12), (13), and (36) 
we get 
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By integration, we obtain 

I {/ - /') + m (g - g') + n (h - //) = constant, 
l(a - a') + m {b - i') + w (^ - c') = constant ; 

but as we are here considering only the results of the 
disturbance, we must suppose /, g, h ; /', g\ // ; a, &o., to 
be all initially zero, and therefore we get (33) and (35). 

265. Propag^atioii of Ught. — If we suppose each point 
of a plane area S to be a centre of disturbance, and draw the 
wave-surfaces of which these points are the centres, and which 
all correspond to the same period of time a plane S', parallel 
to S, whicli touches one of these surfaces will toiicli them all ; 
and if we draw straight lines from the boundary of S to the 
points of contact with S' of the surfaces whose centres are on 
this boundary, the area S' enclosed by this cylinder is made 
up of points at which the disturbances are all in tliis plane, 
parallel to one another, and in the same phase. Consequently, 
X is the waf e-plane at tlie time ti. Outside the cylinder the 
plane S' does not tougli any of the wave-surfaces, and the 
disturbances due to wave-surfaces corresponding to a period 
different from ti are not in the plane S', nor parallel to one 
another, so tliat instead of strengthening they interfere with 
cfich other. Thus the sensible effect is limited to the area 
within the cylinder passing through the boundary of S ; 
accordingly, the light is propagated in a straight line, and 
the direction of the cylindrical beam or ray is that of a line 
drawn from tlie centre of one of the wave-surfaces to its 
point of contact with S'. 

If the medium be isotropic, the wave-surfaces are spheres, 
and the ray of light is perpendicular to the wave-plane. 

If the medium be not isotropic, the ray is in geueral not 
perpendicular to the wave-plane. 

266. Reflexion and Refraction. — When a disturbance 
advancing through a medium reaches the boundary of another 
adjoining medium, the continuity of propagation is inter- 
rupted. The most geueral hypothesis we can make is, that 
disturbances, starting from the boundary, are set up in both 
media. A small portion of the boundary between the two 
media may be regarded as a plane area, and we may suppose 
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a cylindrical ray of disturbance to reach this area. The plane 
containing the wave-normal of the incident ray and the normal 
to the boundary is called the plane of incidence. 

All the plane sections of the cylindrical ray which are 
parallel to the plane of incidence have a common perpendicular 
lying in the tangent plane to the boundary. 

We shall suppose at first that each medium is isotropic. 



Let AB be the line in which the boundary-surface is met 
by that plane of incidence which contains the longest of the 
parallel chords of the cylindrical beam of light. Let lA 
and be the lines of intersection of this plane with the 
cylindrical boundary of the beam. 

Draw AP perpendicular to I A. When the disturbance 
reaches -4, wave-surfaces start from A in each medium ; and 
when the disturbance at P reaches P, the wave-surface start- 
ing from u4 is a sphere having A as centre, and a radius equal 
to PB, There are corresponding wave-surfaces having their 
centres at all the points of the beam which lie on the plane 
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boundary of the two media. If we draw through B a perpen- 
dicular to the plane of incidence, a plane, through this line, 
touching the sphere having A as centre, touches all the wave- 
surfaces, and is therefore the wave-front of the reflected beam. 
A perpendicular to this plane will be in tlie direction of the 
reflected ray. 

If BQ be a tangent to the section of the wave-surface 
starting from A, the reflected ray is in the direction of AQ. 

We see, then, from the equality of the triangles AQB 
and BP Ay that tlie incident and reflected rays make equal 
angles with the normal to the boundary-surface. We liave 
seen above that the reflected ray lies in the plane of incidence. 
We have thus the two laws of reflexion in an isotropic medium. 
The direction of the refracted ray is obtained by a method 
similar to that employed for the reflected. 

Describe, with A as centre, the wave-surface belonging to 
the second medium and corresponding to the period of time 
required by tfee incident ray to travel from P to P. If the 
velocity of propagation in the second medium is less than in 
the first, the sphere in the second medium will have a radius 
AQ; less than PB\ and if i and h be the angles which the 
incident and refracted rays make with the normal to the 
boundary-surface, we see that the refracted ray is in the 
plane of incidence, and that 

sin h AQ' _ Vx 


where V and V\ denote the velocities of propagation in the 
first and in the second medium. 


The ratio is called 


the index of refraction of the two 


media’ ; and if we denote it by jul, we have sin f = /x sin ii. 

When the second medium is crystalline, its wave-surface 
will have two sheets, and two tangent planes can be drawn 
passing through the perpendicular at B to the plane of 
incidence. The corresponding directions of displacement are 
obtained by means of Art. 260, and the lines from A to the 
points of contact of the tangent planes are the directions of 
the rays. In a crystalline medium there is thus double refrac- 
tion, and a single ray of light becomes, in general, two rays. 
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267. Common Ught and Polarized liigbt. — In an 

isotropic medium the direction of displacement may be any 
whatever perpendicular to the ray. In the case of common 
light, the direction of displacement is not fixed, but after a 
few hundred vibrations passes into another direction in the 
wave-plane. In the case of light, some billions of vibrations 
are completed during a second, so that in any appreciable 
length of time we may consider that there are as many 
vibrations in any one direction in the wave-plane as in any 
other. When light is polarized, all the vibrations belonging 
to a given ray are in the same direction. 

We have seen that when light passes into a crystalline 
medium it necessarily becomes polarized. In fact, wlieii the 
direction of the ray is given, the tangent-plane to the wave- 
surface at the point where it is met by tlie ray is the wave- 
front, and the line in which this plane is met by the plane 
containing the ray and the wave-normal is the direclioii of 
vibration. ♦ 

268, Intensity of liight. — Tlie ultimate measure of the 
intensity of light is its effect on the eye, but indirectly we 
can ascertain how it depends on the displacement producing 
the liglit and obtain its mathematical expression. 

It is ascertained experimentally that if light emanate 
from a constant source, the intensity of the illumination of 
a small plane area perpendicular to the direction of tlie light 
varies inversely as the square of the distance from the source. 

We conclude from this that the intensity of light varies 
as the energy of the disturbance per unit of volume. In fact, 
if a disturbance emanates from a source 0 in an isotropic 
medium and spreads equally in all directions, the mean total 
kinetic energy remains constant, and the disturbance atany time 
occupies the space between two spheres whose radii differ by 
a wave-length. Since the wave-length is very small, the space 
occupied by the disturbance is represented by 4irr^X, wliere X 
denotes the wave-length. Hence if thekineticenergy be denoted 

T 

by the energy per unit of volume is g ; and this varies 
inversely as the square of the distance from the source. 
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Tlie simplest form of expression for a periodic disturbance 
producing a plane wave wliose front is perpendicular to the 
axis of X is a cos where 


and V denotes the velocity of wave-propagation. The 
corresponding velocity v of vibration is 

27r V . 

r — a sin fji. 


Hence the mean value of is 



1 

2w 


sin^ip dtp ; 


that is, 



Accordingly, tfie demsity of the medium being constant, the 

• TT* 

kinetic energy per unit of volume varies as or as the 

T 

square of the amplitude if r be assigned. ' If we now suppose 
a small plane surface to be illuminated by two similar sources 
of light, the rays from which are approximately perpendicular 
to the surface, and whose distances from it are equal, the 
disturbance due to one of these sources may be represented 
in any direction perpendicular to the ray by a cos 0, and that 
due to the other by a cos + c). 

The total disturbance is, then, 2a cos cos ((p + Je). In 
a short period € passes from 0 to 27r, and the mean value of 
the square of the amplitude is 



which is equal to 2r4^ Hence, if we suppose that the intensity 
of light is measured by the energy per unit of volume due 
to the disturbance, we find that the illumination given by 
two similar sources of light is double that given by one. 
Thus the conclusion already arrived at is confirmed. 
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269. Knergy due to electromagnetic Disturb- 
ance. — We have seen (Art. 248) that if T be tlie kinetic 
energy of a system of currents in an electromagnetic field, 

(IT 

= I [la + mh + nc)(lS, 


Since I’ is a liomogeneous quadratic function of the strengths 
of the currents, 




, dir dG . 

also rt = — , &c. 

dy dz 


Substituting in (4), and applying Stokes’s theorem, we have 


.dT 

' dl " 




but 


. i/x 

i -- = (tn, 
am 


wliere a is the orthogonal section of the current, and n tlie 
component of its intensity, and adm = Hence we get 

(IT 

5/- - ^l[Fu+ ILe) r/©. 


where the integral is to be taken throughout the whole of 
space. Now, by (13), 


(Ij (Ift 


whence, substituting, we have 
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where the volume-integral is to be taken throughout tlie 
whole of space, and the surface-integral over both sides of 
every surface separating two media, and over a sphere wliose 
radius is infinite. 

Since a, / 3,7 are each at infinity of the order — , wliere 

JLX 

B is infinite, the surface-integral at infinity is zero. Again, 
my - ^?/3 is the magnetic force in the plane of yz perpendicular 
to the normal to the surface S. By (13), such a force, being 
tangential, is continuous in passing from one side of the 
surface to the other, and therefore the corresponding surface- 
integral, when taken over both sides of S, vanishes. 

Hence, if we substitute for 


, &c,. 


ils ^ cm 
dy dz 

their equivalents r/, c, we obtain 

SttT = J (cia 4- + Cy) d^ ^ J (o® + j3* + 7 *) 


Substituting for a* + j3^ + 7 * its value from (31), we 
get 

87rr= 167rV^F*J5V©; 

whence T = 27r (37) 


111 addition to the kinetic energy of the electric currents, 
the disturbance produces potential energy W due to the 
electric displacement. 

If we substitute and for 


_ 1 _ ^ 


and 


1 

'K\ 


by equations (14) and (16), Art. 237, we get 

r = 27r + Bhf + Odi^) f/0 

= 27r J tu + CW) mm - 27r J txj V’^BVm. 

(38) 

Hence, if E denote the total energy per unit of volume, 
we have 


E = 47rwV^D\ 


(39) 
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270. auaiititles to be^^ectromagnetlc Disturb- 

an%t Refraetfon. — When ligl) that if T be the kinetic 
into anotlier, there are four quai^i by 

means of the equations holding good at the boundary. These 
quantities difiEer according to the nature of the media. 

When light passes from one isotropic medium into 
another, the direction, intensity, and line of displacement 
of the incident ray being given, the directions of the reflected 
and refracted rays are known by Art. 266, and we have to 
determine their intensities and lines of displacement. 

When the first medium is isotropic and the second crys- 
talline, the directions and lines of displacement of tlie two 
refracted rays are determined by Arts. 266 and 260, and also 
the direction of the reflected ray. We have, then, to find the 
intensity and line of displacement of the reflected ray, and 
the intensities of the two refracted rays. 

Similarly, when light passes from a crystalline into an 
isotroiuc medium, we have to determine the intensities of the 
two reflected rays, and the intensity and line of displacement 
of the refracted. 

Lastly, when both media are crystalline, we have to 
determine the intensities of tlie two reflected, and of the 
two refracted rays. 

271. Reflexion and Refraction. Isotropic Jlledia. — 

Polarized light passes from one isotropic medium into 
another : determine the intensities and directions of electric 
displacement of the reflected and refracted rays. 

Let D, D\ and A denote the displacements belonging to 
the incident, reflected, and refracted rays ; then we may put 

-D = a cos A = a' cos 0', Di = ai cos 
and we may assume tliat at the surface separating tlie media 
0 = ^' = 01 . 

Again, if we put 

X = \/ 2Trzs Fa, 1' = ^ Zirw Fa', Ji = 2;rwi FiUi, 

wliere F and Fi cfeuote the velocities of propagation of the 
incident and refracted rays, we have i* = F*a* « mean 
value of Ey by Art. 269. 
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where the volume-integral i§ the intensities of the incident, 
whole of space, and the sm 

t\ectriop, xar is sensibly the same, so 
that we may assume wi = w. 

Let the normal to the separating surface drawn into the 
second medium be the axis of X, and the plane of incidence 
the plane of XY. Then, by Art. 266, the axis of Zis the 
line of intersection of the three wave-planes. 



Let 01 be a wave-normal or ray, and the plane OZDSM 
the corresponding wave-plane, OD the direction of displace- 
ment, and OMy perpendicular to OJ3, the corresponding 
direction of magnetic force. 

Equations (21) show that the magnetic force is perpen- 
dicular to the wave-normal and to the displacement, and is 
so directed that seen from it the wave-normal must be turned 
counter-clockwise in order to coincide with the displacement. 

Let the displacements make angles 0, O', and 0i with OZy 
and let the wave-normals make angles i, and ii with OX 
tljen, by Art. 266, wo have 


. . Fi . 

sin *1 = -p sio 


where V and Vi denote the velocities of wave-propagation in 
the first and in the second medium. 
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If X, F, X; X\ Y\Z'\ Xi, YiyZi denote the compo- 
nents of electromotive intensity, and a, /3, 7 ; /8', 7' ; 

«i, /3i, 71 those of magnetic force corresponding to the three 
xays, by Art. 264 wo have 

= i3 + i8' = i3., 7 + 7' = 7i- (40) 

By Arts. 257, 254, we have 

H^iTrVDy F = 47rz^F*2>; 
also, from fig. 3, we see that 

F= - F sin 0 cos Y' = F'sin 0'cos e, Fi = - Fi sin 0^ cos 

^ = F cos 0, X' = F' cos 0', = Fi cos 0^, ) 

(41) 

/3 = - jETcos 0 cos /3' = JT'cos 0' cos /, /3i = - JSfi cos 0i cos 

7 = - jBT sin 0, 7' = - ^T' sin 0', 71 = - i/i sin 0i. ) 

(42) 

As stated above, we may assume tui = w, and if in 
equations (40) the members of the first two be divided each 
by 2 ^ 27 rzjy and the members of the last two be multiplied 

each by » these equations by (41) and (42), when 
2-v/27r 

IS substituted for -rp, become 
sin ti y I 

sin i cos i {I sin 9-1' sin 0') = Ii sin ti cos h sin 0i, 

sin i(I cos 0 + /' cos 0') = Ii sin ii cos 0i, 

cos e (/ cos 0 cos 0') = Ii oob ii cos 0i, 

/ sin 0 + r sin 0' = ii sin Oi. 

From the first and last of equations (43) we obtain 
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And from the second and third we have 


21 cos d = 


sin (i + ii) 
sin i cos i ' 


^ sin (ii- i) T a 

2/ COS0 = — r-A . h cos. 01. 

sin i cos ^ ; 

If the displacement of the incident ray be perpendicular 
to tlie plane of incidence, 0 = 0, and hy (44) we have 
Ii sin 01 = 0, whence Ii = 0, or 0i = 0 ; but if we adopt the 
former alternative, by (45) we have / cos 0 = 0, which is 
impossible. Hence 0i = 0. In like manner, we get 0' = 0 ; 
and we learn that in tliis case all the displacements are 
perpendicular to the plane of incidence. Again, from (45), 
we have 


^ sin 2i j. 

r: -t. 

Sin + ?i) 


I' = sinQ'i-O j 
sin (^ + ii) 


If the displacement of the incident ray be in the plane of 
incidence, 0 = ? , and frdm (45) we see that 


as otherwise, by 44, we should have / = 0. Hence in this 
case all the displacements are in the plane of incidence. 
Also, by (44), we have 

j. ^ 2 sin 2i ^ sin 2i ~ sin 2/i ^ 

^ sin 2i + sin 2ix ’ sin 2^ + sin 2i^ 

III this case, if sin 2/i = sin 2/, we have /' = 0 ; that is, 
there is no reflected ray. When 2/i = tt - 2e, 




= tan i ; 


that is, if juL be the index of refraction, ^ = tan’^ /ul ; and we 
learn that if the tangent of the angle of incidence be equal 
to the index of refraction, there is no reflected ray when the 
displacement of the incident ray is in the plane of incidence. 

The displacements belonging to common light may be 
resolved each into two components, one in the plane of 
incidence, and the other perpendicular to that plane. 
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The whole of the reflected light is produced by the latter 
displacements when the angle of incidence is tan'^^u. Thi& 
light is therefore polarized, and, if it be made to impinge on 
a second reflecting surface so that the second plane of 
incidence is perpendicular to the first, there is no reflected 
ray when the tangent of the angle of incidence at the second 
reflecting surface is equal to the corresponding index of 
refraction. 

The discovery of polarized light was partly based on the 
observation of the phenomenon stated above, and common 
light, when reflected at the angle of incidence tan"^ ju, was 
said to be polarized in the plane of incidence. 

The plane of polarization as thus specified is perpendicular 
to the direction of the electric displacement which produces 
the light. 

272. Reflexion and Refraction. Crystalline 
Redluin. — Polarized light passes from an isotropic into a 
crystalline medium : find the intensity and direction of dis- 
placement of the reflected ray, and the intensities of the two 
refracted rays. 

Adopting a notation similar to that employed in the last 
Article, and putting i^, 2>', A? a^^d i >2 for the displacements 
belonging to the four rays, we may assume that at the 
boundary-surface ^ =r tp^. 

As before, take the plane of incidence for the plane of 
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Then all the wave-planes pass through the axis oi and 
all the wave-normals lie in the plane of XY, 

Let Z0D8 be the wave-plane, OI tlie wave-normal, and 
OB the direction of displacement corresponding to' one of tlie 
refracted rays. Then, by Art. 260, the plane 10 D contains 
OQy the line of direction of the electromotive intensity, and 
QOD is the angle denoted by x Arts. 260 and 263. 
Equations (34) and (36) become, in this case, 

/3 + /3' = + /Bs, 7 + y = 7i + 72 - j 

Also, by Art. 263, we have Fi = 4iiTTsV‘\ Di sec 

Now = Fi cos Qi y ; but (fig. 4) from the spherical 
triangle QJiY we have 

cos Qi y = cos Qi li cos IiY ^ sin Qi h sin h Y cos Qih Y, 
and QJ,=l + \u QJ,Y^-ir-I)OS='^+9,. 

Hence cos Q^Y - - sin %i sin h - cos X\ sin Oi, 

A similar equation holds good for Q.y Substituting in 
the value of Yi, given above, and using /, /i, and I 2 as in 

Art. 271, we get, instead of Yi^ the expression 

- I sin ii sec Xi(shi xi Xi cos i^ sin di) ; 

that is, - Ii (sin n cos Z, sin 9i + sin^/'i tan xi)* 

Again, Zi = Fi cos QiZ = Fi cos xi cos Oi. 

Hence, instead of Z\, we get Ii sin i^ cos Oi, 

The expressions to be substituted for tlie magnetic forces 
are similar to those made use of in the case of isotropic media, 
'i’hus equations (48) become 

sin Z cose (/sin 0 - /'sin 0') > 

= /i (sin ii cos ii sin 9i + sin* Zi tan xi) 

+ 12 (sin ii cos h sin Oz + sin* 4 tan xa) ? 

sin i (/ cos 0 + /' cos 0') = / sin i^ cos 0i •+ /g sin iz cos 02 , 
cos i(lcos0 - 1' COS0') « Ii cose'i cos0i + Iz cos iz cos 0^, 

/sin 0 + /' sin 0' = Ii sin 0| -f /* sin 0a. ^ 
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273. Uniradlal Directlonsi. — When the angle of 
incidence is given there are two directions of the displace- 
ment belonging to the incident ray for which there is only 
one refracted r.iy. 

To find one of these directions we may suppose I 2 zero in 
equations (49), and determine 0 in terms of i, A, and 0i. 

Making I 2 equal to zero, and eliminating T sin 0' from 
the first and last of equations (49), we get 

I sin 2i sin 0 = /i {sin (i + ii) cos (i - h) sin 0i + sin* ii tan xi ) ► 

(50) 

In like manner, from the second and third we obtain 
I sin 2i cos 0 = Ji sin (i + ii) cos 0i. 

Hence, by division, we get 

tan 0 = cos (i - ij) tan 0i + ~ t?fn vt. (51) 

^ ^ sill (e 4 «i) 

The second value of tan 0 is obtained by putting 02 and h 
for 01 and ii in (51). 

274. IJiiiaxal Crystals. — In the case of what are called 
uniaxal crystals, Fresnel’s ellipsoid is a surface of revolution. 
If we suppose c = 6 in the equation of the wave-surface (29), 
Art. 261, that equation becomes 


r* { 4 6* (/f 4 s*) j - a* ft* 4 2 /* 4 s*) - ft^a;* - ft* (y/® 4 «*) 

4rt*ft* = 0; 

that is, (r® - ft®) (y/®.i7® 4 ft® ( 2 /* 4 s®) - a® ft®) =0. (52) 

But (52) is the equation of the surface composed of the 
sphere whose equation is r* = 6*, and the ellipsoid of revo- 
lution whose semi-axis of revolution is ft, and whose other 
semi-axis is «. 

If a > ft, and Fresnel’s ellipsoid is prolate, the ellipsoid 
forming part of the wave-surface is oblate. 
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These conditions hold good in the case of a crystal of 
calcium carbonate, commonly called Iceland spar. This 
crystal is very celebrated in the history of science, ns 
observations of its behaviour led to the discovery of double 
refraction and of polarized light. 

The axis of revolution of Fresners ellipsoid is coincident 
with the line which is called the axis of the crystal. This 
line is the axis of symmetry, and can be determined from 
the geometrical form of the crystal. 

In the case of an uniaxal crystal, all rays inside tlie 
crystal wliose directions of electric displacement are perpen- 
dicular to the axis are propagated witli the same velocity. 
This appears from (25) by making C = B and X = 0 ; then 
= ly and Conversely, if 

cos* Sr + 5* sin* Sr = /i*, 

w'e liave (A^ - JS*)cos*^ = 0, and therefore 3* = For 

these rays the wave-surface is a sphere. 

Again, if the wave-surface be an ellipsoid of revolution, 
since the normal to a surface of revolution meets the axis, 
the ray, the wave-normal, and the axis must be in the same 
plane ; but the plane containing the ray and the wave-normal, 
by Art. 260, contains the direction of electric displacement. 
Hence this direction is in the plane containing the ray and 
the axis. 

Wlien light passing through an isotropic medium is 
refracted at the surface of an uniaxal crystal, one refracted 
ray is refracted in the same manner as if the crystal were 
isotropic, since the wave-surface of this ray is a sphere. 
This ray is called, therefore, the ordinary ray. The other 
refracted ray, whose wave-surface is an ellipsoid of revolution, 
is called the extraordinary ray. 

Eotli rays are polarized, and as a result of experiment it 
is said that the ordinary ray is polarized in the principal 
plane. By the principal plane is meant the plane passing 
through the refracted ray and the axis of the crystal. Hence 
we see again that the direction of electric dis[)lacement is 
perpendicular to the plane of polarization. 
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275. Uniaxal Crystal. Reflexion and Refrac- 
tion. — III the ease of an iiniaxal crystal, since 
equations (49) lieoome 

sin i cos i(I sin 0-1' sin 0') = Ii sin cos ii sin 0i \ 

+ I2 (sin /2 cos h sin02 + sin*/2 tan X2), 
sin /(/cos 0 + i'cos 0 ') = L sin/i cos 0 i + Jg Bin«2COs02A ( 53 ) 
cos /(/ COS 0 - T cos O') = Ii cosn cos 0 i + 7 ^ cos/g cos 62, 

7 sin 0 + 7' sin 0' = 7i sin 0i + I 2 sin 0^, , 

As an example of the use of these equations, we may 
sup])ose light to fall on the surface of an uniaxal crystal cut 
perpendicular to the axis. 

In this case, since the axis of the crystal is the normal to 
tlie surface, the jilane of incidence contains the axis and the 
wave-normal of the extraordinary ray, and, consequently, the 
ray itself. Hence both refracted rays are in the plane of 
incidence; and, by Art. 274, we have 

0. = 0, 0, = 

Accordingly, equations (53) become 

sin i cos i (7 sin 0 - T sin &) \ 

= 72(sin H cos H + siii *?’2 tail X 2 )t 
sin i(I CO80 + r QOS O'] = f sin /i, s. (^ 54 ) 

cos e*(7cos0 - T CO80') = Ii cos fi, 

7 sin 0 + 7' sin 0' = I 2 , ; 

If we now sujipose that the incident light is polarized in 
tlie plane of incidence, 0 = 0; and from the first and last of 
equations (54) we have 

I 2 (sin i cos i + sin h cos ^ + sin */2 tan Xi) = 0. 
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Since the expression by which is multiplied cannot be 
zero, we get = 0, and therefore sin O' = 0. The second 
and third of equations (54) become, then, 

sin i (/+ /') = /i sin A, cos i{I - I') = h cos ; 

whence we get 

2/ sin i cos^ - Ii sin (/i -h «), 2/' sin i cos i = Ii sin (?i - /). 


Finally, we obtain 


^ j. sin 2i 
^ sin (/ + ii) * 


( 66 ) 

sin {t 4- /i) 


Again, if the incident light be polarized in a plane 

TT 

perpendicular to the plane of incidence, 0 , and the 

second and third of equations (54) become 

r sin i cos O' = Ii siji n, - V cos i cos O' = Zi cos i\. 


Hence we obtain h sin {I + A) = 0, and therefore /j = 0 ; 

TT 

whence also cos O' = 0, and O' ^ • 

From the first and last of equations (54) we liave, then, 
sin i cos i [I - 1') = 1% (sin A cos u + sin* A tan ^^ 2 ), / f /' = Iz, 


Wlieiice 


/' = / 


/2 = /- 


. sin i cos i - (sin A cos A + sin* A tan 
sin i cos i + sin A cos A + sin* A tan ^2 
2 sin i cos i 

sin i cos i + sin A cos A + sin* A tan X‘> 


These expressions can be put into a simpler form. 

The angle x is the angle between the directions of electro- 
motive force and electric displacement, and is measured from 
the former towards the latter in the same direction as the 
line of displacement is turned in order to become the wave- 
normal. This appears from the figures and formulae of 
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Arts. 260 and 272. What has been said amounts to this — 
that in equations (53) ^ is to be regarded as positive when 
tlie direction of tlie electromotive force does not lie between 
those of the displacement and the wave-normal, and conse- 
quently the ray does occupy this position. In the present 
case the axis-minor of the wave-ellipse is the normal to the 
surface, and the positive angular direction is from it to the 
refracted wave-normal. The refracted ray lies farther from 
tlie axis than the normal, and consequently does not lie 
between the electric displacement and the wave-normal. 
Hence in (56) the angle x* is negative. 

An exi)ression for tan x «an be found by the geometry of 
the ellij)se. 



In the figure OX represents the axis of the crystal, OQ the 
line of electromotive force, OP that of electric displacement, 
01 the extraordinary wave-normal. Then x is the angle 
QOP; but QOP^EOl and 


tan ROI = 


El 

01 


UI. 01 2 triangle ROI 

OP OP 


Now, if Pi and ]h be the focal perpendiculars on the 
tangent, and ti and h the intercepts on the tangent between 

their feet and the point of contact, ~ and therefore 


(pi +i?2) {ti - ^ 2 ) = {pi - Pi) (ti 4* 12 ) ; 
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but i {pi 4- P 2 ) {ti - ifa) is double the area of the triangle 
-RO/, and \{pi - p 2 ){ti-v U) is double the area of the 
right-augled triangle whose sides are - 6*)sin/a and 

v/(a* - J*) cos /a, the angle XOI being h. Hence, if 01 be 
denoted by p, we have 


tan ROI = 


[a^ - h^) sin 4 cos h 
P' 


In the present case X2 is negative, and we have 

.... . ■ . {a^ - b^) sin* /a \ 

sin ^2 cos ?2 + sm* #2 tan Va = sin 1 % cos /a 1 - ~ r ' « V ~ TT 

V r»*sin*?2 + i*cos*«2/ 

i* sin /‘a cos h ii* sin u cos i% 

a* sin* h + J* cos* h A? sin* h + cos* h 


Again, if V denote the velocity of propagation in the external 
medium, 

F* sin* H = (.4* sin* h + cos* it) sin* i ; 


whence 


siir /2 = 


R* sill* 


8in*F 


cos* /a = 


F* - -4* sin*! 
F^“(A* - J?*) sin 27 ^ 


- 4 * sin* /a + oos*?2 = — j-ji — 

F* - (A*- jB“) snri 

Hence 

^* sin it cos /2 R sin i -v/ F* - A* sin^/ 
- 4 * sin*?i + B cos*/2 

and 

sin/ cose + sin/2 oos/a + sin^h tan xa 

sin / ( F* cos / + -B v^F* - - 4 * sin*/J 

sini cos/ - (sin/j cos/j + siii^/j tauxa) 

sin / ( F® cos / - ^ v/ - - 4 * sin®/| 

_ ^ 


(57) 


( 53 ) 
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Accordingly, 

F* cos i - B v/( F^ - siii^?:) \ 

I* (59) 

2 F* cos I 

^ F* cos i^B^{ F'^ - -4* sill* «•) ‘ ^ 


If tlie value of i be such that /'= 0, the reflected ray, 
when common light falls on the crystal, is polarized in the 
plane of incidence. This value of i is called the polarizing 
angle of the crystal when cut perpendicular to its axis. 

Making = 0 in (69), we have 


whence 


F'(l - sin*0 = B\ F* ~ sin^O ; 


sin^i = 


F* ~ AUP * 


(60) 


276. Reflexion and Refraction at Interior Surface 
of Crystal. — When light passes from the interior of an 
uniaxal crystal into an isotropic medium, there are, in general, 
two reflected rays ; and when the incident ray is an ordinary 
ray, wo have 

sin /i cos ii (Ii sin 6i - 1\ sin 0'i) 

+ (sin i 2 cos sin 9\ + sin’*^'^ tan I'l 

= /a sin fa cos '4 sin 6.i, 

sin /i(/i cos 01 + /'i cos 9\) + 1'% sin i\ cos 0'o, 

r a \ ( 61 ) 

= Jz sin /a cos dz, ^ ’ 

cos /i (/i cos 01 - 1\ cos 0'i) + /'a COS i\ cos 0'a, 

= Ta cos H COS 03, 

/i sin 01 4 1\ sin 9\ 4- Fz sin Q'z = h sin 03 . 
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When the inciden^l^ is an extraordinary ray, the 
equations at the r^ftacting surface become 


{SxLi H cos /2 sin 62 + sin’* u tan X2) I2 ] 

+ I\ sin i\ cos i'l sin 0 \ 

+ (sin i'2 cos i'2 sin 0'^ -f sin^e'a tan X2) ^'2 

= Is sin 4 cos 4 sin 0^^ | 

I2 sin 4 cos 03 + I\ sin i\ cos 0 \ + I's sin i'2 cos O'z ^ 

= Is sin 4 cos 03, 


(62) 


I2 cos 4 cos 02 + I\ cos i\ cos 04 4 74 cos i'2 cos O'z 

= /g C(^s 4 cos O3, 

/2 sin 02 -f /4 sin 04 + 74 sin 04 = 73 sin 03. j 


When the crystal is cut perpendicularly to its axis, 

0, = 04 ^ 0 , o, = 04 = ^. 

In this case, the first and last of equations ( 61 ) become 

(sin /4 cos i'2 + sin’*r2 tan ^ 4 ) 74 = Is sin 03 sin 4 cos 4 , 

7 2 = 73 si n 03. 

Hence 

Is sin 03 { sin 4 cos 4 - (sin i'z cos i'z 4 sin^ /4 tan ^ 4 ) } = 0 ; 

but tlie multiplier of /3sin03 in this equation is nor, in 
general, zero, and therefore we have 

sin03==(), 74 = 0. 


Consequently there is no extraordinary reflected ray, and tlie 
refracted ray is polarized in the plane of incidence. From the 
second and tliird of equations ( 61 ) we then obtain 

^ sin (4 - 4 ) r T 24 

J. \ J. \ . y • • \ > -^3 1 • "/ • » \ * 

sin(?3 + ei) sin(?3 + *,) 


( 63 ) 
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If the ray incident on the 
be the extraordinary ray, and the ^ 

perpendicularly to its axis, since all the rays ni. 
normals are in tlie plane of incidence which outs the wavc- 
ellipsoid in an ellipse whose axis-minor is the axis of the 
orystal and also the normal to the surface, we have 

i'2 = tt - 1 ,, x'2 = - X2 ; 

whence 

sin i'2 cos + sin^2'2 tan x^2 = “ ^2 cos 4 + sin*^ tan Xa)- 

In this case, the second and third of equations ( 62 ) become 
I' I sin i\ - Iz sin ^ cos ^3, cos i\ = cos h cos 63, 


whence Iz cos 0 ^ sin (u - i\) = 0 ; but sin (4 - i'l) cannot be 
^ero, and therefore cos 03 = 0 , and i'l = 0 . Consequently, 
there is no ordinary reflected ray, and the refracted ray is 
polarized in a plane perpendicular to* the plane of incidence. 

The first and last of equations ( 62 ) now become 


(sin h cos h + siii^h tan Xa) (4 - i'z) = Iz sin 4 cos /y, 

I2 + /'2 = / 3 . 


lienee 


jf - / h COS iz + sin^/g tan Xa - sin ^ cos 4 
* * sin iz cos /a + sin* 4 tan Xa + sin is cos is * 

2 (sin I2 cos ?2 + 8111*1*2 tan Xa) 

2 z = I2 ~ : :: . « • 4. i ^ ^ * 

sin t2 cos I2 + sin*i2 tan x* + sin 13 cos 13 


By reductions similar to those effected in the case of 
•equations ( 56 ) we get 


i'2 = /a 


B -fi^Vz - ^ sin^/s) - Vz cos iz 
/y v^( - - 4 * sin* 4 ) + Vz cos 4 ’ 

2 gv/(F 3 *-^*sin*i 3 ) 
B^/{Vz - - 4 * sin* 4 ) + Vz cos 4 


( 04 ) 
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^*^ystal, bounded by faces parallel to each 
the iiieifjia to the axis, be placed in an isotropic 
H JOBS at the bgbt polarized in a plane perpendicular 

’ tuc plane of incidence be transmitted through the crystal, 
tlie incident and emergent rays are parallel, and the plane of 
polarization remains unchanged. Then /s = t, and in virtue 
of equations (5U) equations (64) become 

2 F* { By/( F* - sin^e) - F^ cos / } cos i 

= “{ J57(p - sin*0 + F^s 

j 4F*i?v/( F* - -4* sin*e) cos e 

' Yn^( - A^ miH) + F^ ^j* ‘ 

When the incident ray falls on the first surface of the 
cry still at the polarizing angle, we have 

/ = 0, i '2 = 0, and /a =/, = /. 

In this case, the incident light passes through the crystal 
unchanged in intensity, direction of electric displacement, 
and direction of propagation, 

277. illiigularitles of the Wave-Surface. — The 

equation of the wave-surface. Art. 261, may be put in 
the form 

4 + s* - 6*) ~ ((i^ - b'){b^ - c^)t/ = 0. 

Ifrom tliis equation it appears that if the point of intersection 
of tlie three surfaces 

V + bY + cV - V = 0, ir* + //* + 2* - 6* = 0, // = 0 

be taken as origin, the lowest terms in the equation of the 
wave-surface are of the second degree, and therefore that the 
origin is a double point on the wave-surface at which there 
is a tangent cone of the second degree. 
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If we seek for the coordinates of the points of intersection 
of the three surfaces, we have 

wV + cH^ = riV, y = 0 5 

whence we obtain for the coordinates of the point tlie 
expressions 


= c^ 


- c^ 


, = 0 , ^ w 




( 66 ) 


The equation of the circular sections of Fresnel’s ellipsoid 


IS 


— ^ s" — r- 


= 0 ; 


whence, if ot,, ^ 2 , and wg denote the direction-cosines of the 
perpendicular to a plane of circular section, we have 


OTi* = 




= 0 , 


2 = ~ 


(67) 


From (66) and (67) it appears that a singular ])oint on 
the wave-surface is on a perpendicular to the plane of a 
circular section of Fresnel’s ellipsoid at a distance b from 
the origin. 

The existence of sucli points follows readily from the mode 
of generation of the wave-surface described in Art. 260. From 
thence it appears tliat the perpendicular to each section of 
Fresnel’s ellipsoid meets the wave-surface in two points whose 
distances from tlie centre are equal to the principal semi-axes 
of the section. 

If the section bo circular, every axis is a principal axis, 
and all the corresponding points on the wave-surface coalesce 
into one. 

The perpendiculars on the corresponding tangent-plan gj 
of tlie ellipsoid are, however, not in the same plane ; and thus 
corresponding to the one ray going from the centre to the 
singular point tliere are an infinite number of wave- fronts, 
that is, an infinite number of tangent-planes to tlie wave- 
surface meeting at the singular point. 
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As tlie wave-normals and velocities of propagation are 
different for these fronts, when the ray reaches tlie surface 
of the crystal it is refracted into an infinite number oi rays, 
forming a cone, and the phenomenon exhibited is termed 
wn ica I 1 'pfi 'action . 

From the consideration of the ellipsoid reciprocal to 
FresneFs ellipsoid, it is easy to see that the wave-surface 
must possess singularities of another kind in addition to 
those mentioned above. 

From Art. 260, it appears that the perpendicular to each 
section of the reciprocal ellipsoid is perpendicular to two 
tangent-planes of the wave-surface, and meets them in points 
whose distances from the centre are the reciprocals of the 
semi-axes of the section. If the section bo a circular section, 
every axis is a principal axis, and all the corresponding feet 
of perpendiculars on tangent-planes to the wave-surface 
coalesce into one. 

Tlie centnjl radii of the reciprocal ellipsoid are co-direc- 
tional witli perpendiculars on tangent-planes of Fresnel's 
ellipsoid, which are the* reciprocals of the radii, so tliat all 
the perpendiculars to tangent-planes of Fresnel’s ellipsoid 
which lie in a circular section of the reciprocal ellipsoid are 
equal to the mean semi-axis of Fresnel’s ellipsoid, and corre- 
spond to a single tangent-plane to the wave-surface. The 
corresponding radii of Fresnel's ellipsoid do not, liowever, lie 
in tlic same plane, and are not equal, so that tliere are an 
infinite number of rays corresponding to the same wave-front 
which must therefore touch the wave-surface all along a curve. 
To find the nature of this curve, we may proceed thus. 

Let 7^ denote the length of the central perpendicular on a 
tangent-plane of Fresnel’s ellipsoid, and a, j3, y its direction- 
angles. 

If /? lie in the circular section of the reciprocal ellipsoid, 
we liave p = b, and therefore 

u® cos® a + 6® cos®/3 + o®cos®7 = 6® (cos® a + cos®j3 + cos® 7) ; 
tliat is, (a* - 6®) cos®a - (6® - c®) cos®7 = 0. 

Also, cos® a + cos® 7 = siu®/3 ; 

b^ — c® — ^® 

wlieiice cos®a = sin®/3, cos®7 = j sin®Jl3 (68) 

0 / 0 a 0 
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Let r denote the central radius of FresneUs ellipsoid to the 
point of contact of the tangent-plane perpendicular to 2 ^, then 

cosV/ + ¥ cos*j3 + CO.'S’ y 
^ " jf 

and if p denotes the distance of this point of contact from 
the foot of the perpendicular, p* = r* - //. In the present 
case, 2 ^ - b, and we have 

COS’ a + b* cos^ 13 + 6-^ y - 

p = ^ 

(a^ - b^) COS’ a - (b^ - r}) cos* y 

= p 


Substituting for cos*a and cos^y their values from (68), we get 


p' = 


(a^ - />*) (i* - c*) 
b^ 


sin 


(3. 


It is plain, from the construction in Art. 260, that p is 
the distance from the foot of the perpendicular on the tangent- 
plane to the wave-surface to its point of contact, and that this 
distance is parallel to the corresponding direction of displace- 
ment in the wave-plane. In the present case the wave-plane 
contains the axis of y, and (3 is the angle which the electric 
displacement makes with this axis. Hence (3 is the angle 
which th© line from the foot of the perpendicular to the 
point of contact of the wave-front with the wave-surface 
makes with a parallel to the axis of y in the wave-front*. 

Accordingly, 


P 


~ b^) (6* “ C*) 


sin j3f 


(69) 


is the equation of tlie curve along which the wave-front 
touches the wave-surface. This curve is therefore a circle 
whicii touches the parallel to the axis of y at the foot of the 
perpendicular from the centre, and whose diameter is denoted 
by the expression 

b 
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Corresponding to tL ^ ^ *u.iie we have been considering, 

there are an infinite number of rays which meet the wave- 
front along its circle of contact with the wave-surface. All 
these rays have the same wave-normal, and are propagated 
with the same normal velocity. Hence, when they are re- 
fracted at the surface of the crystal, the emergent rays are 
parallel and form a cylinder. Unless tlie wave-normal be 
normal to the surface, the section of this cylinder made by 
the plane bounding the crystal is an ellipse. 

The remarkable phenomena described above were foretold 
by Hamilton as consequences of ^properties of Fresiiers wave- 
surface discovered by him. They were realized experimentally, 
first by Lloyd, and long afterwards by Fitzgerald. 

278. Total ReHexion. — When light passes from a denser 

into a rarer medium, if the angle of incidence exceed sin"^ -y 

where \x denotes the relative index of refraction of the media^ 
there is no refracted ray. In fact, under these circumstances^ 
a refracted wave-plane is impossible, as it would in the case 
of an isotropic medium, *be a tangent-plane to a sphere drawn 
through a line lying inside the sphere. If both media be 
isotropic, equations (43) seem impossible to satisfy; for, if 
we suppose /i zero, these equations cannot be satisfied unless 
we make I and V each zero. 

Mathematically it is possible to give a solution of equa- 
tions (43), wliich in its final result is physically satisfactory ; 
but it seems impossible to obtain a sound physical basis for 
the equations themselves. 

The mathematical solution is as follows : — Assume 


where 


D = ae"<^, D' = 

£ = 1 , and (j, = 


I)i = ai 

- {/x + my + ?/s)j, 


{Vt - {I'x ^ m'y ^ n'z)], 
A 


= Y- {Vif- (A* + wi// + «i2)l ; 

Ai 

then the differential equations of ^Yave- propagation are 
satisfied, and D, &o., are periodic. 
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If we now suppose the incideat light polarized in the 
•])laiie of iiioidenoe, since 

r, 1 r 1 

Ai Ti At 

ut the origin, where jr, //, and % are all zero, we have 

0 =r 0' = ; 

and as equations (43) mathematically hold good, we have 
sin ix cosi - cos ii sin \ 


/' = V 




sin ii cos ^ + cos ti sm i 
Blit sin ii = fx sin i, cos ii = i sin^/ - 1), 
and therefore 

T fi sin i cos i - i sin i a/(m^ sin* - 1) _ 1 - i ta n € 

I u sin i cos i + i mi i aiiiU - 1) 1 + / tan £ 

v/(ju* sin*/~l) 


= (cos € + L sin f)"* = where’ tan e = 


/X cos i 


Hence 


jD' = s= a^*d<^'+2e), 


and, accordingly, the intensity of the reflected light is equal 
to tliat of the incident; but its phase is increased by ^e. 
Again, 

27r 

fpi ^ {Vit - {hx + nhy + «is) } ; 

Ai 

and since the axis of x is normal to the surface separating 
tlie media, and the axis of z perpendicular to the plane of 
incidence, we have 

/i = cos ix = I \/(/i* sin*e - 1), fUi = fx sin e, tii = 0, 


2iir 


2it, 




1) 2iir ^ 

ai e M . 


( 70 ) 
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In this expression, the power of e whose index is real is 
not periodic ; and since Ai is very small when m is of sensible 
magnitude, this factor tends to become very small. Hence 
at any sensible distance from the boundary 2>i is very small, 
and there is no visible refracted ray. 

In other cases of total reflexion a similar mode of treat- 
ment may be employed. The results obtained above satisfy 
the mathematical conditions holding good when the reflexion 
is not total, and the final result is consistent with the observed 
phenomena ; but the whole investigation can scarcely be 
regarded as having any physical validity. 

279. Absorption of — When a medium is not a 

perfect insulator, an electromotive force produces not only a 
change of electric displacement but also a conduction-current. 

If C be the electric conductibility of the medium, the 
resistance of an element of unit section parallel to the axis 
duo 

of is , ajad the electromotive force for this element is 

Xdoo. Hence, if A dewte the intensity of the conduction- 
current parallel to ir, we have ii = CX. 

The total current is made up of the conduction-current 
and that due to a change of the electric displacement ; 
accordingly, we have 

u = +/= CX -f/, 

and as X = we obtain =/+ 


Substituting for u in terms of the components of magnetic 
force, and for the latter in terms of those of displacement, by 
means of equations (13), (15), and (11), we get 


/ + 


47rg 

K 




The last term in this equation is zero; and if we take 
the normal to the plane of the wave as the axis of z, the 
displacement /is a function of z only, and (71) becomes 


^ ^ df ^ 1 d^f ^ 
dt^ ^ K dt tjsK dz^ 


(72) 


o 
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If U denote the velocity of wave-propagation when there 
is no absorption, we have 

— = TP 

nK ’ 

and putting ttwC? = k, we get 


= UTP. 


Thus (72) becomes 


ilt' <U <lt' 

To solve tins equation, we may assume 


where 


= y-i, 


and m is a quantity to be determined so as to satisfy (73). 
We have, then, 

- n* + 4/c Wni = - m’ ; 


that is. 


- 


Assume m = q - ip, then 


pq=^2kn. 


Eliminating we get 


whence 







Absorption of Light. 


195 


As p is real, must be positive, and therefore 


in' IV »■ y j 


Here h is of the same order as 0, which is of the order 
IC 1 

QV — • 

^ iPd^ 

4*A*i7* . T* 

Hence — ; — is of the form v-^, where v is a numerical 

coefficient depending on C and on the units selected, r the 
time of vibration in the wave of light, and 9' the unit of time. 
In order that C should have any sensible magnitude, 9 must 

be enormously great compared witli r. Hence — — — is a 

small quantity, whose square may be neglected in the 
expansion of the squale root, and we have 




= 4k^U\ 


Substituting q - ip for m, we get 

/ = ae'*”" (76) 

As the wave is advancing in the direction of s positive, 
q is positive ; and since pq is positive, p must be positive. 
Hence 

p ‘ZirzsCU, (76) 

also 

2kn n ‘2 IT 


/ = a e 


(77) 
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It follows, from the expression obtained for f that the 
velocity of wave-propagation is Z7, and is therefore unaltered 
by absorption. In consequence of tlie factor the amplitude 
of/ diminishes as % increases. Since p varies as C, unless C 
be very small, the amplitude of f diminishes rapidly, and the 
medium is practically opaque. 

380. Electrostatic and Electromagnetic Measure. — 

The reader of the foregoing pages may have been struck by 
an apparent inconsistency between the present Chapter and 
Chapter XI. 

In Chapter XL the specific inductive capacity k is of the 
nature of a numerical quantity. In tlie present Chapter, the 
specific inductive capacity K is regarded as the reciprocal of 
the square of a velocity. T'he apparent inconsistency results 
from the fact that in Chapter XI. the various quantities are 
supposed to be expressed in electrostatic measure, whereas in 
tlie present Chapter they are supposed to be expressed in 
electromagnetic. 

We must consider the hypotheses on which tlie two modes 
of measurement are based, and how it is that in reference to 
space, time, and mechanical force, the expression for the same 
physical quality of a body is in one mode of expression a 
quantity of a nature different from what it is in the other. 

Lot e and E denote quantities of electricity expressed in 
electrostatic and electromagnetic measure, X and 3£ the corre- 
sponding electromotive intensities, and / and f the displace- 
ments. Let L denote a linear magnitude, and T a portion 
of time, and let us use the symbol = to mean that two 
expressions denote quantities of the same nature. 

Electrostatic measure is based on the assumption that the 
product of two quantities of electricity divided by the square 

of a line denotes a mechanical force, that is, — = mechanical 
force. ^ 

Electromagnetic measure is based on the assumption that 
the product of the strengths of two magnetic poles divided 
by the square of a line denotes a mechanical force, that is, 
if m denote the strength of a magnetic pole, 

« mechanical force « — ; whence e ^ m. 
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Again, i denoting the strength of a current, E= Ti\ but 
i =y, where / denotes the strength of a magnetic shell, and 
jD = magnetic moment = niL ; whence 


E ^ T 


m 

L 



The electromotive intensity multiplied by a quantity 
of electricity denotes in either system of measurement a 
mechanical force ; accordingly, eX = EH ; but 



and therefore 



Again, iD = E^ and fD = e ; whence 


and therefore* 


f = jf, also = X, 


7 - 1 , in 1 


where V denotes a velocity. Thus A; is a numerical quantity, 
but K the reciprocal of the square of a velocity. 

The magnitude of the unit of electricity differs very much 
in tlie two systems of measurement. 

In the electromagnetic system, two units at the unit of 
distance apart act on eacli other with the unit force. 

Ill the electromagnetic system, two magnet-poles of unit 
strength, at the unit distance apart, act on each other with 
the unit force. 

A circular current of unit strength acts on a unit magnet- 
pole at its centre with a force which is 27r times the unit of 
force, provided the radius of the circle be of unit length. 

The quantity of electricity which passes through a section 
of this circuit in the unit of time is the unit quantity of 
electricity expressed in electromagnetic measure. 

The quantity of electricity contained in the electro- 
magnetic unit is n times the quantity contained in the 
electrostatic. 
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If, then, E and e denote the same absolute quantity of 
electricity expressed, one in electromagnetic, the other in 
electrostatic units, and if i and T denote the units of length 
and of time, we have 


but EH = where 36 and X denote the electromotive 
force corresponding to the quantity of electricity denoted 
by E and e ; whence 


Then, 


3E.!^X. 


rp yi rpi 

E = UkX = 

uL uL ir 


whence 


1 ,// 1 
== ft. — ♦« 

X k 


When the second and centimetre are taken as the units 
of time and length, 


// = 3 X 10^® approximately. 
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NOTE ON THOMSON AND DIKICHLET’S THEOEEM, 

Article 70. 

When the number of given surfaces is reduced to one, 
this theorem is proved by Gauss in the following manner : — 

(1) On a given surface 8 a homogeneous distribution of 
a given quantity of mass is always possible, such that J VadS 
is a minimum. For this distribution, V is constant for all 
occupied parts of the surface, and there is no part unoccupied. 

If r denote the longest distance between any two points 
of /S, and M the total mass, it is obvious that at any point 

of 8 the potential cannot be less than — , since the distri- 

T 

bution is homogeneous, that is, composed of mass having 
everywhere the same algebraical sign. Hence J VcFd8 cannot 
AT" 

be less than — • 
r 

Consequently J Vad8 cannot be diminished without limit, 
and tliere must be a distribution such that J V(jd8 cannot be 
made less. In this distribution V must be constant. For, if 
for an occupied portion Si of the surface, Vi be everywhere 
greater than -4, and for another equal portion Sa of the 
surface Fi be everywhere less than Ay at each point of 2i 
let Sa = - V, and at each point of Sa let So- = -f- y, then the 
total mass remains unaltered, and 

S j Vads = J VSadS + / 8 Vad8 = 2 J VSad8 

--2vf(F,-V,)d8y 

since 8F is the distribution resulting from Scr, and therefore, 
by Art. 51, we have 

fSVad8 « SV&ad8. 
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Accordingly, J VcrdS has received a variation which is 
essentially negative, and consequently cannot he a minimum 
for the distribution ts. Hence, when the integral is a mini- 
mum, V is constant for the occupied part of the surface. If 
there were a part unoccupied by Art. 66, V would be less for 
this part than for the occupied part, and hence as before the 
integral could be made less. Accordingly, in the distribution 
for which ^VadS is least, there is no part of the surface 
unoccupied. 

(2) If Z7 be a given function of tlie coordinates, a 
homogeneous distribution of given mass over 8 is possible, 
such that f(V- 2U) adS is a minimum. For this distri- 
bution F- i7 is constant for all occupied portions of the 
surface. 

If Z7' be the largest value of U on 8, it is clear that 
!(V-2U) <Td8 cannot be less than 



and therefore that there must be a distribution such that 
J ( F - 2 C7) ad8 cannot be made less. For this distribution 
F - J7 is constant at all occupied parts of the surface. 

Let (V-2ir)(rd8y then 

81F- J SVad8+j (F- 2U) S<Td8 = 2 J (F- ?7) B<rdS. 

If V- U he greater than A at every point of an occupied 
portion of surface 2i, and less than A at every point of an 
equal portion ^2 of surface, as in (1), SfF can be made 
negative, and therefore W cannot be the least possible. 

In this case, if part of the surface 8 be unoccupied, 
F- Z7 may be greater on this part than it is on the occu- 
pied part, and therefore in this case we cannot show that the 
whole surface must be occupied. 

(3) Suppose now three distributions of mass on 8. 

1. A distribution whose surface-density is (Tq and potential 
Fo, such that / V<rd8 is the least possible, the total mass 
being M. 
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2. A distribution whose surface-density is ai and potential 
Fi, such that J {V - 2c ?7) ad8 is the least possible, the 
total mass being and c being a given constant. 


3. A distribution whose density is a 2 , and potential F„ 
such that 


(T2 — 


(Ti - (To 


£ 


«an(l Fg 



then the total mass is zero, and 


Vz- U- 


F-CZ7- Fo 
£ 


but this is constant for all parts of the surface occupied by <ti. 

If £ be diminished without limit, the distribution cri 
passes into <to, and in this case tliere is no finite portion 
of the surface S left unoccupied. 

H ence, ^hen c is diminished without limit, Vz - 27 is 
constant for the whole surface 8, 

Let us nt)w superpose on <t 2 the distribution whose density 
is ao-o, where a is constant. Then 

F= a Fo H F, and F- Z7= « Fo + F - U. 

By a proper determination of a the right-hand member of 
this equation can be made zero at all points of 8, 

Accordingly for a single surface Thomson and Dirichlet's 
Theorem is proved. 

This theorem in its most general form can readily be 
deduced from the properties of fluid motion. 

Suppose that the given surfaces 8i, /Sg, &c., are sur- 
rounded by liquid, or incompressible fluid, of unit density, 
extending to infinity. Apply to the liquid at each surface 
an impulsive pressure which at each surface is equal to the 
given value of Thomson's function for that surface. The 
liquid begins to move irrotationally, and the velocity potential 
of the motion is the same as the impulsive pressure, and is 
equal at each surface to the given value of Thomson’s 
function, and satisfies Laplace's equation for the whole of 
space. 

r 
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Thus the truth of Thomson and Dirichlet's Theorem is 
established. 

It is easy to show from Greenes theorem that, if there bo 
a given quantity of mass on each of a number of surfaces, 
this mass may be so distributed that the potential shall be 
constant over each surface. 

Let 2W^^\VadS, Then 


87rir=-SJ 



Hence W is essentially positive, and cannot therefore be 
diminished without limit, and there must be a distribution of 
mass such that W cannot be made less. For this distribution 
V must be constant for each surface. For if V be not con- 
stant, W may be made to receive a variation which is negative 
by transferring positive mass from points on the surface where 
Fis greater than A to points where it is less than A, In this 
case it is not necessary that the part of tlie surface from which 
the transfer is made should be occupied. On the other liand, 
if V be constant for each surface, any change of distribution 
increases W, For, let v be the change in Vy then 


SttA TF = - 2 


dV 


dS- 2 


V do 
dv 


where the surface-integrals are to be taken over both sides 
of the surfaces ; but as Fis constant for each surface, and 

~ dS zero. 


since the total mass is constant, we have 


SttA IF * 



since = 0 throughout the field. Hence the change in W 
is essentially positive, and JFis least when Fis constant over 
each surface. 



Oaiiss’s Method. 


203 


The property of the potential made use of in (1) to show 
that the whole surface must be occupied is perhaps more 
readily established by the method of Gauss than by that 
employed in Art. 66. 

Gauss’s method is as follows : — 

If there be no mass outside a surface S on which the 
potential is everywhere positive, its value at a point O, 
outside S, is positive, and less than A its greatest value on S. 

For if the potential P at 0 be greater than draw 
lines in all directions from O; they meet the surface S or go 
to infinity, and the potential on any one of them passes from 
P to , or to some value less than A. Hence on every line 
there is a point at which the value of the potential is P, 
lying between P and A, All these points form a closed 
surface at which tlie potential is constfiiit ; and as there is no 
mass inside it, the potential has the same value tlirougliout 
the interior of the surface, Art. 61. Hence the value of the 
potential at D is P, and cannot consequently be P as was 
supposed. If P were* negative, we could show in like 
manner that the potential at O must lie between P and zero, 
and could not therefore have the supposed value. 

Again, the potential at O cannot bo A or zero. For if it 
has eithex" ot these values, describe a sphei’e round O as centre. 
At no point of the surface of this sphere can the potential 
be greater than A or less than zero. Hence its mean value 
on this surface cannot bo A or zero, unless it have this value 
for the whole surface of the sphere, in which case it would 
have the same value for the whole of space external to S, 
which is impossible. 

If the potential be everywhere negative on the surface S, 
its value at a point O outside S is negative, and less in 
absolute magnitude than its greatest negative value on S. 
This is proved in a manner pi'ecisely similar to that adopted 
in the case of the positive potential. 
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system, energy due to, 90, 
Magnetism, induced, 102. 
distribution of, 103. 
terrestrial, 107. 



Tonler, 
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Magnetization, intensity of, 81. 
Magnetized ellipsoid, 84. 

body, potential of, 82, 80. 
sphere, 83. 

Magnets, constitution and action of, 
81. 

energy duo to, 89. 

Maxwell’s method of forming solid 
harmonics, 28. 

theory of electric displacement, 
124. 

light, 117. 

Medium, dielectric influence of, 123. 
Moment, magnetic, of body, 95. 


Oblate, ellipsoid of revolution, 51, 50. 


Poisson’s equation for magnetiztd 
body, 83. 

Polarization, electric, 121. 

PolaiL/ed light, 108 . 

Potential, oxpanpion of, in solid har- 
monics, 1. 

of surface •distribution *011 pro- 
late ellipsoid, 50. 
on oblate ellipsoid, 56. 

Prolate, ellipsoid of revolution, 39, 
50. 

Propagation of disturbance in electro- 
magnetic field, 150. 


Reflexion and refraction of light, 165. 
from crystalline medium, 176. 
at interior surface of crystal, 184. 
in isotropic media, 172. 
quantities to be determined in, 
172. 

for uniaxal crystal, 180. 
Reflexion, total, 191. 

Revolving fluid, 10. 


Solenoids, 112. 

Solid Harmonics, 1, 26. 

Spherical harmonics, 2, 3, 18, 23, 24, 
34, 38. 

applications of, 7. 

Spheroid, potential of homogeneous, 8. 

heterogeneous, 8. 

Stokes, 13, 91. 

Surface distribution, potential of, 7, 
50, 66, 71. 

Sylvester, 29. 


Tesseral harmonics, 20. 

Thomson, J. J., 124, 125. 

Thomson and Dirichlet’s theorem, 
199. 

Total reflexion, 191. 


Uniaxal crystals, 178. 

Fresnel’s ellipsoid for, 178. 

wave- surface for, 178. 
reflexion and refraction for, 180, 
Unirad ial directions, 178. 


Vector, potential, of magnetic induc- 
tion, 90, 92. 
for magnetic particle, 95. 


Wave, 152. 

-front, 153. 

-length, 152. 

-plane, 153. 

-surface, 157. 

in crystalline medium, 157. 

equation of, 160, 178. 

singularities of, 187. 


THE END. 






